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Chapter  1:  Introduction 


1.1.  Motivation:  Chemical  Thrusters  for  Micro  satellites 

Advances  in  miniaturizing  electronics  and  hence  miero-avionies  has  resulted  in  the 
building  of  nano/pieo  satellites  (satellites  of  the  mass  range  of  0.1  to  10  Kg).  Smaller 
satellites  are  appealing  for  several  reasons.  Most  of  the  eontent  presented  in  §1.1  ean  be 
found  in  Mieei’s  book  [1],  and  is  summarized  and  itemized  here.  First,  a  eonstellation  of 
smaller  satellites  is  less  subjeet  to  single  point  failures  and  henee  ean  be  more  reliable 
than  a  single  eonventional  size  satellite.  Seeond,  having  smaller  satellites  reduees  the 
payload  weight  and  henee  the  overall  mission  eost.  It  is  estimated  that  about  30%  of  a 
spaee  mission  eost  is  usually  from  the  launeh  eosts  whieh  in  turn  depends  on  the  mass  of 
the  spaeeeraft.  Some  examples  [1]  of  mierospaceeraft  are  MightySat-1  (U.S.  Airforee), 
SNAP-1  (Surrey  Spaee  Centre,  England),  and  Faleonsat  (U.S.  Air  Foree  Aeademy). 

Sueh  small  satellites,  however,  require  small  propulsion  systems.  There  are  a  variety 
of  options  whieh  inelude  ehemieal  thrusters,  eleetro-thermal  thrusters,  eold  gas  thrusters. 
Hall  thrusters,  eolloid  thrusters,  other  Eleetrie  Propulsion  ete.  Eaeh  type  has  its  merits  and 
demerits  and  the  optimum  thruster  ehoiee  ultimately  depends  on  the  partieular  mission. 
As  will  be  shown  towards  the  end  of  §1.1,  ehemieal  thrusters  are  espeeially  desirable  for 
missions  requiring  a  high  thrust  and  relatively  low  propulsion  system  weight. 

l.l.l.  Size  Classes  of  Satellites 

Mierospaeeeraft  ean  be  further  subdivided  into  elasses  based  on  their  mass,  power 
eonsumption,  dimensions,  ete...  Table  1-1  provides  a  eommonly  adopted  sub 
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classification  for  small  scale  satellites.  Generally  speaking,  the  number  of  propulsion 
systems  from  which  to  choose  deereases  with  the  size  of  the  microsatellite. 


Designation 

Mass,  kg 

Power,  W 

Dimension,  m 

Microspacecraft  (A/F  European  definition) 

10-100 

10-100 

0.3-1 

Class  I  microspacecraft  (nanosat) 

5-20 

5-20 

0.2 -0.4 

Class  II  microspacecraft 

1-5 

1-5 

0.1  -0.2 

Class  III  microspacecraft  (picosat) 

<1 

<1 

<0.1 

Table  1-1:  Definition  and  classification  of  microspacecraft  (Adapted  from 
Mueller  [1]) 

1.1 .2.  Missions  for  Microsatellites 

Typical  missions  for  microsatellites  include: 

1)  Moving  from  a  central  distribution  point  to  set  up  a  large-baseline  interferometer 
(satellites  separated  by  tens  or  hundreds  of  km) 

2)  Stationkeeping  within  the  large -baseline  interferometer 

3)  Intercepting  another  satellite  from  a  parking  orbit 

4)  Stationkeeping  while  in  the  parking  orbit 

5)  Attitude  control 

The  next  subsection  shows  that  the  thrust  requirements  are  different  for  each  of  these 
missions.  The  mission  of  interest  in  this  dissertation  is  a  rendezvous/intercept  mission 
which,  as  we  will  see,  requires  a  high  thrust/weight  propulsion  system. 
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1.1 .3.  Thruster  Performance  Requirements 

Interferometry  missions  are  defined  by  their  tight  attitude  eontrol  requirements.  For 
example,  the  Laser  Interferometer  Spaee  Antenna  (LISA)  mission  eonsists  of  three  small 
satellites  orbiting  the  earth  in  a  earefully  designed  formation.  The  thrust  requirements  for 
this  mission  are  in  the  range  of  2  -  20  pN  with  an  allowable  varianee  of  0.1  pN.  These 
numbers  are  typieal  of  formation  type  flying  missions.  Station  keeping  missions  within 
eonstellations  require  a  higher  thrust  of  the  order  of  millinewtons.  Rendezvous  and 
intereept  missions  (whieh  is  the  foeus  here)  are  expeeted  to  have  the  highest  thrust 
requirements  with  AV’s  of  the  order  of  several  thousands  of  meters  per  second.  In 
addition  to  the  thrust  requirement,  intereept  missions  also  require  high  thrust  to  weight  to 
enable  higher  maneuverability  owing  to  the  eomplexity  of  the  orbital  transfers  required 
for  intereept.  Enabling  the  development  of  ehemieal  thrusters  with  sueh  high  thrust  to 
weigh  ratio  requirement  is  the  primary  foeus  of  this  dissertation. 

1.1.4.  Available  Propulsion  Systems 

Microthrusters  can  be  subdivided  into  a  few  broad  categories:  chemical  thrusters 
(reactive  and  cold  gas),  electric  propulsion  based,  and  emerging  technologies  based 
(MEMS  type).  Of  these,  chemical  thrusters  are  most  common  and  include: 
i)  Bipropellant  (liquid  rocket)  engines:  Fuel  and  oxidizer  stored  separately  are  burned  in 
a  combustion  chamber  to  produce  hot  gases  that  are  expanded  out  through  a  nozzle  to 
produce  thrust.  These  are  typically  used  when  high  thrust  is  required. 
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ii)  Monopropellant  (solid  rocket)  thrusters:  Similar  in  concept  to  bi-propellant  engines 
except  that  the  propellant  (usually  a  solid)  contains  both  fuel  and  oxidizer.  These  are 
often  used  in  initial  parts  of  launch  sequences. 

iii)  Cold  gas  thrusters:  Thrust  is  produced  by  expanding  a  stored  cold  gas  through  a 
nozzle.  These  are  usually  used  in  low  thrust  requirement  missions  (orientation 
control). 

iv)  Hybrid  rocket  motors:  Liquid  or  gaseous  oxidizer  is  injected  into  the  combustion 
chamber  which  contains  a  solid  fuel.  The  hot  gases  are  expanded  through  a  nozzle  to 
produce  thrust. 

Examples  of  electric  thrusters  include: 

i)  Ion  engines:  The  propellant  (typically  a  noble  gas)  is  ionized  in  a  plasma  chamber  and 
accelerated  through  a  charged  grid  to  produce  thrust. 

ii)  Hall  thrusters:  Similar  to  the  ion  engine  except  that  thrust  production  is  enhanced  by 
applying  a  magnetic  field.  This  generates  a  Lorentzian  force  on  the  ion  stream 
causing  the  ions  to  bombard  more  incoming  propellant  and  increasing  ion  density. 

iii)  Field  Emission  Electric  Propulsion  (FEEP):  Conceptually  the  same  as  Ion  and  Hall 
thrusters  and  differs  only  in  the  production  of  the  ions  which  occurs  via  field  emission 
from  a  liquid  metal  surface  with  the  metal  being  the  propellant. 

iv)  Pulsed  plasma  Thrusters  (PPTs):  In  PPTs,  the  propellant  is  first  ionized  and  then 
electromagnetically  accelerated  between  two  parallel  electrodes  in  pulses. 

v)  Resistojets:  A  liquid  or  gaseous  propellant  is  heated  by  a  heater  element  and  the  hot 
vaporized  gas  is  expanded  through  a  nozzle  to  achieve  thrust. 
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vi)  Colloid  Thrusters:  An  ionic  liquid  is  and  accelerated  electrostatically  to  produce 
thrust.  If  the  ionic  liquid  is  HAN,  a  colloid  thruster  can  be  paired  with  a  chemical 
thruster  to  make  a  hybrid  thruster  system  able  to  produce  high  thrust  for  repositioning 
and  small  thrust  (high  Isp)  for  stationkeeping  [2], 


Type 

Isp?  S 

Total  system  efficiency 

% 

Mass, 

kg 

Thrust, 

mN 

Cold-gas  Thruster 

75 

95+ 

0.01 

5 

Solid  Rocket  Motor 

185 

90+ 

1.6 

I(P 

Digital  MEMS 

200 

75 

0.04 

50 

Monopropellant  (N2H4) 

220 

95+ 

0.16 

100 

Bipropellant 

290 

- 

0.55 

22000 

PPT 

800 

10 

0.40 

500 

Hall  Thruster 

1500 

35 

1 

10 

FEEP 

8000 

25 

305 

800 

Table  1-2:  Performance  comparison  of  different  types  of  microthrusters  (Adapted 
from  J.H,  Schilling  [1]) 

Other  emerging  technologies  in  the  Micro  Electro  Mechanical  Systems  (MEMS)  area 
include  combining  MEMS  concepts  with  one  of  the  above  technologies,  like  for  example 
a  MEMS-based  FEEP  or  MEMS-based  cold  gas  thrusters.  MEMS-based  technologies 
involve  using  various  chemical  etching  processes  to  form  miniature  versions  -  usually  in 
Silicon  -  of  the  various  thrusters  discussed  above.  One  simple  example  are  digital 
thrusters  which  consist  of  tiny  wells  full  of  monopropellant.  Each  well  can  be  ignited 
individually  to  produce  a  small  impulse  bit.  It  is  beyond  the  scope  of  this  introduction  to 
delve  exhaustively  into  each  of  these  technologies  but  a  good  review  of  these  can  be 
found  in  [1]. 
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1.1 .5.  Advantages  of  Chemical  Propulsion 


Intercept  missions  require  a  high  AV  and  more  specifically  a  high  thrust  to  weight 
(T/W)  thruster.  From  table  1-2,  it  is  clear  that  chemical  thrusters  have  the  highest  T/W  (of 
the  order  of  4  and  above)  and  therefore  would  be  the  best  choice  for  these  missions. 

1.2.  Chemical  Rockets 

In  the  previous  section  we  concluded  that  chemical  propulsion  has  the  advantage  of 
high  thrust  to  weight  ratio  compared  to  the  other  propulsion  options.  In  this  section,  the 
thrust  production  mechanism  in  chemical  rockets  is  explained  briefly  for  solid  and  liquid 
rocket  motors  along  with  the  relevant  engine  cycles  etc.  Consideration  is  also  given  to  the 
performance  metrics  that  characterize  chemical  rockets. 

1.2.1.  Thrust  Generation 

Fig  1-1  is  a  schematic  illustration  of  a  chemical  rocket.  It  consists  of  a  combustion 
chamber  coupled  to  a  convergent-divergent  (c-d)  nozzle.  The  propellant  is  burned  in  the 
combustion  chamber  to  produce  a  high  temperature  gas  that  is  accelerated  as  it  passes 
through  the  c-d  nozzle  to  exit  the  nozzle  at  a  high  velocity  (f/g).  The  pressure  is  highest  in 
the  combustion  chamber  and  decreases  as  the  flow  accelerates  through  the  nozzle. 
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Figure  1-1:  Schematic  sketch  of  a  thrust  chamber  for  a  chemical  rocket  (adapted 
from  [3]) 

The  thrust  can  be  calculated  by  performing  an  integration  of  forcers  over  the  control 
volume  illustrated  by  the  blue  dashed  lines  in  the  figure.  The  final  expression  is  given 
below: 

T'h  =^^e+(^e-^oo)^e  (I'l) 

Th  is  the  total  thrust  produced  by  the  thruster,  m  is  the  mass  flow  rate  through  the 
thruster,  pe  is  the  exit  pressure  at  the  nozzle  exit,  Ae  is  the  exit  area  and  is  the  ambient 
pressure  (in  vacuum  this  would  be  close  to  zero).  The  equation  above  shows  that  for  a 
perfectly  expanded  nozzle  {pe  =  pcc)  the  thrust  is  a  function  of  only  the  exit  velocity  and 
the  mass  flow  rate. 

1.2.2.  Types  of  propellants 

Rocket  engines  differ  from  engines  used  in  other  vehicles  (like  cars  or  airplanes)  in 
that  they  carry  both  the  fuel  and  oxidizer  on  board.  This  is  because  they  often  operate  in 
space  where  oxygen  is  not  available.  Propellants  can  be  broadly  categorized  into  two 
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types:  monopropellants  and  bipropellants.  Monopropellants  have  both  the  fuel  and 
oxidizer  mixed  together  (as  is  the  ease  with  solid  propellants  and  some  liquid 
propellants).  Bipropellants  usually  eonsist  of  a  fuel  and  an  oxidizer  stored  separately  and 
mixed  together  in  the  eombustion  ehamber.  In  this  seetion  we  make  a  slightly  different 
distinction  and  categorize  propellants  as  either  solid  or  liquid. 

i)  Solid  Propellants:  Solid  propellants  can  consist  of  two  or  more  of  the  following  [ref 
Sutton]:  Oxidizer  (Nitrates  or  perchlorates),  Fuel  (plastics  or  metal  hydrides),  Binders 
(compounds  that  combine  fuel  and  oxidizer  together  in  a  matrix).  Additives  (for 
controlling  burning  rate),  and  Inhibitors  (to  restrict  the  burning  surface).  Typically,  solid 
rocket  motors  form  the  initial  stages  of  the  launch  where  high  thrust  is  required.  A 
drawback  of  using  solid  propellants  is  that  they  cannot  be  stopped  once  they  are  ignited. 
Also,  the  control  of  the  combustion  rate  relies  on  inhibitors  and  other  additives  to  the 
propellant  mixture.  Therefore,  combustion  control  is  generally  not  possible. 

Composite  solid  propellants  consist  of  separate  fuel  and  oxidizer  neither  of  which 
would  bum  satisfactorily  independently.  The  oxidizer  is  usually  ground  into  a  small 
particles  and  dispersed  in  a  fuel  matrix.  Homogenous  propellants  are  capable  of 
combustion  in  seclusion.  These  are  often  unstable  chemicals  such  as  nitroglycerin  in 
which  the  oxidizer  occurs  bonded  to  the  fuel.  Some  typical  solid  propellants  and  their 
properties  are  listed  in  table  1-3  below. 
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Oxidizer 

Fuel 

P 

(g/cm^) 

^comb. 

(K) 

c 

(m/s) 

Mwt 

(kg/mol) 

kp 

(s) 

Ammonium 

Nitrate 

11%  binder  and  7% 

additives 

1.51 

1282 

1209 

20.1 

192 

Ammonium 

perchlorate 

78-66% 

polymer  binder  and  4- 

20%  aluminum 

1.69 

2816 

1590 

25.0 

262 

Ammonium 

perchlorate 

84-88% 

12%  polymer  binder 

and  4-20%  aluminum 

1.74 

3371 

1577 

29.3 

266 

Table  1-3:  Theoretical  performance  of  typical  solid  rocket  propellant  combinations 
(Adapted  from  [1]) 

Broadly  speaking  the  manufaeturing  process  for  composite  propellants  involves 
grinding  of  oxidizer  crystals,  mixing  of  appropriate  quantities  of  oxidizer  and  fuel  along 
with  additives  and  inhibitors  (usually  in  a  batch),  casting  of  the  grain  (into  the  desired 
shape)  and  curing  it  usually  by  heating.  For  homogenous  propellants  the  process  involves 
using  a  solvent  in  which  the  propellant  and  additives  may  be  added  and  mixed  together 
and  then  cast  into  the  desired  shape  and  dried  by  evaporation  of  the  solvent. 

ii)  Liquid  Propellants:  Liquid  propellants  can  be  produced  as  mono  or  bipropellants. 
Monopropellant  liquids  have  both  the  oxidizing  agent  and  the  combustible  substance 
together.  An  example  of  this  would  be  hydrogen  peroxide  or  hydrazine.  Table  1-4  gives 
properties  of  some  common  liquid  propellants. 
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Oxidizer 

Fuel 

Avg.  speeifie 

gravity 

Chamber 

T(K) 

Chamber 

e*  (em/s) 

Mwt 

(kg/mol) 

kp 

(s) 

Oxygen 

Methane 

0.81 

3526 

18355 

18.3 

296 

Oxygen 

Hydrogen 

0.26 

2959 

2428 

8.9 

386 

Oxygen 

UDMH 

0.96 

3542 

1835 

19.8 

295 

N2O4 

Hydrazine 

1.20 

3258 

1765 

19.5 

283 

H2O2  (90%) 

RP-1 

1.29 

2760 

21.7 

21.7 

297 

Table  1-4:  Theoretical  performance  of  typical  liquid  rocket  propellant  combinations 
(adapted  from  [3]) 

Monopropellants  are  typieally  stable  at  room  temperature  and  pressure  but 
deeompose  when  heated  or  reaeted  with  a  eatalyst  exothermieally  to  produee  hot  gases. 
Bipropellant  liquids  are  those  in  whieh  the  oxidizer  is  separate  from  the  fuel.  Typieal 
examples  are  liquid  hydrogen/oxygen,  kerosene/oxygen  ete.  Liquid  propellants  are 
typieally  used  in  booster  stages  and  a  few  seeond  stages.  Liquid  propellants  are  more 
attraetive  than  solid  propellants  beeause  the  eombustion  proeess  ean  be  eontrolled  and 
even  stopped  by  varying  the  propellant  flow  rates.  Liquid  propellants  also  have  a  higher 
Isp  than  solid  propellants.  The  desirable  properties  for  liquid  propellants  are  high  speeifie 
gravity  (or  density)  for  eompaet  storing  and  high  kp. 

1.2.3.  Rocket  Cycles 

Solid  and  liquid  roekets  have  vastly  different  eyeles  of  operation.  Solid  roeket  motors 
have  a  simpler  eyele  as  eompared  to  liquid  roeket  motors  (diseussed  next).  Solid  roeket 
motors  usually  eonsist  of  a  ehamber  inside  whieh  the  solid  propellant  is  east  in  a  given 
shape  (star,  eylindrieal  ete),  an  igniter  (eontrolled  eleetronieally)  to  ignite  the  propellant, 
a  eavity  to  aeeumulate  the  hot  gases,  and  a  nozzle  attaehed  to  the  ehamber  through  whieh 
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the  hot  combustion  gases  expand  to  produce  thrust.  The  chamber  is  a  pressure  vessel 
capable  of  withstanding  large  pressures.  A  typical  chamber  pressure  is  -200  bar. 


Figure  1-2:  Sketch  of  solid  rocket  motor:  Lateral  view. 

Solid  rocket  motor  cycle  can  be  broken  into  three  time  sequence  intervals:  1)  Ignition 
time  lag:  This  is  the  time  taken  from  receiving  an  electrical  signal  to  the  igniter  to  the 
actual  commencement  of  ignition,  2)  Flame  spreading  interval:  As  the  name  suggests  this 
is  the  time  taken  for  the  flame  to  spread  throughout  the  burning  surface  till  the  grain  area 
is  starting  to  bum  and  3)  Chamber  filling  interval:  The  time  required  to  fill  the  internal 
cavity  to  reach  initial  equilibrium  chamber  pressure  and  flow  rate.  The  burning  rates  of 
the  solid  propellants  are  usually  characterized  using  the  chamber  pressure  given  by  a 
power  law  as: 

r  =  a  p"  (1-2) 

In  this  expression,  r  is  the  rate  of  burning, is  the  chamber  pressure  and  a  is  an  empirical 
constant  dependent  on  temperature.  As  discussed  previously,  typical  solid  rocket  motors 
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are  single-use.  Onee  the  propellant  is  ignited,  the  eombustion  proeess  does  not  stop  until 
the  propellant  grain  is  eompletely  burned. 

Liquid  roeket  engines  are  more  eomplex  than  solid  roeket  motors  and  require  a 
system  of  pumps  (to  deliver  the  fuel  and  oxidizer  to  the  eombustion  ehamber)  and 
turbines  (to  drive  the  pumps).  These  are  usually  turbo-pumps  that  are  driven  by  the 
exhaust  of  a  separate  gas  generator  whieh  runs  fuel  rieh.  There  are  many  different  roeket 
eyeles  but  they  all  fall  into  two  broad  eategories:  open  and  elosed.  In  the  open  eyele,  the 
turbopump  exhaust  is  either  expanded  separately  in  its  own  separate  nozzle  or  added  to 
the  main  nozzle  in  the  expanding  seetion  after  the  throat.  This  is  a  slightly  less  effieient 
eyele  than  the  closed  cycle  in  which  the  turbopump  exhaust  (which  contains  unburned 
fuel)  is  fed  into  the  combustion  chamber  to  react  completely.  A  few  types  of  open  and 
closed  cycles  are  illustrated  in  Figure  1-3  and  described  below.  Much  more  extensive 
discussions  are  available  elsewhere[4]. 
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Figure  1-3:  Liquid  Rocket  Turbopump  Cycles:  a)  Fuel  rich  staged  combustion 
cycle;  b)  Fuel  and  Oxidizer  rich  with  auxiliary  pumps;  c)  Simple  expander  bleed 
cycle  d)  Expander  cycle  with  internal  heat  exchanger  [4] 

i)  Staged  combustion  cycle:  This  is  a  closed  cycle  and  in  which  a  pre-burner  is  employed 
to  generate  power  to  drive  the  turbines  that  operate  the  fuel  and  oxidizer  pumps.  Fig.  1-3 
a  shows  a  simple  staged  combustion  cycle  involving  a  fuel  rich  pre-burner  which  drives 
two  turbines  which  in  turn  operate  the  fuel  and  oxidizer  pumps.  Fig  1-3  b  shows  a  more 
advanced  version  of  this  cycle  wherein  two  separate  pre-burners  are  employed  (one  fuel 
rich  and  one  oxidizer  rich)  with  an  auxiliary  pump  to  deliver  fuel/oxidizer  to  the  pre- 
burners.  In  both  cases  the  exhaust  from  the  turbines  are  fed  into  the  combustion  chamber 
and  burned  to  extract  all  of  the  remaining  energy  from  it  before  expanding  the  flow 
through  the  main  nozzle. 

ii)  Expander  cycle:  Expander  cycles  are  open  cycles  in  which  the  working  fluid  (usually 
liquid  H2)  is  used  to  cool  the  nozzle  (see  Fig  1-3  c  and  d).  It  is  vaporized  during  this 
process  and  drives  the  turbopump  which,  in  turn,  pumps  the  fuel  and  oxidizer  into  the 
combustion  chamber.  The  exhaust  from  the  turbopump  is  injected  into  the  expander 
portion  of  the  nozzle  to  extract  some  thrust  out  of  it.  However,  part  of  the  energy  in  the 
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working  fluid  remains  untapped  as  unbumed  fuel  is  not  eombusted.  The  advantage  of 
expander  eyeles  is  their  relative  simplieity  and  fewer  eomponents  (no  pre-burner)  whieh 
assist  in  lowering  the  weight.  But,  their  performance  is  poorer  than  staged  combustion 
cycles  as  they  not  extract  the  maximum  energy  from  the  propellants. 

1.2.4.  Principles  of  Operation 

So  far  we  have  seen  some  important  features  of  chemical  rockets  such  as  thrust 
production,  types  of  propellants,  and  operational  cycles  each  in  a  separate  manner.  In  this 
subsection  it  is  attempted  to  present  the  complete  operational  principles  of  a  chemical 
rocket  drawing  upon  the  previous  three  subsections  and  integrating  them.  Let  us  start 
with  a  generic  liquid  rocket,  as  our  motivation  is  to  understand  maneuverable  chemical 
thrusters,  and  sequentially  breakdown  the  entire  operational  procedure. 

i)  Mixing  of  fuel  and  oxidizer:  The  starting  point  in  the  operation  would  be  to  employ 
turbopumps  as  described  in  the  liquid  rocket  cycles  and  bring  fuel  and  oxidizer  from  the 
tanks  in  which  they  are  stored  to  the  combustion  chamber  with  the  desired  flow  rates  to 
achieve  maximum  combustion  efficiency.  Next,  the  fuel  and  oxidizer  need  to  mix  and 
mixing  length  studies  in  the  past  (simple  Burke- Schumann  [5]  for  gas  phase  mixing)  and 
more  recent  liquid  impingement  and  mixing  in  the  presence  of  turbulence  ([6]  and  [7]) 
have  developed  methods  to  achieve  rapid  mixing  to  reach  stoichiometric  proportions  in 
the  shortest  possible  length. 

ii)  Ignition  and  flame  propagation:  Once  the  fuel  and  oxidizer  have  mixed  the  next 
step  is  to  ignite  them  and  achieve  a  stable  combustion  front  to  burn  them  and  extract  the 
energy  stored  in  the  chemical  bonds  of  the  fuel  and  oxidizer  to  produce  hot  combustion 
gases.  Ignition  studies  have  been  undertaken  ([8]  and  [9]),  to  provide  a  detailed 
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understanding  of  the  ignition  proeess.  Combustion  studies  (diseussed  in  literature  review) 
foeus  more  on  the  understanding  of  flame  propagation  onee  ignition  has  been  aehieved. 
In  this  dissertation  the  emphasis  will  be  placed  on  understanding  the  processes 
responsible  for  enhancing  the  burning  rate  and  flame  stabilization  in  microscale 
combustors  such  as  might  be  used  in  a  small  scale  chemical  thruster 

hi)  Flow  acceleration  and  Thrust:  As  seen  in  the  earlier  subsection  the  main  thrust 
producing  mechanism  in  chemical  rockets,  which  is  the  next  step  once  combustion 
produces  hot  gases,  is  to  expand  the  flow  through  a  c-d  nozzle.  Micronozzle  studies 
(discussed  in  literature  review)  analyze  the  performance  of  nozzles  as  the  scale  of  the 
device  is  reduced.  Some  important  considerations  involve  3-D  effects,  boundary  layer 
growth  and  heat  transfer  to  the  nozzles  structure  from  the  hot  gases. 

iv)  Cooling:  Finally  it  is  required  to  have  a  reliable  way  of  cooling  the  nozzles 
structure  to  prevent  structural  disintegration  owing  to  excessive  heating.  One  way  is  to 
extract  the  thermal  energy  using  a  coolant  on  the  outer  surface  and  use  that  energy  to 
drive  turbines  as  is  done  in  the  expander  cycles.  Another  way  to  cool  is  by  film  cooling 
[10].  Film  cooling  is  done  by  injecting  a  coolant  fluid  along  the  internal  walls  of  the 
nozzle  parallel  to  the  flow  direction.  The  injected  film  forms  a  jacket  over  a  certain  finite 
length  and  prevents  the  hot  gases  from  coming  in  direct  contact  with  the  nozzle  walls. 

1.2.5.  Performance  Metrics 

There  are  three  standard  metrics  of  chemical  thruster  performance  that  are  worth 
summarizing  in  preparation  for  the  next  section  which  discusses  performance  scaling. 
The  effective  exhaust  velocity  (c)  is  defined  as  the  average  equivalent  velocity  at  which 
the  propellant  is  ejected  from  the  vehicle.  The  specific  impulse  (4^)  is  defined  as  the  total 
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impulse  per  unit  weight  of  propellant.  The  eharacteristic  veloeity  (c*)  is  a  measure  of  the 
eombustion  effieieney.  Equations  for  these  quantities  are  presented  below  where  Ae  is  the 
exhaust  exit  area,  m  is  the  total  flow  rate,  pe  is  the  exhaust  pressure,  pa,  is  the  ambient 
pressure,  and pi  is  the  eombustion  ehamber  pressure. 

c  =  U^+[p^-p^)Ajm  (1-3) 


-^  =  clgo 
mgo 


(1-4) 


c=p,Ajm  (1-5) 

c  is  a  metric  that  gives  the  overall  thrust  of  the  combustor  per  unit  flow  rate,  but  it  is  still 
a  dimensional  quantity  (dimensions  of  velocity).  In  most  places,  c  is  replaced  with  Pp 
which  indicates  thrust  per  unit  weight  flow  rate  of  fuel.  High  Pp  thrusters  can  produce 
higher  thrust  over  short  time  spans  but  are  usually  less  efficient  than  low  Pp  thrusters,  c* 
characterizes  the  burning  rate  of  the  propellant  and  is  a  characteristic  of  the  combustion 
efficiency. 


1.3.  Scaling  of  Chemical  Thruster  Performance  with  Size: 

There  are  several  factors  that  distinguish  micro-scale  combustion  from  conventional- 
scale  combustion.  These  are  described  in  the  next  three  sections  along  with  the  potential 
impact(s)  on  thruster  performance. 

1 .3  ■  1  ■  What  changes  when  you  miniaturize? 

i)  Increased  Importance  of  Heat  Loss 
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The  energy  release  in  a  eombustor  is  proportional  to  the  volume  of  the  deviee 
whereas  the  heat  losses  are  proportional  to  its  surfaee  area.  As  the  seale  of  the  eombustor 
is  redueed,  the  surfaee  area  to  volume  ratio  inereases  as  the  inverse  of  the  eharaeteristie 
length  seale.  This  is  illustrated  for  eylindrieal  and  reetangular  euboid  eombustors  in 
figure  1-4. 
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Figure  1-4:  Scale  analysis  of  microscale  combustors,  a)  Cylindrical  b)  Cuboid 

Sinee  heat  losses  seale  with  surfaee  area,  it  beeomes  apparent  that  the  heat  losses  inerease 
as  the  seale  of  the  deviee  is  redueed. 

In  early  attempts  at  eharaeterizing  eombustion  in  small  deviees,  it  was  argued  that  the 
heat  losses  would  dominate  and  henee  eombustion  would  not  be  possible  below  a  certain 
characteristic  dimension  (“quenching  diameter”).  The  quenching  diameter  was  defined  to 
be  that  diameter  at  which  the  heat  losses  from  the  surfaces  exceeded  the  amount  of  heat 
released  at  the  flame  thus  preventing  the  propagation/stabilization  of  a  flame.  Critics  of 
micro-combustion  often  reference  the  quenching  diameters  reported  in  the  classic  work 
of  Lewis  and  Von  Elbe  [11]  (where  the  concept  of  quenching  diameter  was  originally 
proposed)  as  the  smallest  possible  length  scale  for  combustion.  What  they  fail  to 
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understand  is  that  the  quenehing  diameter  is  not  a  fixed  property  of  a  fiiel-air  mixture  but 
rather  is  a  property  of  a  burner  that  is  determined  by  how  easily  heat  is  lost  to  the 
environment.  Therefore,  the  quenehing  diameters  observed  by  Lewis  and  Von  Elbe  will 
be  different  from  those  observed  by  others  in  different  burners.  What  really  matters  is  the 
eombustion  mode  and  the  amount  of  heat  loss  to  the  environment, 
ii)  Flame  Stabliziation  via  Heat  Reeireulation: 


Figure  1-5:  Illustration  of  flame  holding  in  combustors  a)  Conventional  Scale,  b) 
Microscale 

Figure  1-5  shows  how  flames  are  stabilized  in  eonventional  and  mieroseale 
eombustors.  In  a  eonventional  seale  eombustor  (Fig.  1-5  a),  the  flame  is  often  stabilized 
using  a  flame  holder  whieh  is  essentially  a  protruding  strueture  in  the  flow.  Reeireulation 
zones  form  behind  it  as  shown  by  the  eireular  arrows  and  provide  a  meehanism  for 
transferring  heat  from  eombustion  produets  to  the  ineoming  reaetants  and  eausing  them  to 
ignite  and  burn.  In  mieroseale  eombustors  (shown  in  Fig.  1-5  b),  the  eombustor  wall 
assists  in  pre-heating  the  ineoming  mixture.  The  eompetition  in  mieroseale  deviees  is 
really  between  how  much  of  the  heat  that  enters  the  structure  of  the  combustor  is  lost  to 
the  environment  and  how  much  is  recirculated. 
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The  thermal  resistanee  of  most  typieal  materials  used  to  build  eombustor  deviees  is 
several  orders  of  magnitude  lower  than  the  gas,  and  the  thermal  eoupling  between  the  gas 
and  strueture  is  mueh  stronger  at  the  miero-seale.  Thus,  the  wall  is  expeeted  to  play  a 
mueh  more  important  role  in  flame  stabilization  than  in  oonventional-seale  systems. 
Figure  1-6  illustrates  the  role  of  the  wall  in  a  mieroeombustor  and  the  eompetition 
between  axial  eonduction  and  heat  loss  to  the  environment. 


Figure  1-6:  Depiction  of  the  role  of  the  wall  in  a  microcomhustor 

It  ean  be  easily  verified  (from  a  simple  energy  balance)  that  in  such  small  combustors 
the  heat  going  into  the  wall  from  the  hot  post  combustion  gases  should  equal  the  sum  of 
the  heat  recirculated  back  to  the  pre-flame  gases  and  heat  losses.  This  is  shown  in  the 
equation  below. 


Q 


gas-wall  ,2 


Qgas-wall,\  Qi 


loss 


(1-3) 


This  was  the  missing  piece  of  information  when  “quenching  diameter”  based 
arguments  were  used  to  show  flames  cannot  be  stabilized  in  microscale  devices.  The  heat 
recirculation  actually  increases  the  flame  speed  which  in  turn  increases  the  heat  release. 
However,  this  increased  flame  speed  comes  with  the  cost  of  flame  broadening  thus 
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increasing  the  length  of  the  combustors  and  hence  increasing  the  weight.  This  will  be 
discussed  later.  This  increased  rate  of  heat  release  is  of  particular  engineering  interest  if 
we  are  able  to  make  use  of  it  to  increase  a  combustor’s  power  density, 
iii)  Low  Peclet  Number: 

Past  research  [12]  and  [13]  has  concluded  that  the  flow  Peclet  {Pe)  number  would  be 
an  important  parameter  for  distinguishing  micro/mesoscale  regimes  from  the  macroscale 
regime.  The  Pe  number  is  defined  as 


Pe  =  Re  Pr 


(  Ud\ 


(1-4) 


where  Re  is  the  Reynolds  number,  Pr  is  the  Prandtl  number,  a  is  the  thermal  diffusion 
coefficient,  v  is  the  coefficient  of  viscosity,  d  is  the  characteristic  length  (channel  width 
or  diameter  of  the  tube),  and  U  is  the  average  flow  velocity  at  the  inlet. 


Figure  1-7:  Contours  of  Pelcet  number  as  a  function  of  flow  velocity  (U)  and  passage 
height  (d)  ([14],  [15],  [16]  and  [17]) 
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Physically,  Pe  is  the  ratio  of  convective  to  conductive  heat  transfer.  In  large  Pe  flows, 
convective  heat  transfer  dominates  conductive  heat  transfer  while  in  low  Pe  flows  the 
opposite  is  true.  As  can  be  seen  from  Eq.  1-4,  the  Pe  number  includes  the  effects  of  the 
flow  velocity  and  the  passage  height  simultaneously.  Another  way  to  think  of  the  Pe  is  to 
think  of  it  as  the  ratio  of  conductive  to  convective  timescales.  Low  Pe  numbers  indicate 
large  convective  timescales  and  therefore  have  either  low  flow  velocities,  small  passages, 
or  both.  Therefore,  in  micro-scale  devices  have  low  Pe  and  heat  conduction  (both  axial 
and  transverse)  plays  an  important  role  in  heat  transfer. 

Figure  1-7  shows  contours  of  Pe  as  a  function  of  the  flow  velocity  and  passage  half 
height.  The  symbols  represent  some  of  the  micro/mesoscale  devices  that  have  been  built 
[14-17].  These  will  be  discussed  later.  Based  on  the  figure,  it  seems  reasonable  to 
conclude  that  Pe  <  30  would  characterize  most  micro-combustors.  For  the  purposes  of 
our  work  we  assume  that  Pe  numbers  below  30  characterize  micro/mesoscale  devices. 

1.3.  2.  Impact  of  Reduced  Scale  on  Thruster  Performance: 

The  concepts  introduced  in  the  previous  subsections  have  shown  that  there  are 
tradeoffs  when  we  attempt  to  scale  down  in  size.  First,  the  problem  with  heat  loss  needs 
to  be  addressed  if  we  are  to  build  viable  thrusters.  Second,  flame  speed  enhancement 
owing  to  heat  recirculation  could  be  beneficial  to  increase  thrust  and  could  offset 
increased  heat  loss.  Third,  we  need  to  have  reliable  metrics  to  show  that  performance  is 
indeed  improved  as  we  scale  down.  Fourth,  we  need  to  have  a  complete  understanding  of 
how  the  tradeoff  between  heat  losses  and  heat  recirculation  occurs  and  identify  the  main 
governing  parameters  that  influence  this  tradeoff. 
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1.3.3  Key  Barriers  to  Developing  Practical  Micro-Chemical  Thrusters 

Some  of  the  prineipal  problems  preventing  the  realization  of  high  performanee 
mieroseale  thrusters  eurrently  are:  1)  Very  short  mixing  length  requirements  2) 
Ineomplete  understanding  of  combustion  process  and  burning  rate  control  in  situations 
where  the  wall  plays  a  prominent  role  3)  Thermal  issues  related  to  structure  4)  The 
performance  of  micro-scale  nozzles  and  5)  Lack  of  instrumentation  to  make 
measurements. 

1.3.4.  Problem  of  interest  in  this  thesis 

This  dissertation  will  focus  on  understanding  flame  propagation  in  small  channels. 

1.4.  Previous  work: 

Flame  stabilization  in  a  microchannel  -  contains  elements  of  three  broad  classes  of 
physical  phenomena: 

1)  Heat  transfer:  It  has  just  been  shown  that  distinguishing  aspects  of  micro-scale 
combustion  are  the  strength  of  conjugate  heat  transfer  between  gas  and  wall  (thermal 
coupling),  heat  transfer  from  one  part  of  the  structure  to  another,  and  heat  loss  to  the 
environment. 

2)  Coupled  heat  and  mass  transfer  with  chemical  reaction:  Heat  and  mass  transfer 
processes  supply  the  energy  and  reactants  necessary  to  initiate  and  sustain  chemical 
reaction. 

3)  Practical  Applications/Limitations:  Ultimately,  our  understanding  of  flame 
stabilization  at  the  small-scale  must  be  used  to  develop  and  demonstrate  practical 
miniature  chemical  thrusters. 
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Apart  from  these,  a  fourth  related  pieee  is  the  development  of  diagnostic  techniques 
suitable  for  making  measurements  in  small  volumes  of  hot,  chemically  reacting  flows. 
Therefore,  it  is  important  to  review  previous  work  in  each  of  these  areas  in  order  to 
characterize  the  current  state  of  the  art  and  identify  the  important  research  challenges. 

1.4.1.  Heat  Transfer  Modeling 

A  simpler  view  of  the  microcombustion  problem  is  the  heat  transfer  problem 
associated  with  a  combustion  wave  in  a  duct.  This  approach  is  useful  because  it  allows  us 
to  investigate  heat  transfer  effects  separately.  Earlier  heat  transfer  analyses  in  ducts  can 
be  traced  back  to  the  classical  paper  by  Graetz  [18]  that  investigated  heat  transfer  to  a  gas 
in  a  duct  with  isothermal  walls.  He  derived  an  expression  for  the  Nusselt  number  (heat 
transfer  number)  as  a  function  of  the  streamwise  coordinate.  However,  Graetz’s 
formulation  did  not  account  for  the  axial  conduction  of  heat  through  the  gas.  This  was 
probably  because  he  was  developing  solutions  for  heat  exchanger  configurations  where 
most  of  the  heat  transfer  occurs  between  the  gas  flow  and  the  wall  in  a  direction 
perpendicular  to  the  flow  direction'.  Sellars  et  al.  [19]  extended  Graetz’s  work  to  account 
for  a  variety  of  boundary  conditions  (constant  wall  temperature,  prescribed  wall 
temperature,  and  prescribed  wall  heat  flux  etc.)  in  both  rectangular  and  circular  ducts. 
However,  axial  conduction  through  the  gas  was  still  not  included.  The  focus  was  on 
obtaining  Nusselt  number  relations  for  various  duct  geometries  with  different  boundary 
conditions.  Following  this,  analyses  by  Ebadian  and  Zhang  [20],  Michelsen  and  Villadsen 
[21],  Telles  et  al.  [22],  Eahjomri  and  Oubarra  [23]  considered  the  extended  Graetzian 

*  Note  that  the  classical  theory  for  flame  propagation  (discussed  in  §1.4.2)  is  based  on  conduction  of  heat  in 
the  streamwise  direction  via  the  gas  so  this  effect  is  important  to  include  in  reacting  flows. 
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problem  with  axial  heat  eonduction  (in  the  gas)  for  various  wall  boundary  conditions.  The 
applications  were  for  heat  exchangers  at  either  the  conventional  or  micro-scales.  A  good 
summary  of  the  work  done  in  this  area  is  given  in  the  book  by  Kakac  [24]. 

The  heat  transfer  discussed  so  far  focused  on  heat  exchangers  which  influenced 
the  types  of  boundary  conditions  imposed  at  the  walls.  In  many  cases  however,  heat 
transfer  between  the  fluid  and  wall  along  with  heat  transfer  within  the  structure  of  the 
wall  itself  are  important.  Jarosinski  [25]  analyzed  heat  transfer  between  hot  combustion 
products  and  a  cold  wall  in  ‘limit  flames’  by  developing  an  analytical  model.  The  author 
defined  ‘limit  flames’  as  those  that  are  close  to  being  extinguished  via  excessive  heat  loss 
to  the  wall.  The  heat  transfer  between  the  wall  and  the  fluid  was  modeled  by  means  of  a 
convective  transport  term.  One  of  the  chief  results  was  the  development  of  a  relationship 
between  the  Nusselt  number  (Nu)  and  the  Peclet  number  for  post  combustion  heat 
exchange  between  hot  products  and  the  cold  wall  for  plug  flow  .  This  analysis  was  quasi 
1-D  as  the  heat  transfer  between  the  gas  and  the  wall  was  modeled  based  on  a  convective 
heat  transfer  coefficient. 


^  Chapter  2  provides  the  Nu  as  a  function  of  various  flow  parameters  (from  our  heat  transfer  models)  for 
parallel  plate  channels  with  a  premixed  laminar  flame  stabilized  inside  the  channel  for  both  pre  and  post¬ 
flame  regions. 
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zPe 

Figure  1-8:  Nusselt  number  at  the  outer  surface  as  a  function  of  zPe  for  pipe 
Poiseuille  flow,  h  =  2,  Pe  =  5,  k/kf  =243  (upright  triangles),  and  k/kf  =  1  (inverted 
triangles).  The  solid  and  dashed  lines  correspond,  respectively,  to  cases  when  axial 
conduction  is  accounted  for  and  cases  when  it  is  neglected  [26] 

Yin  and  Ban  [26]  developed  an  analytical  model  for  heat  exchange  in  small  ducts  and 
channels  considering  plug  and  Poiseulle  flows.  They  solved  the  conjugate  heat  transfer 
problem  by  accounting  for  heat  exchange  between  the  wall  and  the  fluid.  Two  outer  wall 
boundary  conditions,  constant  temperature  and  constant  heat  flux,  were  considered.  The 
governing  equations  were  solved  as  an  eigenfunction  expansion  in  which  the  coefficients 
were  calculated  implicitly.  The  solution  was  used  to  relate  the  outer  boundary  Nusselt 
number  to  the  Peclet  (Pe)  number  and  the  stream  wise  coordinate. 

Figure  1-8  [31]  shows  the  variation  of  the  outer  wall  Nuh  as  a  function  of  the  non- 
dimensional  flow  coordinate  (zPe).  The  figure  shows  that  Nuh  depends  strongly  on  the 
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zPe  number  for  low  zPe  for  conductive  walls  and  that  ignoring  axial  conduction  in  the 
wall  can  give  erroneous  results.  Nm  tends  to  a  constant  value  at  large  iPe  values. 

More  recently  Maranzana  et  al.  [27]  performed  extensive  analyses  of  heat  transfer  in 
micro  and  mini  heat  exchangers  with  passage  heights  of  the  order  of  10  mm  to  a  few  mm. 
They  included  heat  transport  via  the  structure  of  the  heat  exchanger  and  introduced  a  new 
non-dimensional  number  M  to  indicate  the  fraction  of  the  total  heat  transfer  that  occurred 
via  axial  conduction  through  the  walls.  This  was  done  in  order  to  emphasize  the 
importance  of  axial  heat  transfer  through  the  walls  in  small  scale  heat  exchangers. 


Wall  conductivity  (Wm  *  K'*) 

Figure  1-9:  Exchanger  efficiency  as  a  function  of  the  wall  conductivity  [27] 

Figure  1-9  shows  the  efficiency  as  a  function  of  wall  thermal  conductivity  for  1-D  (no 
axial  heat  conduction)  and  2-D  (allowing  for  axial  conduction)  models.  The  numerator  in 
the  definition  of  efficiency  (e)  is  the  change  in  temperature  of  the  cold  fluid  flow  as  it 
passes  through  the  length  of  the  heat  exchanger  and  the  denominator  is  the  starting  point 
temperature  difference  between  the  hot  and  cold  flows.  It  can  be  clearly  seen  that  1-D 
model  over  predicts  the  efficiency.  Thus,  as  we  scale  down  in  size  it  becomes 
increasingly  important  to  account  for  the  two-dimensionality  of  the  heat  transfer  process. 
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Another  important  consideration  is  that  microcombustors  need  to  be  as  small  as 
possible.  Therefore,  flow  within  them  will  probably  not  be  fully  developed  (thermally 
and  or  hydrodynamically)  and  entrance  effects  are  expected  to  be  important.  Hsu  [28] 
studied  a  complex  entrance  effect  problem  that  accounted  for  both  the  transverse  non¬ 
uniformity  in  the  velocity  field  and  axial  conduction  in  the  gas.  In  1992  Silva  [29]  et  al. 
investigated  the  simultaneous  development  of  2-D  temperature  and  velocity  profiles  in  a 
channel  formed  between  two  infinitely  long  parallel  plates.  They  developed  a 
linearization  procedure  for  the  velocity  and  that  decoupled  the  energy  equation  and 
solved  it  using  integral  transforms.  More  recent  studies  (Nield  et  al.  [12])  have  focused 
on  thermally  developing  flow  in  porous  media. 

The  common  thread  in  all  of  the  heat  transfer  studies  described  above  is  that  they  start 
with  the  classical  Graetzian  problem  and  then  adapt  themselves  for  the  particular  problem 
of  interest.  Modern  extensions  to  the  Graetz  problem  also  usually  account  for  axial 
conduction  in  gas  and  transverse  variations  in  velocity  and  temperature. 

1.4.2.  Combustion 

i)  Conventional  Scale: 

One  of  the  earliest  theories  of  flame  propagation  (1883)  was  that  of  Mallard  and  Le 
Chatelier  [30].  It  was  based  on  a  simplified  heat  transfer  model  that  accounted  for  heat 
transfer  from  the  “hot  flame”  to  the  cold  reactants  via  conduction  through  the  gas.  In 
1941,  ZeTdovich  [31]  showed  that  there  was  a  minimum  radius  for  tubes  below  which 
premixed  flames  were  unable  to  propagate  and  extinguished;  this  got  termed  as  the 
“quenching  diameter”.  Following  this,  Spalding  [32]  developed  a  non-adiabatic  theory 
for  flame  propagation  in  ducts  that  accounting  for  heat  loss  to  the  outer  environment  by 
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adding  a  simple  loss  term  to  the  formulation.  He  showed  that  two  flame  speeds  exist  for 
every  “condition”  (namely  mixture  composition  and  heat  loss  value)  and  that  the  higher 
of  the  two  values  was  the  stable  solution.  However,  his  model  only  accounted  for  heat 
loss.  No  heat  recirculation  was  present,  and  the  predicted  flame  speeds  were  always  less 
than  the  freely  propagating  flame  speed.  It  is  noteworthy  that  this  simple  model  still 
managed  to  capture  the  qualitative  behavior  of  flame  propagation  shown  by  later 
extensive  research  (to  be  described  shortly).  In  the  late  1970’s  Takeno  et  al.  ([33]  and 
[34])  studied  excess  enthalpy  flames  first  theoretically  and  then  experimentally 
(described  later).  They  demonstrated  that  heat  recirculation  through  a  combustor’s 
structure  (in  this  case  a  porous  medium  or  a  packed  bed)  was  a  primary  mechanism  for 
stabilizing  the  flame.  “Excess  enthalpy  flames”  were  also  studied  extensively  by 
Weinberg,  a  good  account  of  which  can  be  found  in  his  text  book  [35].  He  observed  that 
heat  recirculation  made  it  possible  to  get  super-adiabatic  flames  where  the  flame 
temperature  is  significantly  higher  than  the  adiabatic  flame  temperature. 

ii)  Microscale: 

Microscale  combustion  research  is  much  more  recent  compared  to  the  classical  flame 
theories  and  experiments  just  described.  Daou  and  Matalon  (2002)  [36]  studied  the 
influence  of  heat  loss  from  the  gas  to  the  wall  on  premixed  flames  in  non-adiabatic 
channels  (i.e.  channels  that  lost  heat  to  the  environment  via  the  outer  wall).  They 
developed  an  extinction  map  that  showed  that  heat  losses  to  the  walls  could  partially  or 
completely  extinguish  the  flame.  In  2003,  Ronney  [37]  developed  a  1-D  model  of  a  heat 
recirculating  burner  and  showed  that  heat  conduction  within  the  structure  is  the  main 
process  governing  flame  stabilization.  Ju  and  Choi  in  2003  [38]  analyzed  the  role  of  heat 
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recirculation  by  studying  oppositely  propagating  flames  in  channels  sharing  a  common 
wall.  This  enabled  them  to  represent  heat  transfer  through  the  wall  using  a  quasi  ID 
model.  They  showed  that  two  flames  were  possible  and  that  the  higher  speed  was  stable. 
Cui  et  al.  [39]  considered  heat  loss  and  showed  that  the  burning  rate  of  flames  stabilized 
in  meso/microscale  channels  no  longer  depends  on  the  Lewis  number. 
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Figure  1-10:  Non-adiabatic  H2-air  combustion  in  a  silicon  micro-cbannel;  d)  =  0,5,  P 
=  1  atm,  benv  =  1  W/m^  K,  Solid  lines,  contours  of  non-dimensional  power  density 
{wolworef)  ^  function  of  non-dimensional  channel  beigbt  and  length.  Dashed 
lines,  contours  of  overall  efficiency  [40], 

Work  by  Leach  et  al.  in  2004  [40]  has  also  shown  that  axial  conduction  of  heat 
through  walls  plays  a  major  role  in  determining  micro-combustor  performance.  They 
modeled  the  coupled  heat  transfer  problem  between  wall  and  gas  using  a  thermo-electric 
analogy  and  solved  the  resulting  heat  transfer  resistance  network.  Figure  1-10  is  one  of 
their  key  results,  which  shows  that  there  is  an  optimum  passage  height  and  length 
combination  for  the  combustor  that  gives  maximum  efficiency  along  with  a  high  power 
density. 
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Wall  Thermal  Conductivity  [W/m/K] 


Figure  1-11:  Critical  velocity  vs  wall  thermal  conductivity.  The  lower  curve  (■^) 
represents  stability  loss  due  to  insufficient  heat  generation.  The  upper  curve  (■B-) 
represents  blowout.  The  shaded  region  allows  stabilized  combustion.  The 
experimentally  determined  laminar  flame  speed  is  plotted  as  a  dashed  line  [41], 

Norton  and  Vlachos  in  2004  [41]  developed  a  2-D  numerieal  simulation  of  premixed 
hydrocarbon  flames  (Methane/ Air  and  Propane/Air)  in  micro-channels  using  Fluent.  The 
simulation  accounted  for  conjugate  heat  transfer  via  the  structure  of  the  combustor  and 
reaction  chemistry  using  a  single  step  model.  The  heat  loss  to  the  outer  environment  was 
modeled  using  Newton’s  law  of  cooling  via  a  single  heat  transfer  coefficient  that 
accounted  for  all  the  heat  losses  together.  They  showed  that  there  exists  an  optimum 
thermal  conductivity  (close  to  that  of  common  ceramics)  that  gives  the  maximum  flame 
speed.  This  optimum  occurs  due  to  the  tradeoff  between  heat  loss  and  heat  recirculation. 
Figure  1-11  shows  that  the  stability  limit  is  broadened  when  the  thermal  conductivity  of 
the  wall  material  is  increased.  Using  a  simple  1-D  numerical  model  (with  full  reaction 
chemistry)  of  a  flame  held  between  two  infinitely  long  plates,  Leach  and  Cadou  [42] 
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showed  that  axial  conduction  through  the  structure  and  distributed  heat  transfer  between 
gas  and  structure  causes  the  reaction  zone  to  broaden  and  the  burning  velocity  to 
increase.  Schoegl  and  Ellzey  [43]  used  a  semi-analytical  model  to  revisit  the  problem  of 
two  flames  propagating  in  two  channels  sharing  a  common  wall  similar  to  earlier  work  by 
Ju.  They  studied  both  co-flow  and  counterflow  configurations  where  heat  conduction 
through  the  wall  separating  the  two  channels  was  accounted  for.  They  concluded  that  the 
counter-flow  configuration  was  superior  to  the  co-flow  configuration  in  that  flammability 
limits  were  considerably  broader  and  higher  and  the  flame  speeds  were  higher. 

While  these  model  problems  are  interesting  and  illuminating,  most  practical 
configurations  feature  a  single  flame  stabilized  in  a  channel  that  loses  heat  to  the 
environment.  These  are  more  challenging  to  model  because  heat  recirculation  and  heat 
losses  occur  together  and  hence  cannot  be  decoupled.  In  2007,  Barrios  et  al.  [44] 
developed  a  numerical  model  employing  two-step  chemistry  to  estimate  the  quenching 
diameter  in  small  channels.  They  concluded  that  the  heat  loss  and  not  radical  quenching 
was  the  primary  reason  for  extinction  of  flames.  More  recently  Kim  (2009)  [45]  studied 
flames  in  conductive  tubes  assuming  that  there  exists  a  critical  length  scale  upstream 
where  thermal  gradients  in  the  radial  direction  are  insignificant.  This  allowed  the 
combustor  to  be  split  into  1-D  and  2-D  regions.  It  was  shown  that  the  propagation 
velocity  of  the  flame  depends  strongly  on  parameters  such  as  the  velocity  profile  adopted 
at  the  inlet  and  the  pre-heating  of  the  inlet  stream.  Plug  and  fully  developed  flow  profiles 
(termed  in  their  work  as  “slip”  and  “no  slip”)  were  examined  along  with  three  different 
temperature  heating  values.  They  concluded  that  the  propagation  velocity  of  the  flame 
increased  as  the  inlet  temperature  was  increased  (the  equivalent  of  pre-heating  by  heat 


31 


recirculation)  and  that  the  effect  was  more  significant  for  fully  developed  flow  profiles. 
This  implies  that  the  heat  transfer  is  impacted  in  the  pre-flame  region  based  on  the  choice 
of  the  velocity  profile. 

On  the  experimental  front,  early  microcombustion  work  dates  back  to  Kotani  and 
Takeno  [46],  who  studied  microcombustion  in  a  stack  of  ceramic  tubes.  They  showed 
that  heat  recirculation  through  the  perforated  tubes  enhanced  flame  stability.  Churchill 
[47]  in  1989  did  further  experimental  studies  that  showed  that  thermally  stabilized 
combustion  was  possible  in  small  passages.  Ahn  et  al.  [15]  built  and  tested  a  swiss  roll 
combustor  and  explored  the  effect  of  using  a  platinum  catalyst  for  flame  stabilization. 
Figure  1-12  shows  the  swiss-roll  combustor  that  was  used.  They  showed  that  the  use  of 
the  catalyst  increased  the  stability  limits  and  allowed  combustion  at  very  low  Re  (-1) 
compared  to  the  situation  without  a  catalyst  (Re~40). 


Figure  1-12:  Plan  view  of  the  Swiss  roll  burner.  Platinum  catalyst  strips  were  placed 
as  shown  when  used.  Thermocouples  1-7  were  located  in  the  exhaust  ‘e’  and  inlet  ‘i’ 
turns  as  indicated  [15], 
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Ronney  et  al.  [48],  [49]  (2001-2002)  have  built  heat  reeireulating  burners  for  use  with 
thermoelectrie  generators.  The  burners  are  a  ‘Swiss  roll’  designs  and  work  by  allowing 
heat  to  be  exehanged  between  produet  and  reaetant  streams.  The  eounter-eurrent  heat 
exehanger  design  greatly  reduees  heat  loss  to  the  environment  as  most  of  the  surface  is 
exposed  to  cool  reactants  and  most  of  the  heat  goes  to  preheating  rather  than  being  lost  to 
the  environment.  Heat  recirculation  also  widens  the  flammability  limits  of  the  burner  by 
preheating  the  incoming  reactants.  While  thermoelectric  materials  are  relatively 
inefficient  in  converting  heat  to  electric  power,  the  Swiss  roll  design  has  the  advantage  of 
high  thermal  efficiency  without  the  need  for  moving  parts. 

Ju  and  Xu  [17]  studied  the  propagation  of  methane/air  and  propane/air  flames  in  a 
tube  reactor  and  measured  flame  speeds  using  a  high  speed  camera.  They  concluded  that 
both  fuels  had  multiple  extinction  limits  and  flame  regimes  and  that  the  flame  speed  had 
a  non-mono  tonic  dependency  on  equivalence  ratio.  They  also  showed  that  the  channel 
influenced  the  type  of  transition  between  slow  and  fast  flames. 
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Figure  1-13:  Quartz  heat  recirculation  reactor  (HRR),  (a)  Catalyst  and  thermal 
shields  position,  (h)  Cross-section  of  the  reactor,  (c)  Flow  direction  in  HRR  mode, 
(d)  Schematic  of  HRR  mode:  heat  flux  direction  and  thermocouple  positions,  [50] 

Scarpa  et  al.  [50]  experimentally  studied  the  effect  of  heat  recireulation  on  self 
sustaining  catalytic  combustion.  They  burned  propane/air  mixtures  in  a  quartz  reactor 
(Figure  1-13)  with  outer  tube  dimensions  of  14.6  mm  x  10.6  mm  and  inner  tube 
dimensions  of  10.6  mm  x  6.6  mm.  The  wall  thickness  of  both  tubes  was  1.3  mm.  They 
eompared  the  results  of  the  heat  recirculating  reactor  (HRR)  with  a  non  heat  recirculating 
reactor  (NRR)  and  showed  that  the  flames  experienced  enhanced  stability  in  the  HRR. 
They  also  showed  that  the  NRR  was  operating  in  the  blow  out  regime  and  was  quenehed 
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owing  to  insufficient  pre-heating.  The  eatalyst  widened  the  stability  limits  by  lowering 
the  eritieal  equivalenee. 

Taken  together,  the  survey  of  the  heat  transfer  and  eombustion  literature  (both 
theoretieal  and  experimental)  has  shown  that  2-D  heat  transfer  models  whieh  ean  aeeount 
for  eonjugate  heat  transfer  between  wall  and  gas  are  required. 

1.4.3.  Microscale  Energy  Conversion  Technologies: 

i)  Mierothrusters 

Some  of  the  previous  miero-thruster  work  found  in  the  literature  on  eoupling  of  heat 
transfer  and  eombustion  suggests  that  mieroseale  ehemieal  thrusters  would  be  strongly 
affeeted  by  heat  loss  to  the  environment.  Rudnyi  et  al  [51]  proposed  that  any  model  for 
the  operation  of  a  mierothruster  should  inelude  an  eleetro-thermal  proeess  for  ignition, 
sustained  eombustion,  membrane  rupture,  and  gas  dynamies.  They  developed  a  strategy 
for  rendering  these  sub-models  eomputationally  traetable  and  henee  praetieal  for  system 
simulation.  Table  1-5  shows  some  of  the  models  they  eonsidered  in  their  modeling 
strategy. 
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Task 

Available  models 

Modellmg  strategy' 

Electio-themial 

Ignition 

1)  Coupled  Poisson  and  heat  transfer  equations, 
Eqs  (1)  to  (4). 

2)  Lumped  resistor  and  heat  transfer.  Eqs  (3) 
and  (5). 

3)  Lumped  heat  transfer. 

The  lumped  resistor  and  3D  heat 
transfer  model  discretized  with 
the  fimte  element  method  and 
followed  by  a  model  order 
reduction. 

Ignition  and 

sustained 

combustion 

1)  Detailed  chemical  kmetics. 

2)  QSllOD. 

3)  Adiabatic  flame  temperature. 

Expenmental  igmtion 
temperature  as  input  to  the 
electro-thermal  igmtion 
problem,  expenmental  bummg 
rate  and  flame  properties  as 
mputs  to  gas  dynanucs. 

Membrane 

rupture 

1)  Molecular  and  multiscale  simulations. 

2)  Fimte  element  method. 

3)  Bulge  test  theor>'. 

Bulge  test  model  and 
expermients. 

Gas  d^manucs 

1)  Direct  Smuilation  Monte-Carlo. 

2)  Navier-Stokes  equations  for  tlie  reactive 
flow. 

3)  Local  equihbna  hypothesis  to  treat 
chenust^^^ 

4)  Fixmg  chenucal  composition  at  the  flame 
front. 

5)  Quasi- ID  approximation. 

6)  Ideal  rocket  model. 

Ideal  rocket  model  followed  by 
empincal  corrections. 

Table  1-5:  Models  for  the  microthruster  operation  [51] 

Kujawa  et  al.  [52]  modeled  a  micro-nozzle  geometry  numerically  using  FLUENT  in  2-D 
and  3-D  for  adiabatic  and  isothermal  walls.  Figure  1-14  shows  their  computational  mesh 
and  Figure  1-15  shows  path  lines  associated  with  the  fluid  flow. 


Figure  1-14:  The  three-dimensional  computational  mesh  used  by  Kujawa  et  al,  to 
simulate  a  micro-nozzle.  It  was  created  using  Fluent  Inc,’s  GAMBIT  software  [52] 
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Their  computations  showed  that  the  2-D  model  compared  well  with  quasi  1-D 
simulations.  The  3-D  simulation  results  showed  that  it  was  important  to  include  conjugate 
heat  transfer  when  modeling  these  nozzles. 


Figure  1-15:  Pathlines  for  the  3-D  model  with  isothermal  boundary  conditions  [53] 

More  recent  work  by  the  same  group  (Louisos  and  Hitt)  comparing  2-D  and  3-D  nozzle 
geometries  show  that  the  thrust  (Figure  1-16)  is  significantly  over  predicted  with  2-D 
models.  Therefore,  3-D  simulations  appear  to  be  required.  However,  their  simulations 
also  indicated  that  heat  loss  from  the  micro-nozzle  flow  increased  the  predicted  thrust  due 
to  Rayleigh  acceleration  and  the  corresponding  density  increase  as  the  flow  is  cooled. 
Heat  transfer  also  raised  the  Mach  number  in  the  vicinity  of  the  expander  walls.  This  was 
helpful  because  it  prevented  the  subsonic  boundary  layer  from  growing  large  as  was 
predicted  by  the  adiabatic  simulations. 
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Figure  1-16:  Thrust  production  in  2D  and  3D  micro-nozzles  for  varying  expander 
half-angles  from  Hitt  et  al,  [53]. 

The  2D  thrust  was  calculated  assuming  a  planar  micro-nozzle  with  depth  of  150pm.  The 
value  of  the  extrapolated  2D  thrust  is  greater  than  3D  owing  to  the  fact  that  the  2D 
simulations  do  not  account  for  the  subsonic  layers  on  the  micro-nozzle  side-walls.  For 
small  half-angles,  performance  degradation  is  due  to  the  growth  of  a  large  subsonic  layer. 
While,  at  large  expander  angles  performance  degradation  is  a  result  of  geometric  losses. 
Note  that  this  trade-off  results  in  an  optimum  expander  half-angle  of  30°  which  becomes 
more  pronounced  in  3D  simulations. 
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Figure  1-17:  Comparison  of  the  measured  data  of  nozzles  with  the  prediction  [54] 

Yang  et  al.  [54]  studied  how  the  pressure  loss  scaled  with  Re  in  micro-diffusers  and 
nozzles  by  means  of  a  simple  model  and  experiments.  Figure  1-17  shows  a  comparison 
between  their  model  and  experimentally  measured  pressure  loss  coefficients  as  a  function 
of  the  Re  number.  Their  results  indicate  that  pressure  losses  are  higher  at  lower  Re  and 
reach  an  asymptote  for  fairly  high  Re  numbers. 
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Figure  1-18:  Regeneratively  cooled  MEMS  thrust  chamber  design  to  produce  15N 
thrust  [55], 

A  bi-propellant  micro-thruster  was  developed  by  London  (Figure  1-18)  at  MIT  that 
was  based  on  the  combustion  of  O2  and  CH4.  In  spite  of  the  device’s  small  size,  viscous 
losses  did  not  prevent  the  generation  of  about  1  N  thrust  at  a  chamber  pressure  of  12 
atmospheres.  Continued  development  is  expected  to  increase  the  thrust  to  around  12-15N. 

ii)  Micro-engines  and  power  technologies: 

Next  we  examine  various  attempts  at  building  compact  power  systems  based  on 
microcombustion.  As  discussed  previously,  interest  in  microscale  heat  engines  has  been 
driven  by  the  need  for  battery  replacements  with  higher  endurance  and  miniature 
propulsion  systems  for  micro  Unmanned  Air  Vehicles  (MAVs).  More  particularly  the 


40 


focus  is  on  the  development  of  heat  engines  that  eonsume  liquid  hydroearbon  fuels  in  a 
miniature  eombustor  with  dimensions  below  1mm. 

The  seminal  program  in  this  area  that  produeed  the  first  miero-engine  was  the  miero 
gas  turbine  [14]  engine  projeet  at  MIT.  The  objective  was  to  mierofabrieate  a  gas  turbine 
engine  with  integral  eleetrie  generator  on  a  silieon  ehip.  It  was  intended  to  produee  10- 
50W  of  power. 
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Figure  1-19:  Hi  demo  engine  with  conduction-cooled  turbine  constructed  from  six 
silicon  wafers  [14], 

The  overall  size  of  the  deviee  (shown  in  Fig.  1-19)  is  about  2.1  em  x  2.1  em  x  0.38 
em.  It  is  eonstrueted  by  fusion  bonding  six  etehed  silicon  wafers  to  form  the  eompressor, 
turbine,  eombustor,  generator  and  other  flow  passages.  Individual  eomponents  of  the 
proposed  engine  have  been  demonstrated  independently  but  a  eomplete  deviee  has  not 
operated. 
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Figure  1-20:  Photograph  of  the  mini-rotary  engine  [56], 


Fernandez -Pello  and  co-workers  [56]  worked  on  a  small  rotary  engine  (Fig.  1-20) 
intended  to  produce  ~10-100  mW  of  power.  The  ‘mini’  version  was  fabricated  using 
electrodischarge  machining  (EDM),  had  a  displacement  of  0.064  mm  ,  and  a  rotor  of  size 
1mm.  The  rotary  design  has  the  advantage  of  simplicity  as  it  eliminates  the  need  for 
valves.  A  key  aspect  of  developing  any  microscale  engine  is  to  avoid  unnecessary 
moving  parts,  sealing  or  lubrication  issues  if  possible  because  higher  complexity  brings 
with  it  a  greater  risk  of  failure.  A  ‘micro’  version  of  the  rotary  engine  that  would  be 
manufactured  using  silicon  micromachining  was  also  planned. 
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Figure  1-21:  Piston-cylinder  assembly.  The  piston  is  on  the  left  and  the  end  plug  is 
on  the  right  [57]. 

Aichlmayr  and  co-workers  [57],  [58]  designed  a  homogenous  eharge  eompression 
ignition  (HCCI)  free  piston  micro-engine  (Fig.  1-21,  photograph  of  the  aetual  piston 
eylinder  assembly)  eapable  of  produeing  lOW.  HCCI  engines  differ  from  eonventional 
eompression  and  spark  ignition  engines  in  that  they  compress  premixed  fuel  and  air  until 
it  ignites  spontaneously.  This  enables  operation  at  very  low  overall  equivalenee  ratio 
whieh  can  reduce  NOx.  It  also  eliminates  the  need  (in  theory  at  least)  for  an  external 
ignition  system.  This  makes  the  design  less  eomplex.  The  free  piston  serves  as  the 
moving  element  of  a  permanent  magnet  generator 
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Figure  1-22:  Assembly  concept  for  all  non-fuel  components  of  the  20W  MICSE 
chemical-to-electrical  energy  conversion  system  indicating  major  system 
components.  System  volume  is  17,6  cc  (20x20x44  mm)  and  mass  is  54  g  [59], 

Dahm  and  coworkers  have  developed  an  oscillatory  rotating  free-piston  engine  called 
MICSE  for  Micro  Internal  Combustion  Swing  Engine.  The  engine  is  intended  to  produce 
about  20  W  with  an  overall  efficiency  of  14  %.  Butane  is  the  fuel  and  the  engine  is 
fabricated  from  steel  using  wire  EDM.  The  diameter  of  the  swing  arm  is  16  mm.  Fig. 
1.18  shows  the  integrated  assembly  drawing  of  the  MISCE  engine.  The  overall  volume  is 
expected  to  be  17.6  cc  (20x20x44  mm)  with  a  mass  of  around  54  g. 
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Figure  1-23:  Base  cavity,  swing-arm,  valve  ports,  and  resulting  four  chambers  (A-D) 
in  the  micro  internal  combustion  swing  engine  (MICSE).  The  resulting  arrangement 
produces  inherent  weight  and  volume  efficiencies  over  conventional  linear  piston 
engines,  while  avoiding  the  mechanical  complexities  of  rotary  Wankel  or  turbine 
engines  [59]. 

Fig.  1-23  shows  the  swing  arm  and  the  valve  ports.  The  meehanical  power  in  the 
form  of  the  oscillatory  motion  of  the  swing  engine  is  converted  to  electrical  power  by 
means  of  an  inductive  alternator  with  either  permanent  magnets  in  the  swing  arm  and 
inductive  coils  in  the  cavity  base,  or  by  an  external  shaft-coupled  inductive  alternator 
optimized  for  the  rotationally  oscillating  motion  produced  by  the  swing  engine.  Of  all 
the  micro-engines  in  existence,  the  swing  engine  is  by  far  the  best  developed  and  the 
closest  to  a  practical  device. 
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1.4.4.  Measurement  and  Non- intrusive  Diagnostic  Techniques 


Making  measurements  is  an  integral  part  of  eharaeterizing  performanee  of  deviees. 
Theoretieal  models  and  eomputational  results  need  validation  from  experimental  work, 
whieh  in  turn  use  various  measurement  methods  and  instruments.  The  more  intrieate  the 
engineering  problem,  the  more  ehallenging  it  beeomes  to  make  measurements. 
Mieroseale  eombustion  deviees  with  their  small  volumes,  high  temperatures  and 
ehemieally  reaeting  environments  are  extremely  ehallenging  to  make  measurements  in. 
Therefore,  there  is  a  need  to  develop  reliable  measurement  teehniques  suitable  for 
making  temperature,  eoneentration  and  veloeity  measurements  in  these  deviees.  Some  of 
the  relevant  work  in  the  literature  is  diseussed  next, 
i)  Conventional  Seale: 

Conventional  seale  measurement  teehniques  ean  be  both  intrusive  and  non-intrusive. 
As  we  will  see  in  the  next  subseetion,  mieroseale  eombustors  require  non-intrusive 
diagnostie  teehniques  owing  to  the  unique  ehallenges  of  the  environment..  One  of  the 
simplest  temperature  measuring  deviees  is  a  thermoeouple.  Shaddix  [60]  developed  a 
heat  transfer  eorreetion  model  to  improve  the  aeeuraey  of  eonventional-seale 
thermoeouple  measurements. 
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Figure  1-23:  Schematic  of  a  typical  thermocouple  arrangement  for  comhustion- 
system  measurements  [60], 

Figure  1-23  shows  the  various  heat  transfer  modes  eonsidered  by  Shaddix  [60]  from 
the  bead  of  the  thermoeouple.  His  analysis  showed  that  without  appropriate  eorreetions, 
measurement  errors  eould  be  fairly  large.  For  e.g.,  the  minimum  error  in  temperature 
measurements  due  to  ineorreet  modeling  of  eonvective  losses  eould  be  of  the  order  of  50- 
100  K.  It  will  be  signifieantly  higher  if  all  modes  of  heat  transfer  are  not  properly 
aeeounted  for.  Unfortunately,  it  is  not  possible  to  probe  miero-seale  flows  using 
eonventional  thermoeouples  as  the  diameter  of  the  bead  is  usually  -Imm  whieh  ean  be 
greater  than  the  flow  passage  height. 

In  applieations  where  it  is  not  feasible  to  insert  a  probe  beeause  the  environment  is 
too  hostile  -  i.e.  too  hot  or  oxidizing-  or  where  we  do  not  want  to  disturb  the  flow  field, 
non-intrusive  teehniques  beeome  necessary.  A  variety  of  non-intrusive  methods  are 
available  for  measuring  temperature,  composition  (species  concentration)  and  velocity 
fields  in  conventional  scale  systems.  Particle  Image  Velocimetry  (PIV)  is  a  flow 
visualization  technique  in  which  the  flow  is  seeded  with  particles  which  are  illuminated 
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by  a  laser  sheet.  Images  of  the  partieles  at  different  times  are  captured  using  a  fast  CCD 
camera  and  the  velocity  field  is  determined  by  computing  the  auto  course  [61]. 

Transmission  Fourier  transform  infrared  (FTIR)  spectroscopy  is  often  used  to  make 
non-intrusive  measurements  of  species  concentration  and  temperature.  It  involves  passing 
a  beam  of  light  through  the  sample  and  measuring  molecular  absorption  as  a  function  of 
wavenumber  using  a  detector  that  records  the  number  of  photons  exiting  the  sample 
volume  for  each  wavenumber.  In  contrast,  emission  FTIR  spectroscopy  collects  photons 
that  are  naturally  emitted  by  the  sample.  Emission  FTIR  is  often  used  in  combustion 
measurements  since  the  high  temperatures  of  the  combustion  environment  excite 
vibrational  transitions  of  most  molecules  leading  to  relatively  strong  emissions  that  are 
easy  to  collect.  Emission  experiments  [62]  have  the  advantage  that  there  is  no  beam  to 
send  through  the  sample  to  be  collected  by  a  detector.  This  makes  the  setup  much  less 
complicated,  reduces  the  number  of  optical  components  required,  and  eliminates 
problems  with  beam  steering.  Also,  this  allows  remote  sensing  measurements  of  radiative 
sources  at  a  distance.  Jaquot  et  al.  [63]  measured  temperatures  of  hot  combustion  exhaust 
gases  using  emission  spectroscopy.  They  used  the  EM2C  narrow  band  model  to  infer  gas 
temperature  from  CO2  spectra.  Measurements  of  temperatures  up  to  1100  K  have  also 
been  made.  Bailly  et  al.  [64]  used  a  high  resolution  spectrometer  (0.02  cm-1)  to  infer  gas 
temperature  from  CO2  spectra  by  only  considering  the  transitions  that  are  excited  at  high 
temperatures.  The  temperatures  were  obtained  with  precisions  on  the  order  of  ±1%  for 
this  method. 

Transmission  FTIR  is  the  most  common  form  of  FTIR  spectroscopy  because  the 
sample  does  not  have  to  be  at  the  high  temperatures  needed  for  emission  measurements. 
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In  this  technique,  a  beam  is  sent  through  the  sample  and  the  molecules  absorb  light.  One 
particularly  relevant  investigation  to  the  work  presented  here  is  that  of  Lindner  and  co¬ 
workers  [65]  who  used  FTIR  measurements  to  infer  gas  temperature,  CO2  concentration 
and  combustion  efficiency.  They  used  an  in-house  wide  band  model  with  constants 
obtained  from  the  literature  to  compute  gas  temperature  and  concentrations.  Both 
emission  and  absorbance  measurements  were  used  in  this  experiment.  Koshland  et  al  [66] 
used  in  situ  transmission  FTIR  to  detect  chlorinated  hydrocarbons  in  combustion  flows. 

While  a  variety  of  different  types  of  spectroscopies  can  be  used  to  measure 
temperature,  all  require  models  for  how  the  spectral  features  observed  depend  on 
temperature.  While  some  models  are  better  than  others,  often  sorting  out  which  is  best  is 
not  easy.  Modest  et  al.  [67]  compared  the  predictions  of  Hitran  [68],  Hi  temp  [69], 
Radcal  [70],  and  EM2C  [71]  using  an  isothermal  drop-tube  burner  and  transmission 
spectroscopy.  Their  results  indicated  that  the  wide  band  models  (Radcal  and  EM2C)  are 
better  able  to  predict  the  band  structure  at  higher  temperatures  (i.e.,  >  1300K)  with  EM2C 
being  slightly  better  than  Radcal  in  this  respect. 

ii)  Microscale: 

The  combination  of  small  volumes,  hot  oxidizing  environments,  and  the  difficulty  of 
providing  optical  access  have  prevented  the  use  of  most  conventional  combustion 
diagnostic  tools  at  the  microscale.  As  a  result,  there  are  far  fewer  examples  of 
measurements  and  diagnostic  tools  in  the  microscale.  Zhang  and  co-workers  [72] 
developed  miniature  thermocouples  that  are  micro-fabricated  into  a  combustor’s 
structure.  While  they  could  be  useful  in  low  temperature  reacting  flows  found  in  catalytic 
systems,  they  are  unable  to  survive  at  temperatures  associated  with  most  gas  phase 
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combustion  processes.  Breuer  [73]  took  advantage  of  silicon’s  transmissivity  in  the 
infrared  regime  to  measure  veloeity  using  Miero-Partiele  Image  Veloeimetry  (MPIV). 
Breuer’ s  teehnique  uses  volume  illumination  and  narrow  depth  of  field  eamera  opties  to 
eliminate  the  need  to  produee  a  sheet  of  light  with  thiekness  of  a  few  mierons. 


Fourier  Transform  Infrared  Spectrometer  (FTIR) 
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Figure  1-24:  Diagram  of  apparatus  [74] 

Heatwole  et  al.  [74]  demonstrated  the  feasibility  of  using  silieon’s  transmissivity  in 
the  IR  to  make  absorption  measurements  using  FTIR  speetroseopy.  The  measurements 
were  earried  out  in  the  silieon  walled  mieroeombustor  shown  in  Fig.  1-24.  Coneentrations 
of  CO,  CO2  and  CH4  were  determined  from  the  integrated  absorbanee  of  eaeh  moleeule’s 
rotational  band.  Gas  temperature  was  determined  from  measurements  of  rotational  line 
strengths  and  eomparison  to  predietions  of  the  HITRAN  database.  Unfortunately,  the 
uneertainties  in  the  individual  line  strengths  in  the  HITRAN  database  led  to  unaeeeptably 
large  uneertainties  in  temperature  (Fig.  1-25). 
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Figure  1-25:  Temperature  versus  downstream  distance  inferred  from  CO 
absorbance  spectra,  [74] 

These  problems  were  overcome  by  following  the  procedure  suggested  by  Modest  [67] 
which  is  to  measure  absorbance  across  the  entire  rotational  band  and  use  the  EM2C 
models  [71]  to  infer  gas  temperature  from  the  shape.  This  was  demonstrated  successfully 
in  a  microburner  [16],  [75]. 


1.5.  Objectives  and  Approach: 

The  main  objectives  of  this  dissertation  are  to  understand  the  influence  of  heat 
recirculation  through  a  combustor’s  structure  on  flame  stabilization  and  burning  rates  in 
microscale  devices,  to  identify  the  critical  combustor  design  parameters,  and  to  exploit 
this  understanding  to  develop  high  performance  microcombustors  and  thrusters.  This  will 
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be  accomplished  through  a  combined  program  of  analytical  modeling  from  hrst 

principles  and  experiments  in  a  simulated  micro-combustor 

The  approach  is  outlined  below: 

1.  Develop  2-D  analytical  models  for  heat  transfer  accounting  for  heat  exchange  with 
the  wall.  Identify  design  parameters  governing  the  heat  transfer  problem  and  explore 
the  influence  of  each  parameter  on  the  temperature  profiles  and  Nu  number  for  heat 
exchange  between  wall  and  gas. 

2.  Incorporate  a  species  transport  model  that  accounts  for  chemical  reaction  and  heat 
release  and  incorporate  it  along  with  the  heat  transfer  model  in  order  to  predict  the 
flame  speed.  Use  the  model  to  identify  the  main  governing  parameter  controlling  the 
flame  speed. 

3.  Improve  the  non-intrusive  diagnostic  technique  based  on  FT-IR  spectroscopy 
developed  in  previous  work  [16]  to  enable  measurement  of  temperature  profiles  along 
a  line  of  sight  and  extend  the  uncertainty  analysis  appropriately. 

4.  Use  the  diagnostic  technique  to  make  detailed  temperature  measurements  at  different 
flow  rates  (velocities),  equivalence  ratios,  and  for  different  fuels  (methane  and 
propane). 

5.  Use  the  temperature  measurements  to  compute  heat  recirculation  and  compare  to  the 
theory  predicted  by  models  from  approaches  1  and  2. 
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Chapter  2:  Heat  Transfer  Analysis 


In  this  chapter  we  look  at  the  heat  transfer  that  oeeurs  for  a  laminar  flame  stabilized 
within  a  eanonical  parallel  plate  reaetor.  The  focus  will  be  on  understanding  the  balance 
between  the  different  heat  transfer  meehanisms  (conveetive,  diffusive  ete)  and  also  to 
understand  the  role  of  the  wall  as  the  reaetor  is  scaled  down  in  size.  The  primary  reason 
for  beginning  this  analysis  (initially)  was  to  develop  a  model  with  which  to  obtain 
approximations  for  the  temperature  profiles  in  the  transverse  direetion  in  order  to  make 
line  of  sight  non- intrusive,  FT-IR  spectroscopy  based  temperature  measurements  by 
fitting  the  absorption  data  with  a  guess  for  the  temperature  profile  [76].  This  led  to  the 
detailed  analytical  modeling  presented  in  this  ehapter.  The  model  was  built  in  inereasing 
levels  of  eomplexity  starting  with  a  simple  problem  of  the  wall  being  at  eonstant 
temperature  and  then  advaneing  to  conjugate  heat  transfer  with  the  wall  aceounting  for 
heat  loss  to  the  environment. 

2.1.  Gas  Energy  Equation/Heat  Release  Model 

2.1.1.  Governing  Equation 

Figure  2-1  illustrates  the  basic  physical  problem.  A  flame  is  stabilized  in  the  infinitely 
ehannel  formed  between  two  parallel  plates.  Reactants  are  fed  into  the  channel  from  the 
left  side  at  exaetly  the  right  velocity  (the  flame  speed)  so  that  the  flame  position  remains 
fixed  in  the  laboratory  referenee  frame. 
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Figure  2-1:  Infinitely  thin  fiame  stabilized  in  a  channel  between  two  parallel  plates 

Hot  combustion  products  in  the  post-flame  zone  heat  the  walls  of  the  ehannel  whieh, 
in  turn,  pre-heat  the  reaetants  entering  the  reaetion  zone.  The  governing  energy  equation 
for  the  gas  (involving  gas  temperature  Tg)  in  this  problem  is  given  in  eq.  2-1  below, 
where  the  heat  release  at  the  flame  is  represented  by  a  eoneentrated  volumetrie  heat 
release  funetion  with  the  zone  of  heat  release  assumed  to  be  infinitesimally  thin. 


pU(y)Cp^  =  k^ 


V 


ax 

dx 


2 


dy^ 


\id 


(2-1) 


Where,  x  is  the  stream  wise  eoordinate  and  y  is  the  transverse  eoordinate.  r  is  the  density 
of  the  gas,  Cp  is  the  speeifie  heat  eapaeity,  kg  is  the  thermal  eonduetivity  of  the  gas,  U(y) 
is  the  veloeity  profile  and  q"g  is  the  volumetrie  heat  release  and  is  assumed  eonstant  in 
the  direetion.  Only  one  half  of  the  ehannel  is  shown  in  fig.  2-1  as  the  problem  is 
expeeted  to  be  symmetrie  about  the  centerline  in  the  absenee  of  body  forees. 

Eq.  2-1  is  re-written  in  terms  of  the  following  non-dimensional  variables: 
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(2-2) 


e  = 


T  -T 

g  CO  ^  + 


X 

~d 


;  j 


U\y^)  =  U{y^)IU 


Where 


U=\u(y^)dy^ 

0 


U;  {u(y^)  =  U  {Plug  flow)} 

=  U^Jl-  y^fl{Hagen  -  Poiseuille  flow^ 


(2-3) 


Note  that  two  velocity  profiles  across  the  channel  are  considered:  A  uniform  velocity 
distribution  (plug  flow)  and  a  parabolic  distribution  (Hagen-Poiseuille  flow).  The  non- 
dimensional  form  of  the  energy  energy  equation  becomes: 


PeUflyfl 


39  dfl  dfl 


+  ^  +  S(xfl 


dx^  dx^^  dy^^ 


(2-4) 


In  this  expression, Pe  is  the  Peclet  number  -  a  ratio  of  convective  to  conductive  heat 
transfer  as  discussed  in  Chapter  1 . 


Pe  =  RePr  = 


pUC  d 


(2-5) 


2.1.2.  Derivation  of  Jump  Conditions  across  the  Flame 

In  §  2.1.1  we  presented  the  governing  energy  equation  for  the  gas.  Now  let  us  look  a 
little  more  closely  at  the  meaning  of  the  concentrated  heat  release  function.  At  the  flame 
itself  we  require  that  the  temperatures  on  both  sides  of  the  flame  be  equal  owing  to  the 
fact  that  the  heat  release  function  is  concentrated  and  therefore  the  temperature  should  be 
continuous.  This  leads  to  the  following  equation: 

e{Q-,yfl^ey)\y^)  (2-6  a) 

Integrating  equation  2-4  over  an  infinitesimally  small  volume  around  the  flame  as  follows 
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(y  )—  =  —-T+—-T  +  ^(x  )\dx 


dx  dx  dy 


gives  the  ‘'jump  condition"  aeross  the  flame  sheet: 


(2-6  b) 


36 


--1 


(2-6  e) 


Equation  2-6  e  ean  be  understood  as  an  energy  balanee.  The  energy  released  at  the 
flame  is  dissipated  via  eonduetion  into  the  pre  and  post-flame  regions.  Together 
equations  2-6  a  and  e  are  known  as  the  “jump  eonditions”  [36,  38]  aeross  the  flame. 


2.2.  Isothermal  Wall  Boundary  Condition 

One  of  the  model  problems  found  in  the  literature  [24]  whieh  eould  be  a  reasonable 
representation  of  a  mieroeombustor  is  heat  transfer  in  a  ehannel  with  eonstant 
temperature  walls.  This  could  be  reasonable  for  a  condition  where  the  wall’s  thermal 
conductivity  will  be  several  orders  of  magnitude  larger  than  the  gas.  Hence,  the  heat 
absorbed  by  the  wall  might  be  distributed  evenly  in  a  rapid  fashion  leading  to  a  constant 
wall  temperature.  While  this  is  indeed  an  over  simplification  of  the  heat  transfer  problem, 
but  is  serves  as  a  useful  starting  point  for  the  purpose  of  comparison  with  traditional 
solutions  to  similar  heat  transfer  problems  in  channels.  With  the  assumption  of  the  walls 
being  at  a  uniform  temperature  (equal  to  the  inlet  and  exit  temperature  for  the  sake  of 
continuity  at  the  exits)  and  along  with  the  temperature  profile  in  the  gas  assumed  to  be 
symmetric  about  the  centerline  we  obtain  the  following  boundary  conditions: 

T(  +  co,y)  =  T^  (2-7  a) 

-  0  (2-8  b) 

y=0 


dy 
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(2-8  c) 


T(x,\)  =  T^=T^ 

Retaining  the  previous  non-dimensionalization  (eq.  2-2)  the  boundary  conditions 
transform  to: 


0{+co,y^)  =  0 

^  =0 

/=0 

e{x\\)  =  0 


(2-9) 

(2-10) 

(2-11) 


Figure  2-2:  Heat  transfer  problem  in  non-dimensional  variables  for  isothermal  BC, 


2.2.1.  Solution 

Figure  2-2  illustrates  the  flame  stabilization  in  an  isothermal  channel  problem.  It 
shows  the  governing  equation  and  the  wall  boundary  condition  cast  in  non-dimensional 
terms.  We  seek  a  separation  of  variables  solution  for  the  temperature  field  in  the  non- 
dimensional  streamwise  and  spanwise  coordinate  directions  andy^  respectively.  Thus, 
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0g(x\y^)  =  X{x^)Y{y^) 


(2-12) 


The  details  of  the  mathematieal  steps  involved  in  obtaining  the  solution  are  presented 
in  Appendix  A  and  ean  also  be  found  in  [77].  The  resulting  solution  is  given  below: 


exp(;L„x^)^(;L„,_y^);x^  <  0  {pre- flame) 


m=\ 


exp(-^„x^)^^(^„,_y^);x^  >  0  {post- flame) 


(2-13) 


Lm=l 


The  mathematieal  proeedure  adopted  to  obtain  An  and  Bn  in  elosed  form  (eqs.  2-14  and 
2-15)  is  presented  in  Appendiees  B  and  C.  While  the  solution  proeedure  is  unique,  it  is 
not  our  main  foeus  and  is  hence  it  is  presented  in  the  appendix. 


fi(-Pet/*(v*)  +  2A.V>* 


(2-14) 


\Wndy^ 

fi(Pet/*(v*)  +  2/?.V>* 


(2-15) 


2.2.2.  Results 

i)  Comparisson  to  CFD: 

The  predictions  of  the  analytical  model  for  Hagen-Poiseuille  flow  were  validated  via 
comparison  to  two-dimensional  CFD  simulations. 
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Figure  2-3:  Representative  computational  grid 


A  commercial  CFD  package  by  ESI  U.S.  R&D  Inc.,  CFD-ACE+  [78]  was  used  to 
solve  the  2D  steady,  incompressible  Navier-Stokes  [79]  equations  (eqs.  2-17  and  2-18) 
and  the  energy  equation  [80]  (eq.  2-16),  on  a  structured  computational  grid  based  on  the 
geometry  illustrated  in  Fig.  2-3.  The  details  of  the  grid  will  be  presented  shortly. 
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(2-16) 


(2-17) 


(2-18) 


In  Eq.  2-16,  kg^is  the  thermal  conductivity  of  the  gas,  T  is  the  absolute  temperature,  u 
and  V  are  the  x-component  and  y-component  velocities,  p  is  the  density,  P  is  the  static 


pressure,  and  t  is  the  viscous  stress,  ho  is  the  stagnation  enthalpy  which  is  defined  as: 


7  •  P 

hr,  —  l  -\ - 1 


p  2 


(2-19) 
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where  i  is  the  internal  energy.  Eqs.  2-16  through  2-19  were  solved  subjeet  to  boundary 
and  initial  eonditions  equivalent  to  those  applied  to  develop  the  analytieal  model. 

The  eomputational  grid  breaks  into  3  domains  as  illustrated  in  Fig.  2-3,  a  100mm  long 
pre-flame  domain  bounded  by  the  inlet  plane  at  the  left,  solid  walls  at  the  top  and  bottom, 
and  a  fluid  interfaee  at  the  right;  a  0.1mm  long  flame  domain  bounded  by  fluid  interfaees 
at  the  left  and  right  and  solid  walls  at  the  top  and  bottom;  and  a  100mm  long  post-flame 
region  with  a  fluid  interfaee  on  the  left,  a  fluid  outlet  on  the  right,  and  solid  walls  at  the 
top  and  bottom.  The  struetured  grid  is  5mm  x  200.1mm  and  eontains  24,745  eells.  It  is 
coarsest  at  the  inlet  and  outlet  and  most  reflned  in  the  center  (where  the  flame  is  located) 
in  order  to  ensure  that  streamwise  temperature  gradients  are  adequately  resolved. 

The  minimum  number  of  cells  required  to  provide  an  accurate  solution  was 
determined  by  performing  a  grid  convergence  test  for  a  Peclet  number  case  (Pe  =10) 
based  on  the  maximum  centerline  flow  velocity,  Umax,  and  peak  flame  zone  temperature, 
Tmax*  It  was  found  that  further  refining  a  grid  of  17,400  cells  causes  Umax  and  Tmax  to 
change  by  less  than  0.1%.  Therefore,  this  resolution  appears  to  be  adequate. 

The  boundary  conditions  at  the  inlet  are  500  K  gas  temperature  and  a  fully  developed 
parabolic  velocity  profile  {umean  =  0.02175  m/s  for  Pe=l  and  Umean  =0.2175  m/s  for 
Pe=10)  associated  with  Hagen-Poiseuille  flow.  No  slip  is  enforced  at  the  walls.  Heat 
transfer  is  allowed  between  the  walls  and  the  fluid  but  the  walls  are  assumed  to  be 
isothermal  at  500K.  The  pressure  at  the  exit  boundary  is  1  atm.  The  flame  is  represented 
as  a  constant  volumetric  heat  source  of  8.5x10  W/m  located  in  the  center  domain.  This 
value  was  identified  by  trial  and  error  as  the  one  that  produced  equal  gas  temperatures  at 
the  inlet  and  outlet  of  the  concentrated  heat  source  volume  that  represents  the  flame.  This 
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was  necessary  in  order  to  ensure  that  the  boundary  conditions  used  to  develop  the 
analytical  solution  (i.e.  that  the  pre  and  post-flame  temperatures  be  equal  at  the  flame) 
were  properly  represented.  The  fluid  density,  specific  heat,  and  thermal  conductivity  of 
the  gas  were  0.706  kg/m^,  1029  J/kgK  and  0.0395  W/mK  respectively  in  all  of  the 
simulations.  All  cells  in  the  computational  domain  were  initialized  to  the  mean  flow 
conditions  at  the  inlet.  The  Peclet  number  was  varied  from  1  to  10  by  changing  the  flow 
velocity. 

Gradients  were  computed  using  a  blended  quasi-second  order  spatial  differencing 
scheme  with  90%  upwinding  and  10%  backwards  Euler.  The  pressure,  velocity,  and 
enthalpy  were  found  using  a  conjugate-gradient-squared  (CGS)  solver  with 
preconditioning.  A  numerical  simulation  was  considered  converged  when  the  L2  norm  of 
all  residuals  was  reduced  to  at  least  1.0x10'^.  Once  a  converged  solution  was  obtained, 
the  temperature  distribution  was  exported  to  CFD-VIEW  where  the  centerline  and  local 
temperature  proflles  were  extracted.  All  simulations  were  performed  using  a  Dell  x86- 
based  PC  with  dual  3391  MHz,  Model  3  Stepping  4  Intel  processors  and  2  GB  of  RAM. 

Figure  2-4  shows  a  typical  temperature  fleld  computed  using  the  analytical  solution 
for  the  HP  flow  situation.  It  shows,  appropriately,  that  the  gas  temperature  peaks  at  the 
flame  location  (x^=0)  and  decreases  as  one  moves  upstream,  downstream,  or  to  the  walls. 
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Figure  2-4:  Non-Dimensional  2-D  temperature  (0/0f)  field  (Pe  =  1), 

Figs.  2-5  and  2-6  compare  the  analytical  and  CFD  solutions.  Fig.  2-5  shows  axial 
temperature  profiles  at  3  different  vertical  positions  in  the  channel  (y^=0,  0.4,  0.8)  for 
two  Peclet  numbers  (Pe  =  1  and  5).  These  Pe  numbers  were  chosen  because  the  flow 
velocities  and  dimensions  generally  associated  with  micro-combustors  and  burners  (our 
interest  here)  lead  to  Pe  in  the  range  of  1  to  10.  In  contrast,  Pe  associated  with 
conventional  scale  devices  are  usually  >35. 
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Figure  2-5:  Comparison  of  axial  temperature  profiles  predicted  by  the  analytical 
model  and  the  CFD  simulations  for  three  transverse  (y^  locations;  Hagen-Poiseuille 
Flow,  a)  Pe  =  1  and  b)  Pe  =  5 


Figure  2-6  shows  transverse  temperature  profiles  at  Pe  =  1  and  Pe  =  5.  The  maximum 
difference  between  the  analytical  model  and  the  numerical  solution  is  less  than  5%  over 
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the  entire  domain  for  both  Pe  numbers  indicating  that  the  analytical  solution  is  valid. 
Interestingly,  the  peak  gas  temperature  does  not  necessarily  occur  along  the  channel 
center  line  as  is  seen  in  Figs.  2-6b  and  2-7b.  The  reasons  for  this  and  the  effect  on  the 
Nusselt  number  will  be  discussed  later. 


(a) 


e/e 


(b) 

Figure  2-6:  Comparison  of  transverse  temperature  profiles  predicted  by  the 
analytical  model  and  the  CFD  simulations  for  three  axial  (x^)  locations  for  Hagen- 
Poiseuille  flow,  a)  Pe  =  1  and  b)  Pe  =  5 
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ii)  Low  Pe  number; 

Fig.  2-7  shows  the  ratio  of  the  average  non-dimensional  temperature  aeross  the 


ehannel  (in  the  transverse  or  y'^  direetion)  to  the  non-dimensional  flame  temperature  as  a 


funetion  of  axial  position  in  the  ehannel  when  the  Peelet  number  is  relatively  small. 


+ 

X 


(a) 


+ 

X 


(b) 


Figure  2-7:  Non-dimensional  temperature  profile  (0avg/0avg,F)  in  low  Fe  limit,  a)  Plug 


Flow  and  b)  Hagen-Poiseuille  Flow 
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The  mean  temperature  profile  is  symmetric  about  the  reaction  zone  and  insensitive  to 
fiow  velocity,  velocity  profile,  or  channel  height  at  low  Pe  because  the  convective 
transport  term  in  Eq.  2-4  is  small  in  comparison  to  the  conduction  terms, 
c)  Moderate  Pe  number: 

Fig.  2-8  shows  the  ratio  of  the  average  temperature  across  the  channel  (in  the 
transverse  or  y+  direction)  to  the  non-dimensional  flame  temperature  as  a  function  of 
axial  position  in  the  channel  for  moderate  values  of  the  Peclet  number.  In  both  the  plug 
and  Hagen-Poiseuille  fiow  cases,  increasing  the  Peclet  number  causes  the  axial 
temperature  distribution  to  become  more  asymmetric  about  the  flame  (heat  release)  zone. 
The  asymmetry  arises  because  the  effect  of  the  convective  term  is  different  in  the  pre  and 
post-flame  regions. 
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Figure  2-8:  Non-dimensional  temperature  profile  (0avg/0avg,F)  for  plug  and  HP  flows 
in  the  high  Pe  limit. 

To  see  this,  return  to  equation  2-4  and  eonsider  two  positions  both  a  distanee  5  from 
the  flame  but  located  in  the  pre  and  post-flame  regions  respectively.  In  the  post-flame 
region,  dO/dx^  <  0  which  means  that  the  convective  term  acts  with  (ie.  in  the  same 
direction  or  tends  to  augment)  the  axial  diffusive  term  {d^Ojdx^^  ).  This  leads  to  a  net 
increase  in  axial  thermal  transport  and  a  temperature  increase  over  the  zero  velocity  case 
at  the  arbitrary  location  x^  =  ^  in  the  post-flame.  In  the  pre-flame  region  however, 
dOjdx^  >  0  and  the  convective  term  acts  against  the  axial  diffusive  term  leading  to  a  net 
temperature  decrease  at  x^  =  -5.  Increasing  the  Peclet  number  increases  the  level  of 
distortion  by  increasing  the  difference  between  axial  diffusive  and  convective  processes 
in  the  pre  and  post-flame  regions.  This  effect  distorts  axial  temperature  profiles  in  both 
the  plug  and  Hagen-Poiseuille  flow  cases. 
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In  the  Poiseuille  flow  case,  however,  the  velocity  profile  causes  the  importance  of  the 
convective  term  relative  to  the  axial  diffusive  term  to  decrease  as  one  moves  away  from 
the  channel  center  line.  In  the  post-flame,  this  just  reduces  the  augmentation  of  the  axial 
thermal  transport  and  the  maximum  temperature  still  occurs  at  the  channel  center  and 
decreases  as  one  moves  toward  the  wall.  In  the  pre-flame,  however,  the  convective  term 
acts  against  the  axial  diffusive  term.  The  reduction  in  the  net  thermal  transport  in  the 
axial  direction  is  greatest  at  the  centerline  and  decreases  as  one  moves  toward  the  wall. 
This  means  that  the  peak  in  the  transverse  temperature  profile  does  not  occur  at  the 
channel  center  -  as  it  does  in  the  post-flame  -  but  at  some  intermediate  point  between  the 
center  and  the  wall  as  shown  in  fig.  2-7b.  Increasing  the  Peclet  number  drives  the  peak 
temperature  farther  from  the  channel  center  line  and  closer  to  the  walls.  For  Pe  >30,  the 
flame  is  essentially  a  discontinuity,  the  problem  is  mostly  one-dimensional,  and  a  two- 
dimensional  solution  is  not  required, 
iii)  Nusselt  number: 

The  Nusselt  number  is  defined  as  the  ratio  of  the  convective  heat  flux  to  the 
conductive  heat  flux.  It  is  usually  written  in  terms  of  a  heat  transfer  coefficient  h,  the 
thermal  conductivity  of  the  gas  k,  and  a  characteristic  length  scale  for  heat  transfer  which 
in  this  case  is  taken  to  be  the  hydraulic  diameter  (dc): 

Nu  = - ^ 

k  (2-20) 

The  exact  solution  for  the  temperature  profile  can  be  used  to  compute  an  overall 
Nusselt  number  for  the  problem  as  well  as  to  develop  a  correlation  for  the  variation  of  the 
Nusselt  number  with  downstream  distance.  The  latter  is  useful  in  one-dimensional 
numerical  simulations  of  reacting  flow  with  multi-step  gas  phase  chemistry  [20]  where  a 
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computationally  efficient  method  for  ineorporating  two-dimensional  heat  transfer  effeets 
is  required.  The  loeal  heat  transfer  eoeffieient  is  given  by: 


8y" 


>’"=1 


d  e. 


avg 


(2-21) 


Inserting  Eq.  2-21  into  Eq.  2-20  and  assuming  a  two-dimensional  planar  passage  {dc  = 
4  d)  gives: 


^80'^ 


Nu  = -4 


8y 


v+=l 


0. 


avg 


In  this  expression,  0avg  is  the  mean  temperature  of  the  flow: 


0avAx^)  = 


joU(y^)0(x\y^)dy^ 

\lU{y^)dy^ 


(2-22) 


(2-23) 
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Fig,  2-9:  Nu  number  vs,  axial  distance  for  different  Pe  numbers,  a)  Plug  flow,  b) 
Poiseuille  Flow 

Fig.  2-9  shows  the  variation  of  Nusselt  number  with  downstream  distance  for  a  range 
of  Pe  numbers.  In  plug  flow,  the  Nusselt  number  peaks  at  the  flame  and  rapidly  decreases 
to  the  value  7i  as  x  goes  to  ±oo.  This  is  consistent  with  other  solutions  for  heat  transfer 
between  parallel  plates  at  constant  temperature  [24].  In  Hagen-Poiseuille  flow,  however, 
the  situation  is  quite  different.  While  the  solution  also  asymptotes  to  tt  as  x  goes  to  qo,  it 
asymptotes  to  different  values  (that  depend  on  the  Peclet  number)  as  x^  goes  to  -oo.  This 
curious  behavior  in  the  pre-flame  is  another  consequence  of  the  effect  of  the  velocity 
profile  on  the  competition  between  convection  and  axial  conduction  in  the  pre-flame 
region;  Increasing  the  Peclet  number  moves  the  peak  in  the  transverse  temperature  profile 
closer  to  the  wall.  This  increases  the  Nusselt  number  by  increasing  the  temperature 
gradient  at  the  wall.  The  Nusselt  number  approaches  a  constant  far  upstream  because  the 
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temperature  gradient  at  the  wall  goes  to  zero  at  the  same  rate  as  the  mean  temperature 
goes  to  zero. 


Fig,  2-10:  Nu  asymptote  as  a  function  of  Pe  for  H-P  flow 

Fig.  2-10  shows  that  while  the  Nusselt  number  in  the  post-flame  is  a  eonstant  that  is 
independent  of  Peelet  number,  it  is  not  a  eonstant  in  the  pre-flame  when  the  Peelet 
number  is  of  order  1  or  greater;  the  horizontal  gray  line  shows  the  asymptotie  value  for 
Nu  from  traditional  heat  transfer  analysis  for  ehannel  flow  with  eonstant  wall 
temperature.  This  result  is  important  beeause  eurrently  aeeepted  praetiee  in  the  simulation 
of  reaeting  flows  in  miero-eombustors  is  to  ehoose  a  single  Nusselt  number  to  represent 
gas-wall  heat  transfer  over  the  entire  eomputational  domain.  Given  that  pre-heating  of  the 
mixture  by  the  strueture  in  the  pre-flame  region  is  a  eritieal  aspeet  of  miero-eombustor 
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performance  [40],  and  that  the  potential  variation  in  pre-flame  Nusselt  number  is  large 
(up  to  3  orders  of  magnitude  for  even  an  over  simplified  isothermal  wall  assumption),  it 
is  clear  that  incorporating  this  effect  will  be  important  in  future  simulations.  Table  2-1 
gives  the  asymptote  values  of  the  pre  and  post-flame  Nusselt  numbers  for  different  values 
of  Pe  for  H-P  flow.  These  values  may  be  used  in  future  simulations  that  model  constant 
wall  temperature  combustors  after  estimating  the  Pe  number  for  the  application.  The 
values  in  the  table  show  that  the  Nusselt  number  increases  rapidly  in  the  pre-flame  region 
beyond  a  Peclet  number  of  1 .  This  makes  sense  since  the  flow  becomes  essentially  one¬ 
dimensional  at  Pe>30. 


Pe 

NUpre 

NUpost 

0.0001 

8.11747 

8.11737 

0.001 

8.11793 

8.11692 

0.01 

8.12246 

8.11244 

0.1 

8.16984 

8.06952 

1 

8.90686 

7.79088 

5 

27.6355 

7.56931 

8 

66.4142 

7.55258 

10 

98.3997 

7.54842 

15 

193.776 

7.54418 

20 

308.596 

7.54267 

25 

440.381 

7.54196 

30 

587.348 

7.54158 

Table  2-1:  Nusselt  number  asymptote  values  for  Pre  and  Post-flame  regions  for 
different  Pe  numbers  for  H-P  flow. 
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2.3.  Conducting  Wall  Boundary  Condition  (Conjugate  Heat 


Transfer) 

In  §  2.2  we  assumed  the  wall  temperature  was  uniform  throughout  the  eombustor. 
While  this  simplified  the  solution,  and  could  be  realistic  in  some  situations,  generally  it  is 
not  physically  realistic.  Therefore,  let  us  revisit  the  problem  when  the  wall  temperature  is 
not  specified.  Fig.  2-11  shows  the  energy  equations  that  need  to  be  considered  both  in  the 
gas  and  the  structure.  Again,  the  figure  only  shows  one  half  of  the  channel  as  the  problem 
is  assumed  to  be  symmetric. 


-I- 

X 


Figure  2-11:  Half  Channel  from  Centerline  to  Outer  Edge  of  Structure, 

The  energy  equations  for  the  gas  and  structure  in  non-dimensional  form  are  given  by: 


PeU\y^) 


de^ 


dx 


+2 


+  ■ 


dy 


+2 


+  J(X^) 


(2-24) 


+  ■ 


dx 


+2 


dy 


+2 


=  0 


(2-25) 
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The  subscripts  g  and  5  stand  for  gas  and  structure  respectively,  where  S  is  the  Dirac 
delta  function.  This  implies  that  all  of  the  heat  release  occurs  at  the  infinitesimally  thin 
flame  location  (x^  =  0).  This  is  commonly  known  as  the  “flame  sheet”  assumption.  The 
non-dimensionalizations  for  the  temperature,  spatial  coordinates,  and  velocity  are  the 
same  as  for  the  previous  problem  (see  eq.  2-2)  but  two  new  non-dimensional  parameters 
K,  and  T  are  added  to  account  for  the  thermal  conductivity  and  thickness  of  the  wall 
(structure). 


ic  = 


w  +  d 
d 


■U\y^)  =  U{y^)IU 


(2-26) 


Where 


1 

U  =  \u{y^)dy^  (2-27) 

0 

The  Peclet  number  is  defined  as  before  (eq.  2-5).  The  first  two  boundary  conditions  on 
the  gas  are  the  same  as  before  (eqs.  2-9  and  2-10): 


=  0 


(2-28a) 


de^ 


=  0 


(2-28b) 


The  gas  temperature  equals  the  ambient  temperature  at  the  ends  of  the  channel  (eq.  2- 
28a)  which  is  equivalent  to  assuming  that  the  ambient  temperatures  are  the  same  at  the 
channel  inlet  and  exit  and  the  temperature  field  is  symmetric  about  the  center  line  (2- 


28b). 
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Now  that  the  structure  is  conductive,  it  is  no  longer  isothermal  and  additional  boundary 
conditions  are  required.  The  structure  temperature  is  assumed  to  go  to  the  environment 
temperature  at  +/-  infinity.  Therefore: 

=  0  (2-28c) 


The  structure  also  loses  heat  to  the  environment.  This  is  represented  by  specifying  a 
Nusselt  number  for  heat  transfer  from  the  outer  wall  to  the  environment: 


de. 


dy 


Nue  O^r) 


E_^s^ 

K 


(2-28d) 


The  Nusselt  number  for  heat  loss  to  the  environment  is  given  by:, 


Nu^=^- 


(2-29) 


This  heat  loss  model  is  based  on  Newton’s  law  for  cooling  [78].  Setting  the  heat 
transfer  coefficient  (/?£)  equal  to  zero  gives  the  adiabatic  outer  wall  condition.  Setting  hg 
>  0  incorporates  heat  loss  to  the  environment.  Finally,  continuity  of  temperature  and  heat 
flux  at  the  gas-structure  interface  requires: 


dy^ 


de^ 


=  K- 


y*=l 


dy^ 


y*=l 


(2-30a) 

(2-30b) 


Here  /ris  the  structure  to  gas  thermal  conductivity  ratio.  As  before,  we  require  that  the 
temperature  on  both  sides  of  the  flame  match  but  the  gradients  do  not  match  because  of 
the  heat  production  at  the  flame.  Integrating  over  an  infinitesimally  small  volume  around 
the  flame  gives  the  'yump  conditions'": 
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0AO-,y^)  =  0siO\y^) 


(2-3  la) 


0+ 

=  -l  (2-3  lb) 

Temperatures  in  the  strueture  must  also  mateh  at  the  flame  loeation  butintegrating  over  a 
similar  volume  gives  no  diseontinuity  in  the  strueture  temperature  beeause  there  is  no 
heat  produetion  in  the  strueture. 

0AO-,y^)  =  0siO\y^)  (2-32a) 

^  =0  (2-32b) 

dx^  0- 


2.3. 1.  Solution 

The  solution  to  the  governing  equations  is  obtained  by  using  the  same  separation  of 
variables  proeedure  used  in  the  first  problem  and  deseribed  in  Appendix  A  exeept  that 
there  is  now  are  now  two  solutions:  one  for  the  gas  and  one  for  the  strueture.  As  before, 
we  ehoose  (j)  as  the  pre-flame  solution  and  y/  as  the  post-flame  solution.  The 

eorresponding  eigenvalues  in  the  two  regions  are  given  by  X  and  /3  respeetively. 

00 

exp(2X)4(2„,/);x^  <0 

^,(x^/)=  (2-33) 


^,(x^/)=  t  (2-34) 

^n  =  \ 

Where, 
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(2-37) 


2.3 .2.  Results 

In  §2.2,  the2D  solution  for  a  flame  stabilized  in  a  eonstant  wall  temperature  ehannel 
showed  that  the  Peelet  number  was  the  only  parameter  influeneing  the  temperature 
profiles.  In  the  eonjugate  heat  transfer  problem  eonsidered  here,  four  parameters  now 
govern  the  temperature  distribution  beeause  the  wall  eonduets  heat  axially  and 
transversely.  These  are:  the  Peelet  number  (Pe),  the  ratio  of  thermal  eonduetivities 
between  wall  and  gas  (k),  the  geometrie  ratio  of  eross  seetional  areas  for  eonduetion  (t), 
and  the  heat  loss  eoeffieient  to  the  outer  environment  (Nue).  The  influenee  of  eaeh 
parameter  is  explored  below. 

i)  Effeet  of  Pe 

The  Pe  number  as  diseussed  previously  is  the  ratio  of  eonveetive  to  eonduetive  heat 
transfer.  Large  Pe  numbers  are  assoeiated  with  large  passage  heights,  large  flow 
veloeities  or  a  suitable  eombination  of  the  two.  Peelet  numbers  less  than  30  are 
representative  of  the  miero  or  meso-seale  flow  regimes  usually  assoeiated  with 
miniaturized  power  applieations. 
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Figure  2-12a:  Effect  of  Changing  Pe  on  the  Non-Dimensional  Temperature  for: 

(Plug  Flow,  K  =  10,  X  =  2,  Nue  =  10) 

(i)  Pe  =  0,1,  (ii)  Pe  =  1,  (iii)  Pe  =  5,  (iv)  Pe  =  10 

Fig.  2- 12a  shows  non-dimensional  temperature  distributions  in  the  gas  and  structure 
associated  with  four  values  of  the  Peclet  number  lying  in  the  sub  meso-scale  regime.  The 
velocity  profile  is  uniform  (plug  flow)  and  the  values  of  the  other  governing  parameters 
are  held  constant. 
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Figure  2-12b:  Effect  of  Changing  Pe  on  the  Non-Dimensional  Temperature  for: 

(Hagen-Poiseuille  Flow,  k  =  10,  x  =  2,  Nue  =  10) 

(i)  Pe  =  0,1,  (ii)  Pe  =  1,  (iii)  Pe  =  5,  (iv)  Pe  =  10 

Fig.  2- 12b  shows  the  non-dimensional  temperature  distributions  for  plane  Poiseuille 
flow.  Figures  12-12a  and  12-12b  show  a  significant  range  of  behaviors.  At  low  Pe,  the 
temperature  fields  are  spread  out  and  the  structure  temperature  is  relatively  uniform.  As 
the  Peclet  number  is  increased,  convective  heat  transfer  becomes  more  important  and 
structural  heat  conduction  is  not  able  to  smooth  out  the  temperature  distribution  in  the 
structure  as  effectively.  At  the  same  time,  the  temperature  field  in  the  channel  becomes 
more  ID  (ie.  less  variation  in  the  y^  direction).  Therefore,  the  model  shows  that  it  is 
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important  to  consider  2D  effects  when  designing  miero  and  mesoscale  eombustors  that 
operate  at  small  Pe. 

The  temperature  eontours  plotted  in  the  X-Y  plane  below  the  figure  illustrate  how 
heat  is  transported  through  the  mieroeombustor.  Sinee  heat  flux  is  the  negative  of  the 
produet  of  the  temperature  gradient  and  thermal  eonductivity,  heat  flows  in  a  direetion 
perpendicular  to  the  eontour  lines.  The  faet  that  the  axial  components  of  the  pre-flame 
temperature  gradients  in  both  the  gas  and  strueture  inerease  with  inereasing  Pe  suggest 
that  more  heat  reeireulation  is  required  to  stabilize  the  flame  as  the  flow  veloeity  is 
inereased.  The  chief  difference  between  the  plug  and  Poiseuille  flow  cases  is  the  location 
of  the  peak  temperature.  For  plug  flow,  the  peak  oecurs  along  the  eenterline  while  for 
Poiseuille  flow,  the  peak  oecurs  between  the  ehannel  eenter  and  the  wall  interfaee.  The 
consequenees  of  this  will  be  diseussed  in  more  detail  later. 

ii)  Effeet  of  k 

Another  important  parameter  governing  the  development  of  the  gas  and  strueture 
temperature  distributions  is  the  ratio  of  the  thermal  eonduetivity  of  the  strueture  to  that  of 
the  gas  (/f).  This  parameter  determines  the  preferred  heat  eonduetion  path  as  intuition 
tells  us  that  if  the  strueture  is  a  few  orders  of  magnitude  more  eonduetive,  then  the  heat 
transfer  path  of  least  resistanee  eould  be  through  the  strueture.  Figs.  2- 13a  and  b  show 
that  at  eonstant  Pe,  inereasing  /rdeereases  temperature  gradients  in  the  strueture.  In 
praetice,  it  is  hard  to  find  materials  whose  thermal  conductivities  are  much  less  than  the 
gas  (air  +  fuel)  although  /f~  1  is  possible  if  the  struetural  material  is  a  solgel. 
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Figure  2-13a:  Effect  of  Changing  k  on  the  Non-Dimensional  Temperature  for: 

(Plug  Flow,  Nue  =  10,  X  =  2,  Pe  =  5), 

(i)  K  =  0,1,  (ii)  K  =  1,  (iii)  k  =  10,  (iv)  k  =  100 


Even  when  k  is  only  moderately  large  (/t=10),  Fig.  2-13  shows  that  temperature 
gradients  in  the  strueture  are  redueed  signifieantly.  In  the  non-adiabatie  case,  the  structure 
temperature  is  also  reduced  as  the  heat  entering  the  structure  from  the  hot  gases  near  the 
flame  gets  conducted  more  efficiently  to  the  incoming  reactants  and  the  environment. 
Again,  the  chief  difference  between  Figs.  2-13  a  and  b  is  the  occurrence  of  the  peak 
temperature  at  a  different  location  from  the  centerline  in  the  case  of  of  HP 
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flow. 


Figure  2-13b:  Effect  of  Changing  k  on  the  Non-Dimensional  Temperature  for: 
(Hagen-Poisuille  Flow,  Nue  =  10,  t  =  2,  Pe  =  5), 

(i)  K  =  0.1,  (ii)  K  =  1,  (iii)  k  =  10,  (iv)  k  =  100 


iii)  Effect  of  t 

Figs.  2- 14a  and  b  show  that  while  increasing  the  wall  thickness  also  enables  greater 
heat  conduction  through  the  wall  and  hence  enhanced  thermal  interaction  between  gas 
and  structure,  changing  the  wall  thickness  is  not  equivalent  to  changing  the  thermal 
conductivity  of  the  wall. 
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(iii)  (iv) 

Figure  2-14a:  Effect  of  Changing  x  on  the  Non-Dimensional  Temperature  for: 


(Plug  Flow,  Nue  =  10,  k  =  10,  Pe  =  5) 

(i)  X  =  1.5,  (ii)  T  =  2,  (iii)  T  =  3,  (iv)  t  =  11 

This  is  because  of  the  two-dimensional  nature  of  the  problem:  Thicker  walls  reduce  the 
temperature  gradient  between  the  inner  and  outer  walls  with  the  result  that  less  heat  gets 
lost  to  the  environment  and  more  heat  gets  transported  axially  and  recirculated. 
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Figure  2-14b:  Effect  of  Changing  x  on  the  Non-Dimensional  Temperature  for: 
(Hagen-Poisuille  Flow,  Nue  =  10,  k  =  10,  Pe  =  5) 

(i)  X  =  1.5,  (ii)  T  =  2,  (ill)  T  =  3,  (iv)  x  =  11 


iv)  Effect  of  Nue 

The  final  important  parameter  governing  temperature  distributions  in  micro  and 
mesoscale  channels  is  heat  loss  to  the  environment.  Prior  work  has  shown  that  the 
balance  between  the  heat  loss  to  the  environment  and  the  heat  recirculated  upstream  via 
the  combustor’s  structure  determines  the  feasibility  of  achieving  continuous  and  stable 
chemical  reaction  [41]. 
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(iii)  (iv) 

Figure  2-15a:  Effect  of  Changing  Nue  on  Non-Dimensional  Temperature  for: 

(Plug  Flow,  X  =  2,  K  =  10,  Pe  =  5), 

(i)  Nue  =  0,1,  (11)  Nue  =  1,  (ill)  Nue  =  10,  (iv)  Nue  =  100 

Figs.  2- 15a  and  b  show  non-dimensional  temperature  distributions  in  a  micro-channel 
for  different  heat  losses  to  the  environment  described  using  the  Nusselt  number  (Nue).  As 
Nue  is  decreased  (or  heat  loss  is  reduced),  the  structure’s  temperature  increases,  the 
thermal  interaction  between  gas  and  structure  is  enhanced  in  the  pre-flame  region  and 
two-dimensional  effects  become  more  important.  While  operation  at  low  Nue  is 
advantageous  from  the  point  of  view  of  power  system  thermodynamics,  the  figure  shows 
that  low  Nusselt  number  operation  also  leads  to  steep  temperature  gradients  in  the  wall 
which  could  lead  to  thermal  degradation  of  the  structure  material. 
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Figure  2-15b:  Effect  of  Changing  Nue  on  Non-Dimensional  Temperature  for: 
(Hagen-Poiseuille  Flow,  x  =  2,  k  =  10,  Pe  =  5) 

(i)  Nue  =  0,1,  (ii)  Nue  =  1,  (iii)  Nue  =  10,  (iv)  Nue  =  100 

v)  Nu  between  structure  and  gas: 

The  solution  developed  here  can  also  be  used  to  develop  an  analytical  expression  for 
the  Nusselt  number  associated  with  heat  transfer  between  the  gas  and  the  structure  . 
Such  an  expression  is  useful  for  gaining  further  insight  into  the  heat  transfer  processes 
occurring  in  micro-combustors  as  well  as  for  developers  of  numerical  simulations  with 


^  This  is  an  important  difference  with  the  previous  quasi- ID  representation  requiring  the  prescription  of  a 
Nu  number. 


86 


multi-step  chemistry  that  need  to  save  computational  time  by  incorporating  ID  heat 
transfer  models.  The  Nusselt  number  based  on  the  hydraulic  diameter  (4  d)  is  given  by: 


Nu  +  = 

X 


dy^ 


where, 


0,{x^)  = 


iv{y')dy' 


(2-39) 


(2-40) 


As  noted  in  the  problem  formulation,  there  are  three  modes  of  heat  transport  in  this 
problem:  convection  arising  from  the  gas  flow,  conduction  through  the  gas,  and 
conduction  through  the  structure.  All  influence  the  gas  temperature  profiles  and  so  all 
effects  are  considered  simultaneously  in  the  present  work.  The  influence  of  these  various 
modes  of  heat  transfer  can  be  tracked  by  monitoring  the  Nu  between  the  gas  and  the  wall. 

Fig.  2-16  shows  that  the  axial  variation  in  Nu  is  a  strong  function  of  the  Peclet 
number.  At  low  Pe  (Pe~0.1),  convective  heat  transfer  between  the  gas  and  the  wall  and 
conduction  through  the  structure  are  much  smaller  than  axial  conduction  in  the  gas. 
Since  the  wall  is  cooler  than  the  gas  in  both  the  pre  and  post-flame  regions,  heat  transfer 
is  always  from  the  gas  to  the  wall  indicating  that  no  heat  recirculation  takes  place. 
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Figure  2-16:  Axial  variation  in  the  Nusselt  number  for  various  Pe  numbers, 
(Plug  flow,  K=  10,  T=  2,  Nue  =  10) 

At  somewhat  higher  Peclet  numbers,  a  discontinuity  appears  at  the  axial  location 
where  the  slope  of  the  temperature  gradient  at  the  wall  changes  sign.  In  the  case  of  Pe=l 
shown  in  the  figure,  this  sign  change  occurs  in  the  -4  <  <  -3  range.  Physically,  this 

means  that  the  gas  to  the  left  of  the  discontinuity  is  cooler  than  the  wall  and  hence  is 
being  pre-heated.  To  the  right  of  the  discontinuity,  the  wall  is  cooler  than  the  gas  and 
conductive  heat  transfer  through  the  gas  from  the  post-flame  to  the  pre-flame  remains  the 
dominant  mode  of  heating.  It  should  be  pointed  out  as  an  aside  that  the  discontinuity  is 
also  present  in  the  low  Pe  case  considered  earlier  but  does  not  appear  in  the  figure 
because  it  occurs  at  x"^  <  -5.  So,  some  heat  recirculation  is  still  going  on  but  it  is  very 
small. 

As  Pe  is  increased  further,  the  convective  term  in  Eq.  2-24  drives  the  discontinuity 
closer  to  the  flame.  At  Pe  >  10,  the  part  of  the  pre-flame  region  where  axial  conduction 
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through  the  gas  dominates  is  very  small  making  it  diffieult  to  distinguish  between  the 
loeation  of  the  diseontinuity  and  the  flame  loeation.  At  this  point,  effeetively  all  of  the 
gas  in  the  pre-flame  region  is  pre-heated  by  heat  reeireulation  through  the  structure. 

The  figure  also  shows  that  the  asymptotic  value  of  the  Nu  in  the  pre-flame  changes  is 
a  function  of  Pe  number.  For  comparison,  the  results  of  previous  heat  transfer  analyses 
[24]  of  forced  convective  flows  in  channels  with  two  types  of  boundary  conditions  are 
presented  below: 


Nu 


00,  Plug  Flow 


Nu 


00,  H-P  Flow 


J  9.87,  Constant  Wall  T 
[12.00,  Constant  heatflux 
{ 7.60,  Constant  Wall  T 
[8.23,  Constant  heatflux 


(2-41) 


Figure  2-17:  Comparison  of  the  axial  variation  of  the  Nusselt  number  in  plug  versus 
fully  developed  flow  for  two  different  Pe  numbers  {k=  10,  r=  2,  Nue  =  10) 
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In  plug  flow,  the  values  of  Nucc  predieted  by  the  present  model  lie  between  the  values 
for  the  eonstant  wall  temperature  and  eonstant  heat  flux  eases.  This  is  to  be  expeeted  as 
the  post-flame  region  resembles  the  model  problem  eonsidered  in  [24]  where  hot  gas 
enters  the  ehannel  and  loses  heat  to  the  walls.  In  the  pre-flame  region,  however,  the  role 
of  the  eonveetive  heat  transfer  is  different  beeause  the  axial  temperature  gradient  points 
in  the  opposite  direetion  (i.e.  has  the  opposite  sign).  This  leads  to  the  eompetition 
between  axial  and  transverse  heat  transfer  diseussed  earlier.  Fig.  2-17  eompares  the  axial 
variations  in  Nusselt  number  for  plug  and  Hagen-Poiseuille  (H-P)  flows.  Both  exhibit 
similar  behaviors  exeept  that  the  asymptote  values  for  the  H-P  flow  oases  are  lower  than 
those  of  the  oorresponding  plug  flow  oases. 


S.  No. 

Pe 

K 

X 

Nue 

Nu-oo  (Plug) 

Nuoo  (Plug) 

Nu-oo  (HP) 

NUoo  (HP) 

1 

1 

10 

2 

10 

12.09 

11.04 

3.38 

8.28 

2 

2 

10 

2 

10 

12.77 

10.80 

7.13 

7.96 

3 

5 

10 

2 

10 

14.86 

10.61 

8.71 

7.78 

4 

10 

10 

2 

10 

18.07 

10.56 

10.12 

7.74 

5 

5 

10 

2 

0.1 

13.20 

11.93 

8.57 

8.21 

6 

5 

10 

2 

1 

13.63 

11.50 

8.67 

8.06 

7 

5 

1 

2 

10 

16.24 

11.44 

7.42 

8.04 

8 

5 

0.1 

2 

10 

16.15 

11.92 

4.75 

8.21 

9 

5 

1 

2 

1 

15.74 

11.66 

8.36 

8.12 

10 

5 

10 

1.1 

10 

15.73 

10.30 

2.72 

7.70 

11 

5 

10 

3 

10 

14.04 

10.88 

8.73 

7.86 

12 

5 

10 

11 

10 

12.66 

11.62 

8.45 

8.10 

Table  2-2:  Nu  asymptotes  for  different  operating  conditions  of  the  burner 

This  is  because  the  low  speed  flow  near  the  wall  in  the  Hagen-Poiseuille  case  inhibits 
heat  transfer.  Table  2-2  gives  Nu  asymptote  values  for  some  ‘typical’  micro/mesoscale 
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burners.  The  behaviors  of  H-P  and  plug  flow  are  similar  in  the  post-flame  region  exeept 
at  low  Pe  (between  1  and  5)  where  the  H-P  asymptote  Nu  is  lower  than  that  predieted  by 
traditional  heat  transfer  analysis. For  Pe>=5  (not  shown),  the  H-P  asymptote  again  falls 
between  the  values  of  heat  transfer  analysis  for  eonstant  wall  temperature  and  eonstant 
heat  flux  boundary  eonditions.  As  diseussed  previously,  the  asymptote  value  for  the 
Nusselt  number  is  quite  different  in  the  pre-flame  region  for  both  types  of  flows.  This  ean 
be  seen  in  Table  2-2  for  Nu^. 
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Chapters:  Flame  Speed  (Eigenvalue)  Solution 


In  Chapter  2,  the  eonjugate  heat  transfer  model  aeeounting  for  heat  transfer  between 
wall  and  gas  in  2-D  for  a  pre-mixed  laminar  flame  stabilized  in  a  ehannel  was  presented 
along  with  the  analytieal  solution  for  two  veloeity  profiles  (Plug  flow  and  Poiseuille 
flow).  The  ehemieal  reaetion  term  was  modeled  as  a  eoneentrated  heat  release  funetion, 
thus  enabling  the  foeus  to  lie  on  the  pure  heat  transfer  problem.  This,  however,  over 
simplifies  the  reaetion  ehemistry  and  thus  makes  the  Pe  number  or  equivalently  the  flame 
speed  a  user  speeified  parameter.  In  real  flames,  the  ehemieal  reaetion  and  speeies 
transport  is  eonneeted  with  the  heat  transfer  and  thus  the  flame  speed  is  set  by  the 
interlinking  making  it  an  eigenvalue  of  the  problem.  Therefore,  in  this  ehapter  we 
ineorporate  the  speeies  transport  equation  with  the  previously  developed  eonjugate  heat 
transfer  model  to  obtain  the  flame  speed. 


3.1.  Governing  Equations 


Figure  3-1  is  a  sehematie  diagram  of  the  model  problem:  a  pre-mixed  flame  stabilized 
between  parallel  plates.  The  problem  is  made  non-dimensional  by  defining  the  following 
variables: 


.  -  r  „ .  r/ 

T^-T.P  '  y,.-. 

k,  w  +  d 

k  =  —;t  = - ; 


/(y)  =  k^;S,=C/;S;=^ 


+  +  y 

T  =—  y  =  — 
d  d 


L,free 


(3-1) 
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The  governing  equations  are  eonservation  of  energy  and  speeies  in  the  gas  phase: 


Pefr..  S,  f(f) 


dd„ 

+■ 


dx^  dx 


+2 


dy 


+2 


+  coS(x'^J 


(3-2) 


.  ,  dr. 

P^free  f(y  ) 


1 


dx^  Le 


dx^^  dy^ 


-coS(x^  ) 


(3-3) 


where  S  is  the  Dirae  delta  funetion  and  Le  is  the  Lewis  number  (kg/a).  The  treatment  of 
the  flame  is  the  same  as  in  Chapter  2  and  the  use  of  the  delta  funetion  implies  that  all  the 
heat  release  oecurs  at  the  infinitesimally  thin  flame  loeation  (x^  =  0).  This  is  commonly 
known  as  the  “flame  sheet”  assumption.  The  flame  is  represented  as  a  planar  structure 
with  no  curvature. 


Figure  3-1:  Schematic  illustration  of  flame  stabilized  in  a  micro-channel. 

Of  course  this  is  an  extremely  simplified  representation  of  a  two  dimensional  flame  but  it 
is  consistent  with  other  research  [38  and  43]  and,  as  we  will  see  in  the  results  section, 
this  only  affects  heat  transfer  locally  in  the  near  vicinity  of  the  flame  (from  gas  to  the 
wall).  It  does  not  affect  the  determination  of  the  eigenvalue  (flame  speed)  which  is  the 
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focus  of  this  work.  The  main  difference  between  the  model  presented  here  and  that 
presented  in  Chapter  2  is  the  inclusion  of  the  species  transport  equation  (eq.  3-3)  which 
makes  the  flame  speed  the  eigenvalue  of  the  problem  the  species  equation.  The  gas  phase 
energy  equations  are  coupled  by  the  reaction  term.  Conservation  of  energy  in  the 
structure  is  given  by: 


=  0 


The  zeroth  order  chemical  reaction  term  is  given  by: 


CO 


Ze  Pe^  ze{()p-i) 

free  ^i+y(0,-i) 


(3-4) 


(3-5) 


Pcfree  is  the  ratio  of  convective  to  conductive  heat  transfer  based  on  the  freely  propagating 
flame  speed  Si/ree  while,  Ze  and  y  are  the  Zeldovich  number  and  heat  release  parameter. 
These  are  given  by: 


Pefr,,=RePr  = 


PSlJreeCpd 


^Sijree  =  exp 


^  7^  ^ 


V 


j  _  Pg  (Pgd  P-x )  _  jPgd  P-cc) 

m2  ’  f  m 


^  ad 


^  ad 


(3-6) 


All  temperature  gradients  are  assumed  to  go  to  zero  at  +/-  infinity  leading  to  the 
following  boundary  conditions  on  the  gas  and  structure  energy  equations: 


dy^ 


=  0 


C=o 

dO, 


dy 


=  Nuj,e^{T) 


y=T 


(3-7) 


Nue  is  the  Nusselt  number  for  heat  transfer  from  the  walls  to  the  outer  environment: 
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(3-8) 


nr  hr,  d 

Nu^  =^— 


If  Nue  is  zero,  then  the  walls  are  adiabatie.  The  temperature  and  heat  flux  must  mateh  at 
the  wall-gas  interfaee  which  leads  to  the  following  additional  boundary  conditions: 


Since  the  flame  is  modeled  as  an  infinitely  thin  sheet  of  heat  release,  gas  temperatures 
on  both  sides  of  the  flame  must  match  but  the  gradients  do  not  because  of  the  heat 
production  at  the  flame.  Integrating  over  an  infinitesimally  small  volume  around  the 
flame  gives  the  ‘jump  conditions"  [38,  and  81]  across  the  flame. 


de. 


dy~ 


=  K 


de. 


/=i 


dy" 


/=i 


o,{o-y) 


de^ 

eA^\y^\  ^ 


dx 


x=0+ 


x=0 


-d) 


(3-10) 


Similarly,  we  require  that  the  temperatures  in  the  structure  match  at  the  flame  location. 
Integrating  over  a  similar  volume  gives  no  discontinuity  in  the  structure  temperature 
because  there  is  no  heat  production  in  the  structure.  The  formulation  of  the  energy 
equations  is  identical  to  that  in  Chapter  2. 


dx^ 


x=0+ 


x=0 


=  0 


(3-11) 


The  boundary  conditions  for  the  species  equations  are: 


YA-d,,yA  =  \  =  0 


dY, 


dy~ 


=  0 


r  =0,1 


(3-12) 
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These  arise  from  the  faet  that  at  the  inlet  (x^  =  -oo),  the  fuel  speeies  has  not  reaeted  or 
diffused.  At  the  flame  sheet,  we  assume  the  fuel  reaets  eompletely.  Henee,  Yf  is  zero  at 
the  flame  loeation  (x^  =  0)  and  remains  so  for  x^>  0.  At  the  eenterline  and  at  the  wall-gas 
interfaee,  the  speeies  gradients  go  to  zero  due  to  the  symmetry  and  non-permeable  wall 
boundary  eonditions.  Integrating  the  speeies  equation  over  an  infinitesimal  volume  over 
the  flame  gives  jump  eonditions  for  the  speeies  equation. 


Y^{0-,y^)  =  Y^{0\y^)  =  0 


(3-13  a) 


dx^ 


x=0 


=  -Le  CO 


(3-13  b) 


Like  the  temperature  gradient,  the  speeies  gradient  at  the  flame  is  also  set  by  the  reaetion 
rate. 


3.2.  Flame  Speed  Solution 

Equations  3-land  3-2  are  solved  using  separation  of  variables.  The  equations  are 
eoupled  through  the  temperature  and  heat  flux  matehing  boundary  eonditions.  The 
solutions  for  the  energy  equations  are  the  same  as  those  presented  in  Chapter  two. 
However,  in  order  to  obtain  the  solutions  in  the  present  ease,  we  need  to  assume  an  Sl  . 
The  temperature  fields  are  given  by: 


X  A  exp(;i„x^ )  (;i„ , / );  <  0 


n=\ 


(3-14) 


n=\ 
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exp(-/^X)</^s(A’y^y^^^  >0 


(3-15) 


where 


A.^  = 


_ _ 

£  (-  ^^/ree  K  /(3^^)  +  <^s,Jy' 


(3-16) 


= 


[i^^free  f(y^)  +  ^Pn\l,ndy^  +  ^  P ^In^y" 


(3-17) 


The  prineipal  differenee  from  the  heat  transfer  analysis  in  Chapter  2  is  that  the  flame 
speed  is  now  an  eigenvalue  of  the  problem.  A  similar  separation  of  variables  and 
applieation  of  boundary  eonditions  leads  to  the  following  expression  for  the  fuel  speeies 
distribution: 


-1  (3-18) 

0;x^  >0 


where 


C  =  — 

-1 


Leco^^i^ndy^ 

l(-Pe,„LeS-J{r)^2^,);ldy^ 


(3-19) 
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The  model  developed  here  can  accommodate  plug  flow  (U  =  constant)  or  fully 
developed  flow  profiles.  In  both  cases,  the  general  solution  procedure  is  as  follows:  1. 
Assume  a  value  for  Sl*.  2.  Compute  the  flame  temperature  Op  from  the  solution  to  the 
species  equation.  3.  Compute  the  flame  temperature  from  the  solution  to  the  energy 
equation.  4.  Adjust  Sp  and  repeat  steps  2  and  3  until  they  produce  the  same  Op. 

In  the  case  of  plug  flow,  the  solution  to  the  species  equation  reduces  from  the  series 
solution  to  a  simple  one  dimensional  form: 


1  -  exp(P£?^^^^SL  Le  x^);  <  0 

0;x^  >0 


(3-20) 


Applying  Eq.  3-20  to  the  “jump  conditions”  for  species  (Eq.  3-13  b)  and  substituting  for 
the  reaction  rate  term  gives: 


s 


* 

L 


ZePe 


free 


2 


Ze{6p-\) 


(3-21) 


Eq.  3-21  relates  ^^to  Sl  (the  flame  speed/eigenvalue). 

Moving  on  to  the  energy  equation,  the  flame  temperature  is  deflned  as  the  velocity- 
weighted  average  temperature  across  the  channel  at  the  flame  location  {x"=0): 


(3-22) 


Substituting  the  solution  to  the  energy  equation  (Eq.  3-14)  into  eq.  3-22  gives: 


uu  1  D 

I, —1  CO 


n=\ 


(3-23) 
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Note  that  cl>  is  a  known  function  of  6f  (Eq.  3-5),  Bn  are  given  by  eq.  3-17,  and  everything 
inside  the  summation  in  eq.  3-23  can  be  computed. 

In  the  case  of  Poiseuille  flow,  a  closed  form  expression  like  eq.  3-21  is  not  obtained 
but  the  approach  remains  the  same.  Applying  eq.  3- 13a  gives  an  equation  for  which 
has  Sl  embedded  in  the  terms  of  the  series  solution.  It  is  important  to  note  that  in  both 
plug  and  Poiseuille  flows,  the  Of  obtained  from  the  species  equation  has  no  dependency 
on  any  of  the  thermal  parameters  such  as  the  thermal  conductivity  ratio,  heat  loss,  or  wall 
thickness  ratio.  Rather,  it  corresponds  merely  to  the  value  required  to  achieve  this 
particular  value  of  the  flame  speed  from  the  point  of  view  of  the  chemical  reaction.  This 
fact  will  be  used  later  in  §3.3.  Finally,  the  parameters  of  the  chemical  reaction  are 
selected  to  match  experimental  conditions  in  a  Methane-Air  micro-combustion  test 
burner.  Ea  is  1.25  x  10^  J/mol-K  (Ta  =  EJRu),  Tad  is  2200  K,  the  inlet  temperature  is  300 
K,  Pcfree  is  taken  to  be  1  and  Le  is  assumed  to  be  1  for  simplicity. 

3.3.  Results 

3.3. 1.  Temperature  Profiles 

Figure  3-2  shows  a  typical  non-dimensional  temperature  field  in  the  burner.  Only  one 
half  of  the  channel  is  shown  because  of  the  symmetry  of  the  flow  problem.  Plots  such  as 
these  are  typical  of  computational  approaches  (which  can  be  performed  in  2-D)  and  it  is 
noteworthy  that  a  2-D  analytical  model  allows  us  to  plot  the  entire  temperature 
distribution  across  the  domain  while  keeping  the  computation  time  low.  This  allows  us  to 
explore  the  design  space  very  rapidly  and  to  identify  interesting  regimes. 
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X 

Figure  3-2:  2-D  Non-dimensional  temperature  plot  as  a  function  of  the  spatial 
coordinates  {k=  10,  r=  2,Nue  =  0,0375),  a)  Plug  flow;  b)  Poiseuille  flow. 

The  discontinuity  at  the  gas-structure  interface  (7=  1)  in  Figure  3-2  arises  because  of  the 

discontinuity  in  thermal  conductivity  across  the  interface.  The  temperature  gradients  in 
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the  wall  indicate  that  the  wall  participates  in  heat  transfer.  The  temperature  contours  are 
plotted  in  the  z  =  0  plane  below.  The  heat  fluxes  are  normal  to  these  constant  temperature 
contours  at  every  point.  Thus  we  can  observe  that  the  wall  has  thermal  gradients  that 
allow  both  axial  and  transverse  conduction  of  heat.  There  is  also  a  difference  in  the  shape 
of  the  temperature  distribution  between  the  Plug  and  Poiseuille  flow  cases  -  the  latter 
having  a  more  top  hat  like  shape.  This  is  the  result  of  the  competition  between  convective 
transport  and  conductive  transport  of  heat  in  the  governing  equations  as  was  discussed  in 
Chapter  2.  Figure  3-3  shows  non-dimensional  temperature  profiles  in  the  gas  and 
structure  at  various  axial  stations  in  the  pre  flame  region  for  the  same  operating 
conditions  as  in  Fig.  3-2. 


Figure  3-3:  Pre-flame  non-dimensional  temperature  distributions  in  gas  and 
structure,  a)  Plug  Flow;  b)  Poiseuille  Flow 

Far  upstream  of  the  flame,  the  gas  temperature  increases  from  the  channel  centerline  to 
the  wall  indicating  that  the  wall  is  pre-heating  the  incoming  mixture.  Near  the  flame, 
however,  the  trend  reverses:  the  peak  gas  temperature  occurs  in  the  center  of  the  channel 
because  direct  conduction  through  the  gas  is  the  most  efficient  path  for  transferring  heat 
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from  the  flame.  The  strueture  temperature  is  highest  at  the  flame  loeation  and  deereases 
in  the  deereasing  x+  direetion  beeause  of  heat  loss  to  the  environment.  The  temperature 
profiles  are  shallower  in  the  strueture  beeause  its  thermal  eonduetivity  is  mueh  greater 
than  the  gas. 


Figure  3-4:  Post-flame  non-dimensional  temperature  distributions  in  gas  and 


structure,  a)  Plug  Flow;  b)  Poiseuille  Flow 

Figure  3-4  shows  temperature  profiles  downstream  of  the  flame  for  the  same 
operating  eonditions  as  in  Figures  3-2  and  3-3.  The  peak  temperatures  remain  in  the 
eenter  of  the  ehannel  for  Plug  flow  indieating  that  heat  is  always  being  transferred  from 
the  hot,  post  eombustion  gas  to  the  strueture.  For  Poiseuille  flow  we  see  the  same  off 
eenter  temperature  peak  present  in  the  pure  heat  transfer  models  developed  in  Chapter  2. 


3.3.2.  Flame  Speed  Contours 

Figure  3-5  shows  the  non-dimensional  laminar  flame  speed  as  a  funetion  of  the 
Nusselt  number  for  heat  loss  to  the  environment  (Nue)  for  fixed  k  and  t.  The  baekward 
‘C’  shaped  eontour  shows  that  two  flame  speeds  are  possible  as  long  as  the  Nusselt 
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number  is  low  enough.  This  is  eonsistent  with  Spalding’s  [32]  eoneept  of  multiple  flame 
speeds  and  with  the  quenehing  and  blowoff  modes  observed  by  Ronney  [37]. 


Figure  3-5:  Flame  speed  as  a  function  of  heat  loss  Nusselt  number  (k  =  10,  x  =  2) 

The  figure  shows  that  for  this  partieular  eombustor  design  (k=10  and  x=2),  wall 
eonduetion  ean  inerease  the  flame  speed  by  almost  a  faetor  of  10  if  the  walls  are  well 
insulated.  Sinee  most  struetural  materials  have  thermal  eonduetivities  that  are  at  least  50 
times  larger  than  eombustion  gases,  mueh  higher  flame  speeds  appear  to  be  possible  in 
other  systems.  The  upper  and  lower  solutions  meet  at  the  extinetion  Nusselt  number 
(point  D)  whieh  represents  the  maximum  heat  loss  that  ean  be  tolerated  before  the  flame 
is  extinguished.  Only  the  upper  solution  is  stable  for  reasons  that  will  be  explained  later. 

Additional  insight  ean  be  gained  by  computing  the  heat  recirculation  through  the  gas 
and  structure  at  points  A-F  on  figure  3-5.  The  total  heat  recirculation  is  defined  as 
follows: 
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The  first  term  in  the  numerator  aeeounts  for  the  heat  reeireulation  that  oeeurs  via 


eonduetion  through  the  gas.  The  seeond  term  in  the  numerator  aeeounts  for  the  heat 
reeireulation  via  the  wall.  The  heat  reeireulation  is  made  non-dimensional  by  dividing  the 
numerator  by  the  heat  eondueted  via  the  gas  for  a  freely  propagating  flame.  Re-writing 
eq.  3-24  in  terms  of  the  non-dimensional  parameters  defined  in  eq.  3-1  gives: 
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Recirc,  structure 


(3-25) 


The  horizontal  axis  of  figure  3-6  shows  the  relative  eontributions  of  heat  reeireulation 
via  the  strueture  and  the  gas  for  eaeh  of  the  points  marked  in  fig.  3-5  while  the  vertieal 
axis  shows  the  non-dimensional  flame  speed  at  these  points.  It  shows  that  inereasing  the 
total  heat  reeireulation  inereases  the  flame  speed. 
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Figure  3-6:  Flame  speed  as  a  function  of  Heat  Recirculation  (k  =  10,  t  =  2) 

The  fraction  of  the  heat  recirculation  that  occurs  via  the  structure  also  increases  as  the 
flame  speed  increases.  This  confirms  that  heat  recirculation  through  the  structure  plays  an 
important  role  in  increasing  flame  speed.  Figure  3-7  shows  the  corresponding  heat  loss  in 
the  pre-flame  region  to  the  outer  environment.  The  maximum  heat  loss  occurs  at  the 
extinction  point  (point  D). 


105 


^Loss,  Pre-Flame 

Figure  3-7:  Flame  Speed  as  a  function  of  the  Heat  Loss  (k  =  10,  t  =  2), 

We  are  now  in  a  position  to  understand  why  only  the  upper  branch  of  the  Sl  vs  Nue 
curve  corresponds  to  stable  flames.  Consider  first  a  flame  stabilized  in  the  channel.  It 
would  be  located  somewhere  along  the  contour.  Now  let  us  perturb  the  flame  by 
increasing  the  heat  loss  (Nue)  instantaneously  by  a  small  amount.  We  would  expect  the 
increased  heat  loss  to  decrease  the  heat  recirculated  back  via  the  structure 
instantaneously.  This  should  cause  the  flame  speed  to  decrease.  If  the  flame  starts  on  the 
upper  branch,  increasing  the  heat  loss  does  indeed  drive  the  flame  speed  lower  and  the 
new  stability  point  can  lie  on  the  contour.  This  is  consistent  with  our  observation  that 
decreasing  the  heat  recirculation  decreases  the  flame  speed  and  it  will  be  possible  to 
establish  a  new  equilibrium.  If  the  flame  starts  on  the  lower  branch  however,  increasing 
the  heat  loss  would  require  that  the  new  equilibrium  point  have  a  higher  flame  speed  (in 
order  to  lie  on  the  contour).  This  is  inconsistent  with  our  observation  that  decreasing  the 
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heat  recirculation  decreases  the  flame  speed  and  it  will  not  be  possible  to  establish  a  new 
equilibrium.  Therefore,  this  flame  is  unstable. 

Having  considered  the  flame  speed  contour  for  a  particular  combustor  in  detail,  we 
now  consider  the  effects  of  changing  the  combustor  design  by  changing  r  and  k.  Figure 
3-8  shows  the  effect  of  changing  the  thermal  conductivity  ratio  for  a  fixed  value  of  the 
wall  thickness  ratio  (r=  2). 


Figure  3-8:  Burning  velocity  as  a  function  of  heat  loss  to  the  environment  for  four 
wall  thermal  conductivity  ratios  (t  =  2,0), 

As  before,  Nu^  equal  zero  corresponds  to  an  adiabatic  wall  and  the  point  of  infinite  slope 
on  each  curve  corresponds  to  the  Nusselt  number  for  extinction.  Increasing  the  thermal 
conductivity  of  the  structure  decreases  the  Nusselt  number  associated  with  extinction  but 
increases  the  magnitude  of  the  flame  speed  enhancement  by  permitting  more  heat 
recirculation  through  the  structure.  Even  when  the  wall  is  less  conductive  than  the  gas  (a 
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very  unlikely  situation  in  praetiee),  enhaneed  flame  speeds  are  still  possible  beeause  of 
the  additional  heat  reeireulation  through  the  structure. 


Figure  3-9:  Burning  velocity  as  a  function  of  heat  loss  to  the  environment  for  three 
wall  thickness  ratios  (k=10). 

Figure  3-9  is  similar  to  fig.  3-8  except  that  the  contours  correspond  to  different  wall 
thickness  ratios.  Increasing  the  wall  thickness  with  respect  to  the  flow  passage  height 
increases  the  burning  velocity  in  both  the  high  and  low  speed  modes.  It  also  increases  the 
extinction  Nusselt  number.  This  is  because  increasing  the  wall  thickness  makes  it  more 
difficult  for  heat  to  be  lost  to  the  environment  which  enables  more  heat  to  be  transferred 
upstream  to  the  pre-flame. 

3 .3 .3.  Heat  recirculation 

The  intersecting  contours  in  figs.  3-8  and  3-9  show  that  the  same  flame  speed  can  be 
achieved  using  a  variety  of  different  burner  designs  and  suggests  that  k  and  r  are  only 
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important  in  terms  of  their  effeet  on  the  net  heat  reeireulation.  To  test  this  assertion,  all 
of  the  points  used  to  ereate  figs.  3-8  and  3-9  are  plotted  against  the  associated  non- 
dimensional  heat  recirculation  in  fig.  3-10.  The  fact  that  they  appear  to  fall  on  a  single 
curve  regardless  of  whether  the  velocity  profile  is  flat  (plug)  or  fully  developed 
(Poiseuille)  and  that  the  symbols  corresponding  to  the  intersections  in  figs.  3-8  and  3-9 
overlap  each  other  in  fig.  3-10  suggests  that  heat  recirculation  is  the  primary  parameter 
that  determines  the  flame  speed. 


Figure  3-10:  Burning  velocity  as  a  function  of  the  non-dimensional  total  heat 
recirculation. 

The  primacy  of  heat  recirculation  can  be  proven  more  rigorously  as  follows:  Begin 
by  examining  each  term  of  the  total  heat  recirculation  (eq.  3-25) 
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The  inner  integral  ean  eomputed  direetly.  The  seeond  term  in  the  integral  is  zero  from 
the  definition  of  the  Dirae  delta  funetion.  Applying  eq.  3-22  gives: 

Hi,r.=s:i(f(f)0,w=s:  0,  (3-30) 

Sinee  Sl  is  a  funetion  of  Of  (see  end  of  §3.2),  the  eonverse  is  also  true  and  therefore  we 
have  proven  that  H* Redrc  only  depends  on  Sl  .  In  the  ease  of  plug  flow,  eq.  3-21  ean  be 
used  to  express  Of  as  a  funetion  of  Sl  direetly  to  get  an  analytieal  expression  relating 
H* Redrc  and  Sl  from  eq.  3-30.  This  is  the  dashed  eurve  in  figure  3-10. 

Therefore,  the  total  heat  reeireulation  appears  to  be  the  universal  governing  parameter 
that  sets  flame  speed.  This  finding  is  quite  helpful  beeause  it  ean  be  used  to  provide 
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practical  guidance  for  initial  micro-combustor  designs.  This  illustrates  why  analytieal 
models  -  though  not  as  aeeurate  as  CFD  based  numerieal  models  -  ean  still  be  very  useful 
tools  for  identifying  and  understanding  the  essential  underlying  physieal  proeesses. 


3.4.  Practical  Applications  of  the  Model: 

The  model  developed  here  ean  be  used  in  praetieal  systems  to  identify  the  range  of 
operation  and  to  make  preliminary  estimates  of  potential  improvement  of  thrust-to-weight 
ratio  over  eonventional  sized  thrusters.  The  model  ean  also  be  used  to  develop  new  non- 
intrusive  diagnostie  teehniques.  These  applieations  are  deseribed  next. 


3.4.1.  Identifying  Design  Ranges: 


We  start  by  re-writing  Eq.  3-25  as: 
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After  performing  the  inner  integral  of  the  seeond  term  and  some  more  algebraie 
rearranging  it  ean  be  shown  that: 
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The  first  term  in  parenthesis  gives  the  total  heat  flux  entering  the  pre-flame  region 
through  both  the  gas  and  strueture.  The  seeond  term  gives  the  heat  loss  to  the  outer 
environment.  It  has  already  been  shown  that  the  left  hand  side  {H* Redrc)  depends  only  on 
Sl  .  Therefore,  it  is  possible  to  ehoose  values  for  Sl  (range  of  operation)  and  one  other 
design  parameters  and  use  eq.  3-33  to  solve  for  the  third.  This  would  give  very  quiek 
‘first  eut’  design  guidanee  for  estimating  what  materials  or  wall  thieknesses  might  be 
used  to  achieve  operation  over  a  certain  range  of  flame  speeds. 

3.4.2.  Temperature  Profile  Guesses  for  Non-intmsive  Diagnostics: 

There  might  be  diagnostic  applications  (like  the  one  presented  in  Chapter  4)  or  other 
practical  system  applications  that  might  require  one  to  know  the  functional  form  of  the 
temperature  profiles  a-priori.  The  solution  presented  for  the  conjugate  heat  transfer  paper 
involves  a  Fourier  series  sum  over  infinite  eigenvalues  and  involves  complicated 
eigenfunctions  (hypergeometric  functions  etc)  and  therefore  we  seek  to  represent  these 
profiles  by  simpler  polynomial  profiles.  Previous  work  [16],  [76]  demonstrated  the 
feasibility  of  fitting  2"*^  and  4*  order  profiles  to  the  temperature  profile  of  a  much  simpler 
heat  transfer  model  for  flames  in  a  channel. 

Og(y)  =  +b^y^  +c^y^  +  d^y  +  e^  (3-34) 

We  take  the  same  approach  to  the  more  sophisticated  conjugate  heat  transfer  flame  model 
presented  here  and  in  particular  for  the  Poiseuille  flow  velocity  profile  model  which  has 
an  off  center  peak  in  the  temperature  profile.  The  idea  is  to  see  if  these  can  be 
represented  by  a  4*  order  profile  (as  the  off  center  peak  cannot  be  captured  with  a  2"*^ 
order  profile).  Figs.  3-1 1  and  3-12  compare  the  exact  solution  to  the  4**^  order  polynomial 
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fit  in  the  pre  and  post-flame  regions  of  the  flame  {xfiame  =  0).  The  temperature  profile 
appears  to  be  well  approximated  by  the  4*  order  profiles  as  ean  be  seen  visually  in  both 
pre  and  post-flame  regions. 


a 


b 


Figure  3-11:  Comparison  of  transverse  temperature  profiles  in  the 
a)  Pre-flame  and  b)  Post-flame,  for  Poiseuille  flow  {k=  10,  r=  2,Nue  =  0,0375) 
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Table  3-1  shows  the  eomparison  between  the  model  and  the  fit  showing  that  the  R 
for  the  fit  is  either  1  or  very  elose  to  1  for  all  the  profiles  that  were  fit.  Thus,  this  appears 
to  justify  the  use  of  the  4*'^  order  fit  to  approximate  the  temperature  profile.  However,  the 
slope  at  the  wall  from  whieh  heat  fluxes  would  be  estimated  had  a  larger  deviation  in  the 
near  flame  region  (around  20%).  This  suggests  that  higher  order  profiles  might  be 
neeessary  to  improve  aeeuraey  of  heat  reeireulation  estimates.  Inereasing  the  order  of  the 
profile  poses  new  problems  as  it  inereases  the  number  of  unknowns  (i.e.,  eoeffieients  for 
the  profile).  This  might  not  be  good  for  experiments  like  the  one  deseribed  in  the  next 
ehapter  where  finding  unknown  eoeffieients  is  already  quite  ehallenging  and  relatively 
few  data  points  are  available  for  fitting."^ 
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Slope  (%) 

(4“*  order  fit) 
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1 

-1 
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1.6831 
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-0.1 

-1.41753 
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19.6299 

0.9997 

1 

-0.39055 

-0.37471 

4.0558 
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-0.06098 

-0.05962 

2.2302 

1 

Table  3-1:  Comparison  between  tbe  analytical  model  and  tbe  4*'*  order  temperature 
profile  fit  for  Poiseuille  flow  {k=  10,  r=  2,  Nue  =  0,0375) 

^  The  number  of  wavelengths  over  which  absorption  is  measured  must  equal  or  exceed  the  number  of  free 
parameters  in  the  temperature  profile  in  order  for  it  to  be  possible  to  find  a  solution.  In  the  experiment 
described  in  the  next  section,  interferences  and  the  spectral  resolution  of  the  instrument  limit  the  number  of 
wavelengths  to  8.  This  means  that  the  highest  possible  fit  order  would  be  6.  This  will  be  discussed  in 
greater  detail  in  the  next  chapter. 
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3.3.4.  Thrust  Calculation 


One  of  our  main  interests  in  developing  microcombustors  is  as  part  of  a  chemical 
thruster  for  small  satellites.  Therefore,  we  would  like  to  use  the  model  presented  above  to 
make  first  order  predictions  of  thruster  performance.  In  §3.3.3  we  have  seen  that  wall 
conduction  can  significantly  enhance  the  flame  speed  in  micro/mesoscale  combustors. 
The  enhanced  mass  flow  could  mean  that  a  microcombustor/nozzle  combination  could 
produce  more  thrust  than  a  conventional  device  that  cannot  exploit  structure  thermal 
feedback. 


Nozzle  exit 


Figure  3-12:  Schematic  illustration  of  a  combustor  and  nozzle  unit. 

As  presented  previously  in  §1.2.1  the  thrust  is  computed  using  a  simple  momentum 
balance  on  the  nozzle: 


(3-34) 


In  this  expression,  the  subscript  e  denotes  exit  conditions.  The  main  aim  here  is  to 
establish  how  micro-thrusters  perform  in  comparison  to  conventional  sized  thrusters.  In 
order  to  do  this  we  write  eq.  3-34  for  both  conventional  sized  (using  freely  propagating 
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flame  values)  and  microsized  thrusters  and  take  the  ratio  of  their  respective  thrusts.  Let  us 
consider  only  a  perfectly  expanded  nozzle  (the  pressure  term  in  eq.  3-34  drops  out)  to 
achieve  maximum  thrust. 


T,  ^  pS,dU^ 

'^hjree  P  ^LJree  e,  free  j  e,  free 

In  this  expression,  d  and  D  are  the  combustors  width  for  the  micro  and  conventional  scale 
device  respectively.  If  both  combustors  have  the  same  total  pressure  and  the  nozzles  are 
perfectly  expanded  we  would  get  the  exit  velocities  to  be  the  same,  i.e.  L4  =  Ue,  free- 
Accounting  for  the  dynamic  pressure  using  the  flame  speed  can  be  expected  to  change  the 
total  pressure  in  the  combustor  for  microscale  combustors  where  flame  speeds  can  be  10- 
20  times  as  high  as  the  freely  propagating  flame  speed.  However,  the  dynamic  pressure 
component  is  still  small  compared  to  the  static  pressure  (of  even  say  1  atm)  and  this 
effectively  means  we  can  assume  that  the  static  and  total  pressures  in  the  thrust  chamber 
are  equal  in  both  micro  and  conventional  sized  thrusters.  Applying  this  to  eq.  3-35  gives 

(3-36) 

The  ratio  d/D  can  be  calculated  if  we  assume  the  wall  thickness  to  be  d  in  both  cases  and 
do  a  stress  analysis  based  on  the  hoop  stress  for  the  conventional  sized  device.  The  hoop 
stress  for  the  combustors  wall  is 


(3-35) 


(3-37) 


Where,  (Jcr  is  the  critical  stress.  The  usual  practice  is  to  add  a  factor  of  safety  (/^)  to  the 
critical  stress  value  to  allow  for  a  safety  margin.  Using  this  we  get  d/D  as. 
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^  2 

7)  “7  f/x  (3-38) 

If  we  are  interested  in  obtaining  the  thrust  ratio  sealing  (eq.  3-36),  we  ean  use  eq.  3-38  to 
represent  the  ratio  of  widths  in  eq.  3-36.  However  we  are  interested  in  the  sealing  of  the 
thrust  to  weight  ratio  and  sinee  we  are  preserving  the  shape  of  the  thruster  shape  as  we 
seale  down,  it  ean  be  shown  that  the  weight  seales  as  the  eombustors  width.  Therefore  we 
get  the  weight  ratio  between  eonventional  seale  and  mieroseale  thruster  to  be  D/d.  Using 
this  in  eq.  3-36  gives 

{T:/Wt')=Sl  (3-39) 

Thus,  for  a  perfeetly  expanded  nozzle  the  thrust  to  weight  ratio  seales  direetly  with  the 
enhaneement  in  flame  speed.  In  fig.  3-10  it  was  seen  that  it  is  possible  to  get  non- 
dimensional  flame  speeds  of  10  -  20  for  praetieal  ehoiees  of  eombustor  material  and 
thiekness.  Thus,  equation  3-39  indieates  that  we  should  expeet  inereased  thrust  to  weight 
by  a  similar  faetor. 
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Chapter  4:  Experiment:  Theory  and  Apparatus 


4.1.  Physics  of  Absorption: 

The  non-intrusive  diagnostic  technique  developed  in  this  work  is  based  on  absorption 
of  infrared  (IR)  light  by  carbon  dioxide  (CO2)  molecules  that  are  naturally  produced  in 
hydrocarbon-air  combustion.  To  understand  how  we  can  infer  important  physical 
quantities  such  as  temperature  and  concentrations  of  species  from  an  absorption 
spectrum,  we  first  need  to  understand  the  physics  behind  absorption  of  light  by 
molecules.  A  lot  of  the  details  presented  in  this  chapter  may  be  found  in  [16]  and  [81]. 
However,  they  are  presented  here  in  a  condensed  form  for  completeness  along  with  more 
recent  enhancements  to  the  technique  that  were  made  after  the  above  references  were 
published. 

4.1.1.  Qualitative  Description  of  Absorption: 

The  total  internal  energy  in  a  molecule  is  quantized  or  stored  in  discrete  bundles  that 
are  associated  with  various  vibrational,  rotational,  electronic,  and  nuclear  quantum 
mechanical  states.  The  highest  energies  are  associated  with  nuclear  states  followed  by 
electronic,  vibrational  and  rotational  states.  Electronic  states  refer  to  the  configuration  of 
the  outer  electron  shell,  vibrational  states  refer  to  the  various  vibrations  between 
constituent  atomic  nuclei,  and  rotational  states  refer  to  rotations  that  the  nuclei  make  with 
respect  to  each-other  The  overall  quantum  mechanical  state  (and  hence  total  energy)  of  a 
particular  molecule  is  described  by  its  nuclear,  electronic,  vibrational  and  rotational 
quantum  numbers  (N,  E,  V,  and  J). 
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When  a  photon  strikes  a  molecule,  it  will  only  be  absorbed  if  the  photon’s  energy 
matches  the  energy  associated  with  an  allowed  transition  to  a  higher  nuclear,  electronic, 
rotational,  and  vibrational  state.  Figure  4- la  is  a  schematic  illustration  of  the  potential 
wells  associated  with  the  ground  and  first  excited  electronic  state  for  a  typical  diatomic 
molecule.  The  blue  horizontal  lines  in  each  well  mark  the  vibrational  energies  associated 
with  different  vibrational  quantum  numbers  and  the  vertical  arrows  illustrate  two  possible 
transitions.  Energies  associated  with  rotational  states  are  not  shown  to  make  the  figure 
clearer  and  because  they  are  usually  much  smaller  than  the  vibrational  energies.  Figure 
4-lb  shows  the  energies  and  wavelengths  associated  with  various  types  of  transitions. 
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Figure  4-1:  a)  Vibrational  and  electronic  transitions  for  a  diatomic  molecule  b) 
Electromagnetic  Spectrum  and  corresponding  transitions  [83] 

Since  molecules  only  absorb  energy  in  discrete  ‘quanta’  of  energy  that  depends  on  the 
molecular  structure,  each  molecule  has  a  unique  absorption  signature.  Therefore,  it  is 
possible  to  determine  a  sample’s  chemical  composition  by  measuring  its  absorption 


119 


spectrum.  IR  light  being  “low  energy”  (or  longer  wavelength)  light  primarily  excites 
vibrational  and  rotational  transitions  owing  to  their  low  energy  states  (see  Fig.  4-lb). 
Absorption  of  higher  energy  photons  in  the  ultraviolet  region  is  usually  associated  with 
electronic  in  addition  to  vibrational  and  rotational  transitions. 

The  strength  and  shape  of  a  molecule’s  absorption  spectrum  is  related  to  its 
concentration  and  temperature.  The  concentration  sets  the  number  of  molecules  available 
to  absorb.  The  temperature  sets  the  distribution  of  initial  energy  states  through  the 

-BJ(J+1) 

Boltzmann  distribution  {^^  =  (2J  +  \)e  ,  see  [84]  pp.  445-449  for  more  details). 

Figure  4-2  shows  that  raising  the  gas  temperature  raises  the  average  rotational  quantum 
number.  This  means  that  samples  at  higher  temperatures  will  have  more  molecules  in 
higher  energy  states. 
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Figure  4-2:  Relative  populations  of  the  rotational  levels  of  HCl(g)  at  three  different 
temperatures,  (Adapted  from  [84]) 
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Since  the  absorption  speetrum  depends  on  the  distribution  of  energy  states  in  the  gas, 
it  is  a  direet  indieator  of  temperature.  Therefore,  measuring  the  absorption  spectrum 
provides  information  on  both  temperature  and  eoneentration.  A  more  detailed  deseription 
of  the  absorption  physios  oan  be  found  in  [84], 


Figure  4-3:  Rotational  bands  in  CH4  spectrum. 

Fig.  4-3  shows  the  absorption  speetrum  of  CH4.  The  spectral  resolution  of  the  plot  is 
low  whioh  means  only  the  outer  envelope  of  the  vibrational  band  ‘head’  (which  is 
composed  of  individual  rotational  lines)  is  visible.  The  rotational  lines  oorrespond  to 
transitions  [84]  and  [85]  with  AJ  =  -1,  0  and  1  respeetively  where  J  is  the  rotational 
quantum  number.  These  are  called  respeetively  the  P,  Q,  and  R  branehes. 

In  this  work,  we  are  primarily  eoneerned  with  the  infrared  speetrum  of  the  CO2 
moleeule  which  arises  from  the  four  vibrational  modes  illustrated  in  figure  4-4.  CO2  is  a 
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linear  molecule  with  two  stretching  (symmetric  and  asymmetric)  and  two  bending 
(degenerate)  modes.  The  symmetric  stretching  mode  is  excited  by  photons  at  2350  cm'\ 


4.1.2.  Review  of  Absorption  Spectroscopy: 

The  Beer-Lambert  law  [86]  relates  the  change  in  intensity  of  a  beam  of  light  passing 
through  a  gas  column  of  length  dx  to  the  concentration  C  of  the  absorbing  gas  molecules 
in  moles/liter: 

^^  =  -k{v)Cdx  (4-1) 

I{v) 

In  this  equation,  /  is  the  initial  light  intensity;  dl  is  the  change  in  transmitted  intensity  due 
to  absorption  by  the  sample  along  the  optical  path  of  length  dx,  v  is  the  optical  frequency 
(typically  in  cm'^  for  infrared  measurements)  and  k{y)  is  an  absorption  coefficient  that 
depends  on  the  optical  frequency.  Integrating  over  the  total  length  of  the  optical  path  L 
and  converting  natural  logarithms  to  base  10  leads  to  the  definition  of  the  molecular 
absorbance  A(v): 
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logio  =  4^  )  =  )CI 

I[v) 


(4-2) 


In  this  expression,  the  absorption  eoeffieient  {k)  has  been  replaeed  with  the  molar 
absorptivity  s  (^(k)  =  2.306A:(v))  whieh  is  also  a  funetion  of  optieal  frequeney.  The  faet 
that  the  absorbanee  is  proportional  to  the  path  length  L  presents  a  ehallenge  at  the  miero- 
seale  where  path  lengths  are  redueed  to  a  few  millimeters  or  less.  This  is  offset,  however, 
by  the  faet  that  many  of  the  speeies  of  interest  in  this  work  (CO,  CO2,  H2O,  and  CH4)  are 
present  in  relatively  large  eoneentrations. 

While  Eq.  (4-2)  eould  be  used  to  infer  eoneentration  from  an  absorption  measurement 
made  at  a  single  wavelength  assoeiated  with  a  single  ro-vibrational  transition,  more 
reliable  results  are  obtained  if  measurements  of  absorption  over  many  transitions  are 
eonsidered.  As  a  result,  eoneentration  measurements  are  usually  made  by  measuring  the 

integrated  absorbanee  ,  whieh  is  simply  the  integral  of  the  absorbanee  over  the 
entire  absorption  band: 

Aband  =  I s{y)CL  dv  =  sCL  (4-3) 

hand 

Raising  the  gas  temperature  broadens  the  Boltzmann  distribution  of  energies  in  the  gas 
(ie.,  a  wider  range  of  energy  levels  are  signifieantly  populated)  but  the  area  under  the 
distribution  remains  eonstant  beeause  total  number  of  moleeules  is  fixed.  Therefore,  the 
total  area  under  the  speetrum  ( )  is  insensitive  to  temperature  beeause  it  is  a  funetion 
of  the  number  of  moleeules  present  or  the  eoneentration  but  the  shape  of  the  absorption 
band  (  A(v) )  ehanges  appreeiably  with  temperature  beeause  the  shape  of  the  Boltzmann 
distribution  ehanges  appreeiably  with  temperature.  This  work  exploits  the  faet  that  the 
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shape  of  the  absorption  band  is  a  known  function  of  temperature  to  infer  the  gas 
temperature  from  the  absorption  spectrum. 


Figure  4-5:  Sample  hot  versus  cold  CO2  spectra,  measured  at  16  cm'^  resolution 

There  are  two  approaches  to  modeling  the  temperature  dependence  of  the  absorption 
spectrum.  High  resolution  models  like  HITRAN  [68]  and  HITEMP  [69]  provide  data  for 
the  location  and  temperature  dependence  of  individual  spectral  lines  associated  with 
particular  ro-vibronic  transitions.  These  individual  lines,  when  added  together  and  viewed 
at  relatively  coarse  resolution  (Av  ~16  cm’'),  form  the  bands  as  shown  in  Fig.  4-4.  In 
contrast,  low  resolution  models  like  EM2C  [71]  describe  the  temperature  dependence  of 
the  overall  band  shape.  While  high  resolution  models  provide  much  more  spectral  detail, 
they  require  detailed  knowledge  of  the  temperature  dependence  of  individual  spectral 
lines.  Unfortunately,  reliable  data  is  not  available  for  most  lines  at  temperatures  greater 
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than  about  1200K.  Recent  work  [67]  has  shown  that  narrow  band  models  are  more 
reliable  for  inferring  temperature  from  infrared  spectra  in  combustion  flows,  and  so  that 
is  the  approach  taken  here. 

The  narrow  band  model  is  implemented  using  EM2C,  a  FORTRAN  program  with 
associated  spectral  databases  developed  by  EM2C  laboratories  to  compute  the  radiative 
intensity  and  transmissivity  of  a  non-isothermal,  non-homogeneous  gas  column 
consisting  of  mixtures  of  gases  like  H2O,  CO2,  CO  etc.  These  calculations  are 
accomplished  by  discretizing  the  optical  path  into  n  homogeneous  isothermal  elements 
using  the  Malkmus  [87]  statistical  narrow  band  model  and  the  Curtis-Godson 
approximation  at  a  spectral  resolution  of  25  cm'^  The  temperature  fitting  procedure  is 
described  in  detail  later. 

The  Malkmus  model  assumes  Lorentz  line  shapes  averaged  over  a  wide  spectral 
interval  and  that  the  individual  column  elements  are  homogenous  and  isothermal.  For  a 
given  column  length  L  (cm),  total  pressure  p  (atm),  and  mole  fraction  X  of  absorbing  gas, 
the  transmissivity  f^''  over  the  spectral  interval  Av  (in  this  case  25  cm'^)  is  given  by: 

=  exp 

In  this  expression,  k  (cm"'  atm"'),  (7(cm'^)  and  y(cm'^)  are  the  model  parameters  for 
the  particular  molecular  species  under  consideration.  The  Curtis-Godson  approximation 
involves  replacing  k  and  //Sin  eq.  4-1  with  average  values  of  these  parameters  over  the 
particular  spectral  interval  being  considered.  The  input  parameters  for  the  EM2C  model 
come  from  the  HITRAN  database  and  supplementary  data  for  k,  1/5,  and  y  obtained  by 


■2^ 
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r) 


-1 


(4-4) 
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EM2C  laboratories.  In  the  work  reported  here,  we  foeus  on  absorption  by  CO2  sinee  it  is 
naturally  produeed  in  the  flame  and  gives  a  relatively  strong  signal. 

4.1.3.  Gas  Temperature  Measurement  Principle: 

Gas  temperature  measurements  are  made  by  fltting  y4(v)for  CO2  predieted  by  the 
EM2C  eode  to  A(v)  measured  from  the  experiments.  This  is  aeeomplished  in  MATLAB 
using  a  Newton-Gauss  algorithm  and  a  least  squares  fitting  routine  Tsqeurvefif .  The 
Newton-Gauss  algorithm  ealls  the  EM2C  FORTRAN  program  with  initial  guesses  for  the 
gas  temperature  profile  and  CO2  eoneentration.  The  residual  is  eomputed  based  on  the 
differenee  between  the  experimental  speetrum  and  the  model  predietion.  The  routine 
perturbs  the  values  slightly  to  find  the  largest  gradient  for  the  input  parameters  of 
temperature  and  eoneentration,  and  uses  the  direetion  of  the  gradient  to  iterate  forward  to 
find  the  values  of  eoneentration  and  temperature  that  best  fit  the  experimental  speetrum. 
The  temperature  eorresponding  to  the  best  fit  is  taken  to  be  the  gas  temperature.  The 
resolution  of  the  measurements  is  16  em'^  while  the  resolution  of  the  EM2C  model  is 
only  25  om'\  As  a  result,  the  output  of  the  EM2C  model  must  be  interpolated  onto  the 
experimental  data  before  the  residual  is  eomputed. 
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4,2,  Experimental  Setup: 


Figure  4-6:  Schematic  diagram  of  experiment. 

The  experimental  set  up  consists  of  1)  a  simulated  microcombustor  2)  an  external 
optical  setup  for  non-intrusive  diagnostics  and  3)  mass  flow  controllers  to  control  flow 
rates  and  inlet  gas  composition.  The  entire  set  up  has  been  described  in  detail  else  where 
([16]  and  [82])  and  is  presented  very  briefly  here.  Fig.  4-6  is  a  schematic  diagram  of  the 


127 


overall  experiment.  Two  UFC-3000A  mass  flow  controllers  (MFCs)  from  Unit 
Instruments  control  the  flow  of  fuel  (methane)  and  oxidizer  (air)  into  a  mixing  chamber 
consisting  of  a  1-1/4  inch  diameter  steel  pipe  filled  with  steel  wool.  The  gas  mixture 
exiting  the  chamber  passes  through  a  15  pm  millipore  filter  that  serves  as  a  flash  arrestor 
before  it  enters  the  burner  plenum.  The  MFCs  are  calibrated  during  each  run  using  a 
Drycal  DC-Lite  positive  displacement  flow  meter. 

4.2.1.  Simulated  Microcombustor: 

The  burner  (Fig.  4-7)  itself  consists  of  two  silicon  plates  attached  to  movable  supports 
which  allow  the  gap  between  the  plates  to  be  varied  from  500  pm  to  50  mm. 


Figure  4-7:  Simulated  micro-combustor  mounted  on  a  vertical  traverse. 
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Two  vertical  guides  support  the  plates  on  the  outside  and  a  set  of  4  vertically-oriented 
stainless-steel  tubes  (two  at  each  wafer  edge)  support  the  wafers  from  the  inside.  The 
diameter  of  the  internal  stainless-steel  tubes  (2  mm  in  this  work)  fixes  the  minimum 
wafer  spacing.  Compliant  packings  of  glass  wool  are  held  in  place  along  the  outside 
edges  of  the  plates  by  loose-fitting  ‘C’  shaped  metal  retainers  (see  Fig.  4-7).  Together, 
the  stainless  steel  tubing  and  the  glass  wool  insulation  provide  an  effective  barrier  to 
leakage  and/or  entrainment  from  the  edges  of  the  plate  as  the  flow  passes  vertically 
upward  between  them.  The  silicon  plates  are  constructed  from  ordinary  silicon  wafers 
used  for  MEMS  micro  fabrication.  Hence,  we  refer  to  the  burner  as  a  “simulated” 
microcombustor  as  it  has  many  of  the  important  characteristics  of  a  MEMS  type 
combustor  (material,  passage  dimension)  but  is  made  using  conventional  fabrication 
processes.  Two  types  of  wafers  are  used  in  this  experiment:  a  p-type  wafer  doped  with 
boron  that  has  both  sides  polished,  and  an  n-type  wafer  doped  with  phosphorous  with 
only  a  single  side  polished.  Both  wafer  types  are  described  as  Tightly’  doped  by  the 
manufacturer  but  the  specific  level  of  doping  is  unknown.  The  wafer  electrical  resistivity 
is  10-100  Q-cm  for  both  types  of  wafers,  the  thickness  is  500  pm,  and  the 
crystallographic  orientation  is  100. 
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4.2.2.  External  Optical  System: 


Focusing  Mirrors 
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Figure  4-8:  External  optical  arrangement  (IR  beam  path  in  red). 

As  illustrated  in  Fig.  4-6,  the  beam  from  a  Thermo  Nieolet  Nexus  870  Fourier 
transform  infrared  speetrometer  (FTIR)  is  brought  to  a  foeus  inside  the  simulated 
mieroeombustor  using  a  190  mm  foeal  length  off-axis  parabolie  mirror.  The  off-axis 
mirror  reduees  the  beam  diameter  from  approximately  38  mm  to  10  mm.  Fig.  4-8  is  a 
photograph  of  the  external  optieal  system.  The  IR  beam  from  the  FTIR  speetrometer  is 
illustrated  in  red.  This  beam  passes  through  the  simulated  mieroeombustor,  is  eolleeted 
by  a  parabolie  mirror  with  a  foeal  length  of  50.8  mm,  and  is  foeused  onto  the  FTIR’s 
Mereury  Cadmium  Telluride  (MCTA*)  deteetor  that  has  been  removed  from  the 
instrument  eabinet  and  attaehed  to  the  optieal  beneh.  The  theory  of  operation  of  the 
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spectrometer  is  described  in  detail  elsewhere  [82],  It  suffices  here  to  say  that  the  FTIR 
spectrometer  uses  a  Michelson  arrangement  to  create  an  optical  path  difference  between 
the  two  light  beams  which  are  then  recombined.  The  optical  path  difference  in  the 
recombined  beams  creates  an  interference  pattern  (interferogram)  whose  resolution  is 
controlled  by  the  speed  of  the  moving  mirror.  The  entire  recombined  beam  is  passed 
through  the  sample  (in  this  case  the  microcombustor)  which  absorbs  light  at  different 
wavelengths  according  to  the  type,  concentration,  and  temperature  of  the  species  present. 
The  total  radiant  flux  that  over  all  wavelengths  that  makes  it  through  the  sample  is 
collected  by  a  detector  that  records  the  intensity  as  a  function  of  time.  This  time  series  is 
transformed  to  frequency  (or  wavenumber)  space  by  performing  a  numerical  Fourier 
transform. 


Detector 


Figure  4-9:  Side  view  of  optical  system. 
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As  mentioned  previously,  the  IR  beam  emerging  from  the  spectrometer  is  focused 
down  to  a  beam  of  about  10  mm  diameter  by  the  first  focusing  mirror.  However,  better 
spatial  resolution  is  necessary  to  resolve  temperature  gradients  in  the  streamwise 
direction.  Spatial  resolution  is  improved  in  this  experiment  using  a  spatial  filter  (mask 
with  a  1mm  slit)  as  shown  in  Fig.  4-9.  This  reduces  the  size  of  the  beam  so  that  it  only 
extends  1  mm  in  the  flow  direction  and  5  cm  in  the  spanwise  direction.  The 
microcombustor  is  mounted  on  a  vertical  jack  that  enables  different  streamwise  locations 
in  the  burner  to  be  probed  without  adjusting  the  external  optics. 

4.2.3.  Flow  Control  System: 

Fig.  4-10  shows  the  flow  control  apparatus.  As  discussed  earlier,  it  consists  of  an 
MKS  digital  controller  which  controls  the  MFC’s  which  in  turn  are  calibrated  using 
DryCal  DC-Lite  flow  meters. 


Figure  4-10:  Flow  Control  System 
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The  flow  control  apparatus  enables  control  of  the  mixture  composition  (or 
equivalence  ratio)  and  volumetric  flow  rate  (or  velocity).  When  the  flame  is  stabilized 
inside  the  burner,  the  flame  speed  is  equal  to  the  velocity  of  the  reactants. 

4.3.  Experimental  Procedure: 

The  flow  control  system  (MKS  unit)  is  started  after  opening  the  fuel  cylinder  valve 
and  air  supply  (drawn  from  an  external  tank).  Initially,  only  air  flows  into  the  burner  and 
the  fuel  control  unit  is  set  to  zero.  The  MCT  A*  detector  (connected  externally  to  the 
spectrometer)  is  cooled  by  filling  its  dewar  with  liquid  N2.  The  FT-IR  spectrometer  is 
switched  on  and  the  external  optical  system  is  adjusted  to  yield  maximum  signal  strength. 
A  background  scan  is  obtained  and  stored.  Then,  the  flow  is  seeded  with  a  trace  quantity 
of  CO2  (1-2  %)  and  spectra  are  collected  to  verify  that  the  CO2  spectrum  is  distinctly 
visible.  This  trace  quantity  of  CO2  aids  in  collection  of  good  spectra  in  the  pre-flame 
regions  where  CO2  concentrations  are  fairly  low.  The  CO2  flow  is  maintained  throughout 
the  experiment  for  reasons  that  will  be  discussed  in  the  next  section. 

The  fuel  flow  rate  is  set  to  the  value  associated  with  the  desired  equivalence  ratio  and 
the  premixed  fuel/air  (either  methane  or  propane)  mixture  enters  the  plenum  below  the 
plates.  The  equivalence  ratio  and  the  velocity  of  the  flow  are  maintained  by  the  fuel  and 
air  flow  controllers.  A  stainless  steel  screen  stretched  across  the  discharge  slot  of  the 
plenum  flattens  the  velocity  profile  of  the  flow  entering  the  plate  gap  and  serves  as  a 
flashback  arrestor.  The  mixture  is  ignited  with  the  silicon  plates  at  the  desired  spacing  (~ 
2mm).  The  flame  typically  stabilizes  on  the  outer  lip  of  the  burner  at  the  start  up  owing  to 
the  burner  being  at  room  temperature.  The  silicon  plates  get  heated  by  conduction  from 
the  flame  and  sometimes  are  also  heated  externally  by  a  flame  lighter  or  an  external  heat 
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gun  which  raises  the  temperature  of  the  plates  and  allows  the  flame  to  stabilize  inside  the 
ehannel.  The  external  heating  is  then  stopped  and  the  eombustor  is  allowed  to  equilibrate 
to  its  final  state  with  a  flame  standing  inside  the  ehannel  between  the  two  silieon  plates. 

Onee  the  eombustor  reaehes  equilibrium  (about  20  minutes  after  the  external  heating 
is  stopped)  the  FT-IR  spectrometer  begins  eollecting  absorption  speetra  (using  the 
previously  stored  baekground  sean  as  the  baseline)  at  various  vertieal  positions  in  the 
burner  i.e.  in  the  pre-flame,  flame,  and  post-flame  regions.  Typieally,  200  seans  are 
eolleeted  per  position  whieh  takes  about  2  minutes. 

An  Omega  0.01”  diameter  k-type  ungrounded  thermoeouple  with  a  stainless  steel 
sheath  measures  the  outer  surfaee  temperature  of  the  silieon  plates  at  a  single  point 
through  direet  eontaet.  This  is  used  to  ealibrate  measurements  of  the  surfaee  temperature 
distribution  aeross  the  plates  made  using  a  FLIR  Systems  SC3000  infrared  eamera.  The 
IR  eamera  colleets  an  image  of  emitted  intensity  in  the  8-12um  range  from  the  silicon 
wafer.  The  wafer  surfaee  temperature  distribution  ean  be  obtained  from  this  image  using 
the  Stefan-Boltzmann  equation  for  radiation  [80] 

q"  =  C>  -T^)  (4-5) 

In  this  expression,  q"  is  the  ineident  radiation  flux  per  unit  area,  <j  is  the  Stefan- 

Boltzmann  eonstant  and  s  is  the  emissivity.  The  above  equation  is  embedded  in  the 
software  that  interfaees  the  IR  eamera  to  the  eomputer,  however,  the  software  has  a  built 
in  starting  guess  for  the  emissivity  whieh  will  not  be  aeeurate  for  all  operating  eonditions. 
The  single  point,  direet  eontaet  measurement  made  with  the  thermoeouple  is  used  to 
eorreet  for  this  by  recaleulating  the  emissivity  (using  in  built  funetions  in  the  software). 
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Applying  the  emissivity  correction  to  the  IR  camera  image  gives  the  final  surface 
temperature  distribution. 

4.4.  Temperature  Measurement  in  the  Microburner: 

The  optical  path  between  the  two  silicon  wafers  is  split  into  n  equal  cells.  Each  of  the 
n  cells  is  assigned  a  temperature  based  on  the  kind  of  profile  we  choose  to  fit  (2"‘*  or  4* 
order  profiles)  for  the  temperature  between  the  two  wafers. 


Assumed  Temperature 


Si  Wafer 


Plate  spacing  (W) 

Figure  4-11:  Illustration  of  discretized  beam  path  with  an  assumed  temperature 
distribution. 

Figure  4-11  illustrates  the  optical  path  for  absorption  between  the  two  silicon  wafers 
along  with  a  cartoon  of  the  assumed  temperature  profde.  Previous  work  [82],  [76] 
showed  that  both  2"‘^  and  4*  order  profiles  represent  the  temperature  profile  with 
adequate  accuracy.  The  temperature  gradient  at  the  wall  was  captured  to  within  1%  by 
the  4*  order  profile  and  somewhat  less  accurately  by  the  2"‘*  order  profile.  However,  this 
conclusion  was  based  on  less  sophisticated  modeling  that  assumed  plug  flow  in  the 
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channel.  Subsequent  modeling  efforts  (as  seen  in  ehapter  3)  have  shown  that  under  some 
eonditions,  the  temperature  aetually  peaks  off  the  ehannel  eenterline.  Capturing  this 
behavior  requires  a  4*  order  profile  at  a  minimum  and  henee  this  is  what  we  use  here. 
The  4*  order  temperature  profile  is  given  by 


T(y)  =  (4-6) 

In  this  expression,  y  varies  from  zero  to  2  d,  the  passage  width,  whieh  in  the  eurrent 
experiment  is  around  2.10mm.  The  temperature  profile  is  foreed  to  be  symmetrie  about 
the  ehannel  eenterline  by  requiring  that: 


T(0)  =  T(2d) 


(4-7) 


dT 


dy 


=  0 


y=d 


(4-8) 


dT 


dT 


dy 


y=0 


(4-9) 


Subjeeting  eq.  4-6  to  these  eonditions  gives 

=  -4d  «4  and  d^  =  -2d  C4  +  <34  {Id)^  (4-10) 

Eq.  4-10  relates  two  of  the  fitting  parameters  to  the  other  three.  Therefore,  there  are  only 
three  independent  temperature  fitting  parameters  (namely  a^,  C4  and  64)  even  though  we 
fit  for  a  4*  order  profile. 
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4.5.  Challenges  in  Implementation  and  Remedies: 


4.5. 1.  Interference  from  Silicon  and  CO: 

Fig.  4-12  compares  absorption  spectra  of  CO2  and  CH4,  with  the  transmission 
spectrum  of  silicon  over  the  same  range.  It  shows  that  the  silicon  wafer’s  transmissivity  is 
approximately  60%  except  for  a  small  region  around  the  CO2  band  where  it  varies  from 
55%  to  60%. 


Figure  4-12:  Comparison  of  CO2,  CH4  (absorption)  and  Si  (%  transmittance) 
spectra. 

This  interference  prevents  capture  of  the  strongest  part  of  the  CO2  band  in  the  2300- 
2400  cm'^  range.  The  complete  nature  of  this  interference  is  unclear  but  the  ‘dip’  in  the 
absorption  spectrum  of  CO2  near  the  peak  appears  to  be  related  to  spectral  features  of 
both  silicon  and  CO2.  Silicon  has  an  atomic  absorption  line  [88]  at  -2340  cm'^  but 


137 


experimental  measurements  [89]  of  CO2  absorbance  that  are  not  made  with  silicon  in  the 
optical  path  (see  figures  4-4,  4-8  and  4-9  of  ref  [89])  also  show  a  ‘dip’  in  this  area  -  one 
that  is  not  predicted  by  EM2C.  Regardless,  the  interference  is  a  function  of  temperature 
that  is  not  easily  subtracted  out  and  it  causes  problems  when  trying  to  fit  the  CO2  band. 

The  EM2C  program  which  is  used  to  fit  the  experimental  data  requires  input  of 
species  concentrations  of  H2O,  CO  and  CO2,  the  path  length  over  which  the  absorption 
spectrum  is  to  be  calculated,  and  a  guess  for  the  gas  temperature. 


Figure  4-13:  Absorbance  spectra  of  CO2  showing  tbe  relative  influence  of  each 
species  given  as  an  input  to  tbe  EM2C  program  in  tbe  spectral  region  of  interest  (T 
=  900  K,  Optical  Path  Length  (2  d)  =  2,15  mm) 

The  spectral  (wavenumber)  region  over  which  the  absorption  of  CO2  occurs  without 
interference  from  other  species  can  be  predicted  by  the  EM2C  model.  Fig  4-13  shows  the 
relative  absorbance  of  each  of  the  species.  The  mole  fraction  values  assumed  for  each 
species  correspond  to  a  CHVAir  premixed  flame  burning  at  stoichiometric  proportion 
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calculated  using  CHEMKIN.  The  figure  indicates  that  we  ean  expect  interference  from 
CO,  whieh  is  likely  to  be  present  in  hydroearbon  flames,  in  the  wavenumber  region  of 
2150-2250  em'^  Therefore,  in  order  to  avoid  the  region  where  this  interference  occurs, 
the  experimental  data  is  fitted  in  the  range  of  2250  -  2300  cm'^  wavenumbers.  The 
solution  implemented  here  is  to  fit  only  the  region  of  the  band  (2250  -  2300  em'^)  that  is 
not  affected  by  the  interference  from  silieon  and  CO2  and  from  CO  on  the  other  side  of 
the  band.  This  is  illustrated  in  Fig.  4-14. 


Figure  4-14:  Fit  of  CO2  spectra:  (|)  =  0,85,  U  =  50  cm/s,  2,12  mm  wafer  spacing,  21 
mm  downstream  from  inlet 


4.5.2.  High  Temperature  Effects: 

Apart  from  the  interference  in  the  CO2  band  by  silicon,  additional  complications  arise 
because  the  IR  transmissivity  of  silieon  is  also  a  relatively  strong  funetion  of  temperature 
[90].  This  neeessitates  making  measurements  of  the  temperature  dependence  of  the 
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transmissivity  of  the  silicon  plates.  This  was  accomplished  by  heating  the  plates  with  a 
heat  gun,  recording  the  plate  temperature  with  a  sheathed  k  type  thermocouple,  and 
recording  the  absorbance  measured  by  the  FTIR. 


Figure  4-15:  Percent  transmittance  of  double  side  polished  Si  wafer  vs,  temperature 
at  different  wave  numbers  around  the  CO2  regime  (2250  -  2450  cm'^). 

The  measured  transmittance  of  Si  as  a  function  of  temperature  at  2250  cm'^  and  2450 
cm"'  is  shown  in  Fig.  4-15.  Transmittance  in  the  spectral  region  around  the  CH4  band 
varied  in  a  similar  manner.  The  significant  scatter  in  this  data  arises  primarily  from 
inconsistencies  in  the  thermal  contact  between  the  thermocouple  and  the  silicon  plates  as 
the  thermocouple  was  repeatedly  pressed  against  the  plate  surface.  The  line  drawn 
through  the  data  points  in  this  figure  is  a  best  fit  polynomial.  The  data  show  that  the 
transmittance  of  the  silicon  wafers  drops  dramatically  with  increasing  plate  temperature. 
Alone,  this  does  not  pose  a  problem  for  the  diagnostic  technique  because  a  baseline  shift 
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and  scaling  is  performed  on  every  speetrum  before  it  is  co-added.  However,  the  deerease 
in  transmitted  signal  deereases  the  signal/noise  ratio  of  the  transmitted  intensity 
measurement  and  therefore  increases  the  uneertainty  in  eoneentration  and  temperature 
measurements.  As  a  result,  uneertainties  in  eoneentration  and  temperature  inerease  with 
operating  temperature  beeause  of  the  inereasing  opaeity  of  silieon 

In  the  results  presented  in  Chapter  5,  the  peak  temperature  of  the  silieon  plate  reaches 
over  350  °C  when  operating  at/elose  to  stoiehiometrie  eonditions.  This  eauses  a  dramatie 
drop  in  the  transmissivity,  a  deerease  in  signal  to  noise  ration  and  henee  very  poor  quality 
data.  Therefore,  we  are  only  able  to  obtain  good  data  in  leaner  ((j)  <  0.9)  and  richer  flames 
(cj)  >  1.15),  where  the  silicon  wall  temperature  stays  below  350  “C.  Nevertheless,  we 
would  like  to  be  able  to  make  measurements  over  the  entire  flammability  range  ineluding 
stoiehiometrie  eompositions. 

A  fix  was  attempted  by  replaeing  the  silieon  wafers  with  sapphire  plates  of  the  same 
thieknesses.  Sapphire  has  thermal  and  material  properties  similar  to  silieon  but  has  mueh 
higher  IR  transmissivity  than  silieon  so  we  were  hopeful  that  this  would  solve  the 
problem.  However,  the  sapphire  wafers  turned  out  to  be  mueh  more  sensitive  to  thermal 
gradients  than  silieon  as  they  eraek  almost  immediately  when  they  eome  in  eontact  with  a 
flame.  This  eould  not  be  cireumvented  even  by  increasing  the  plate  thiekness.  The  failure 
of  the  sapphire  windows  foreed  a  return  to  silieon.  Therefore,  a  means  of  limiting  the 
flame  temperature  is  required.  This  eould  be  aecomplished  by  diluting  the  flow  with 
exeess  N2  and  is  diseussed  in  more  detail  in  the  future  work  (Chapter  6). 
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4.5 .3.  Lack  of  CO?  in  pre-flame  region 


The  temperature  measurement  teehnique  relies  on  the  presenee  of  CO2  in  the 
absorption  path. 


Figure  4-16:  Comparison  of  CO2  absorbance  spectrum  at  tbe  same  location  for 
seeded  and  unseeded  flows  in  tbe  pre-flame,  <j)  =  0,8,  U  =  50  cm/s,  z  =  12  mm. 

Near  the  flame  and  in  the  post-flame  regions,  there  is  a  suffieient  eoneentration  of 
CO2  owing  to  the  produetion  of  CO2  in  the  reaetion  zone.  However,  in  the  pre-flame 
region  only  baekground  quantities  (384  ppm)  of  CO2  are  present.  This  reduees 
signal/noise  and  signifieantly  affeets  the  quality  of  the  data.  Fig  4-16  shows  the  pre-flame 
absorption  speetrum  in  the  CO2  absorption  region  for  the  two  oases  (unseeded  and 
seeded).  The  equivalenoe  ratio  is  0.8  and  the  flow  velooity  is  50  om/s  at  a  pre-flame 
looation  of  z  =  12  mm  from  the  slot.  As  oan  be  seen  from  the  figure,  in  the  absenoe  of 
seeding  (solid  line)  the  absorption  speetrum  in  the  pre-flame  region  is  not  differentiable 
from  the  noise.  Attempting  to  fit  this  using  the  EM2C  program  will  result  in  large  error 
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bars  and  an  effectively  meaningless  gas  temperature.  Seeding  the  flow  with  a  trace 
quantity  (1-2%)  of  CO2  by  volume  (dashed  line)  improves  signal  strength  dramatically 
allowing  us  to  make  quantitative  absorption  measurements  in  the  pre-flame.  A  few 
sample  calculations  were  made  in  CHEMKIN  for  pre-mixed  flames  with  a  2%  addition 
of  CO2  to  the  reactant  mixture  and  the  variance  in  flame  speed  (from  unseeded)  was 
found  to  be  less  than  5%. 


4.6.  Gas  Temperature  Uncertainty: 


Several  factors  (instrumental  errors,  temperature  fitting  errors  etc.)  contribute  to  the 
final  uncertainty  of  the  gas  temperature  inferred  using  the  diagnostic  technique  developed 
here.  In  order  to  make  meaningful  measurements,  it  is  very  important  to  quantify  the 
various  sources  of  uncertainty  and  to  understand  how  each  contributes  to  the  final 
uncertainty  of  the  temperature  measurement.  A  rigorous  analysis  of  the  uncertainty  in  the 
predicted  temperature  is  presented  elsewhere  [16]  so  the  focus  here  will  be  on  what  is 
new:  incorporating  the  uncertainty  introduced  by  the  choice  of  the  functional  form  of  the 
temperature  profile.  However,  a  brief  review  of  our  previous  work  will  be  useful.  The 
standard  error  is  given  by: 


/=1 


In  this  expression,  Rsq  is  the  square  of  the  residual  from  the  fitting  procedure,  Aexj  is  the  i* 
point  of  the  experimentally  measured  absorbance  spectrum,  Afitj  is  the  point  of  the  best 
fit  to  the  absorbance  spectrum,  and  N  is  the  total  number  of  points  defining  the  spectrum. 
The  basic  premise  of  the  uncertainty  analysis  is  to  express  the  absorbance  observed  by 
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the  instrument  as  the  sum  of  the  ‘real’  absorbanee  AreaU  (ie.  the  absorbanee  that  one 
would  measure  in  the  absenee  of  instrument  error)  and  the  error  introdueed  by  the 
instrument: 


A  —A  -\- 

'^real ,i  rms 


(4-12) 


A  similar  approaeh  is  taken  with  the  absorbanee  speetrum  fitting  proeedure: 


A  =  A  + 

fit,i  perfect  fit,i 
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(4-13) 


In  this  expression,  Aper/ect  fu.i  is  the  absorbanee  at  point  i  in  the  absenee  of  fitting, 
measurement,  and  instrument  errors  while  the  seeond  term  represents  the  errors  that  arise 
from  uneertainties  in  the  inputs  to  the  fitting  proeedure.  In  this  latter  term, /is  the  EM2C 
funetion  used  to  eompute  the  absorbanee  as  a  funetion  of  gas  temperature  and 
eoneentration  and  SXcoi  is  the  uneertainty  in  the  gas  mole  fraetion.  The  latter  is  treated 
as  a  known  uneertainty  in  this  analysis.  Substituting  eqs.  4-12  and  4-13  into  4-11  leads  to: 
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The  solution,  derived  from  solving  Eq.  (4-14)  [7],  is: 
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i=i 


In  this  expression,  a  and  ji  are  eonstants  that  ean  be  eomputed  using  5rms,f  and  SKcoi- 
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(4-16) 

n 

Pi,j  -  (4-17) 

There  are  two  solutions  to  Eq.  4-11  and  the  larger  is  taken  as  the  uneertainty  in  the 
temperature.  The  effeet  of  the  temperature  profile  ehoiee  is  ineorporated  as  follows:  A 
fourth  order  temperature  profile  is  represented  as: 

T(yj)  =  «4  y]  +  ^4  y]  +  ^4  y]  +  d,  yj  +  ^4  (4-i8) 

where, 7  is  the  mid-point  of  eaeh  of  the  n  eells  and  the  plate  spaeing  is  2d. 

IdQ-V.) 

Q  =  \:n);  y.= - ^  (4-19) 

n 

The  EM2C  funetion  depends  on  the  eoeffieients  of  the  temperature  profile  and  the  CO2 
eoneentration 

f  ~  f  5  ^4  ’  ^4  ’  ^C02  )  (4-20) 


Now,  looking  baek  at  Eq.  4-17,  we  ean  express  fdij  as 
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The  index  i  represents  the  wavenumber  spaee  over  whieh  the  absorbanee  speetrum  is 
eolleeted  and  the  index  j  represents  the  loeation  in  the  ehannel  between  the  two  plates  at 
whieh  the  uneertainty  is  being  eomputed.  Direet  differentiation  of  eq.  4-18  gives: 
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(4-22) 
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dTj 

Substituting  eqs.  4-22,  4-23  and  4-24  in  eq.  4-21  and  then  into  eq.  4-15  gives  the 
uneertainty  in  temperature  at  different  loeations  in  the  ehannel  while  aeeounting  for  the 
ehoiee  of  the  fitted  temperature  profile. 


4.7.  Fundamental  limitations: 

The  basie  diagnostie  teehnique  deseribed  above  was  developed  and  verified 
experimentally  in  earlier  work  [16].  The  prineipal  eontributions  of  this  thesis  in 
diagnosties  area  lie  in  the  ineorporation  of  the  temperature  profile  along  a  line  of  sight. 
While  it  has  been  found  experimentally  that  the  approaeh  deseribed  in  §4.4  works  quite 
well  -  i.e.  it  produees  measurements  that  are  eonsistent  with  eonservation  of  energy  and 
other  observations  of  overall  flame  behavior  -  some  important  fundamental  limitations 
have  also  been  diseovered. 

In  order  to  understand  their  signifieanee,  how  they  arise,  and  how  they  ean  be  dealt 
with,  it  is  useful  to  briefly  step  baek  and  review  the  basie  physies  of  the  problem:  The 
degree  of  light  attenuation  by  absorption  is  given  by  Beer’s  law: 

(4-25) 
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In  this  expression,  Iq  is  the  incident  flux,  /  is  the  transmitted  flux,  s  is  the  molar 
absorptivity  of  the  absorbing  species,  X  is  the  wavelength,  T  is  the  temperature,  C  is  the 
concentration  of  the  absorbing  species  and  L  is  the  optical  path  length.  Absorption 
measurements  at  two  different  wavelengths  are  required  to  determine  both  concentration 
and  temperature  of  a  mixture  of  known  gases.  Equation  4-25  assumes  that  the 
composition  and  temperature  are  uniform  across  the  line  of  sight  but  this  may  not  be  the 
case  in  practical  systems.  To  account  for  this,  one  can  break  the  absorbing  path  into  a 
series  of  adjacent  cells  each  with  its  own  temperature  and  concentration.  The  overall 
absorbance  of  a  series  of  N  such  cells,  is  given  by: 


-w=n 


(4-26) 


At  a  minimum,  measurements  of  EIo  at  2N  wavelengths  are  required  to  solve  this  system. 
The  system  is  given  by: 


(4-27) 


There  are  a  number  of  factors  that  complicate  attempts  to  solve  eq.  4-27  for  7)  and  Q 
-  the  temperature  and  concentration  distributions  along  the  line  of  sight: 

1 .  Since  the  absorbances  appear  as  a  product,  T,  and  C,  can  appear  in  any  order  in  the 
product  without  changing  the  total  observed  absorbance.  This  means  that  spatial 
information  is  lost  when  solving  eq.  4-27  for  T,  and  C,  directly  and  individual  cells 
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can  be  shuffled  and  arranged  in  any  order  to  give  the  same  intensity  ratio  (or 
absorbanee) 

2.  The  2N  wavelengths  required  to  solve  the  system  need  to  be  spaeed  far  enough 
apart  so  that  the  values  of  I/Io  are  substantially  different.  In  the  eurrent 
experiment,  interferenees  from  silieon  and  CO  eombined  with  the  speetral 
resolution  of  the  FTIR  limit  the  resolution  of  C02’s  band  head  to  only  8  points. 
This  means  that  the  temperature  profile  across  the  ehannel  ean  only  be 
represented  using  4  points  -  a  rather  eoarse  profile. 

3.  Inereasing  the  spectral  resolution  of  the  measurement  so  as  to  inerease  the  number 
of  I/Io(A,i)  available  to  the  fitting  routine  may  actually  make  things  worse.  The 
reason  is  that  inereasing  the  resolution  of  the  measurement  would  mean  that  the 
rotational  lines  will  start  to  emerge  in  the  measured  speetra.  These  features  are 
not  part  of  the  EM2C  model  which  is  only  for  the  rotational  band  envelope  and 
the  presenee  of  these  fluetuations  presenee  would  aetually  make  it  more  diffieult 
to  mateh  the  model  to  the  measurement. 

4.  The  absorbanee  is  a  highly  non-linear  funetion  of  wavelength,  temperature,  and 
coneentration.  This  makes  it  ehallenging  to  solve  eq.  4-27  for  Ti  and  C,. 

a.  The  most  eommon  approaeh  to  this  problem  -  the  least-squares  method  -  is 
not  suitable  for  solving  eq.  4-27  for  Ti  and  C,.  The  reason  is  that  the  least- 
squares  method  seeks  Ti  and  C,  that  minimize  the  differenee  between  the 
predieted  and  target  funetion  over  a  user-defined  fitting  interval.  This  has  the 
effeet  of  biasing  the  fitting  procedure  toward  absorbanee  at  wavelengths  with 
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larger  absorbance  when  matching  at  wavelengths  with  smaller  absorbances  is 
equally  important.  What  is  required  is  a  technique  in  which  a  fitting  tolerance 
can  be  specified  for  each  I/Io(A,i)  point. 

b.  Most  non-linear  systems  have  multiple  solutions.  Methodologies  for 
identifying  the  physical  meaningful  solution  have  not  been  developed. 

It  turns  out  that  our  knowledge  of  the  functional  form  of  the  temperature  profile  - 
established  earlier  from  the  modeling  efforts  -  is  the  key  to  mitigating  many  of  these 
problems.  Instead  of  solving  for  Ti  and  Ci  explicitly  at  every  point,  we  assume  that  the 
concentration  is  approximately  uniform  across  the  channel  and  that  the  general  shape  of 
the  temperature  profile  is  well  represented  using  a  polynomial  (something  shown  earlier). 
This  has  several  effects: 

1 .  It  provides  a  mechanism  for  retaining  the  shape  of  the  temperature  profile  across 
the  channel 

2.  It  reduces  the  number  of  unknowns  in  the  problem  thereby  making  it  much  easier 
to  solve  the  system.  In  this  experiment  with  absorbance  measurements  at  8 
discrete  wavelengths,  representing  the  temperature  distribution  using  a  symmetric 
4*  order  polynomial  reduces  the  number  of  unknowns  from  8  to  4.  This  greatly 
increases  computational  efficiency. 

3.  It  provides  a  mechanism  for  excluding  at  least  some  physically  meaningless 
solutions. 
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This  approach  is  tested  in  the  next  section  by  performing  the  least-squares  fitting 
procedure  using  a  variety  of  initial  guesses  for  the  temperature  profiles  in  different  plaees 
in  the  combustor.  It  is  found  that  in  most  cases,  the  resulting  temperature  profile  is 
insensitive  to  the  initial  guess  provided  a  sensible  initial  guess  is  made.  For  example, 
guessing  that  the  temperature  profile  far  upstream  of  the  flame  peaks  at  the  center  of  the 
channel  as  opposed  to  the  wall  is  not  physieally  reasonable  and  will  produee  a  solution 
with  a  similar  non-physical  character. 

4.8.  Initial  Guess  for  T  fitting  parameters: 

As  diseussed  earlier  in  §  4.4  the  optieal  path  between  the  two  silicon  wafers  is 
diseretized  into  n  eells  for  the  purpose  of  eomputing  the  integrated,  “line  of  sight” 
absorbance  using  the  EM2C  model.  A  4*  order  profile  is  ‘guessed’  for  the  temperature 
between  the  two  wafers  based  on  the  results  of  the  analytieal  modeling.  The  parameters 
of  the  profile  along  with  a  concentration  guess  for  CO2  are  given  as  inputs  to  the  EM2C 
program  to  eompute  the  absorbanee  as  a  function  of  the  wavenumber  and  fit  the 
experimental  data.  The  implementation  of  the  actual  fitting  is  as  follows:  First,  a  matlab 
function  was  developed  that  takes  the  starting  guess  for  the  temperature  fit  parameters 
(namely:  a^,  C4  and  for  eq.  4-6)  and  a  guess  for  the  CO2  mole  fraction  as  input  variables 
and  modifies  the  input  files  used  by  the  EM2C  program.  It  then  invokes  the  EM2C 
program  to  compute  the  corresponding  absorption  speetrum.  Second,  a  matlab  script  was 
created  whieh  loads  the  experimental  data  files  onto  the  workspaee  and  then  calls  the 
above  developed  matlab  function  via  an  internal  optimization  routine  (‘Isqcurvefit’) 
based  on  non-linear  least  squares  fitting  and  optimizes  the  fit  parameters  to  give  a  best  fit. 


150 


The  critical  requirement  is  to  be  able  to  provide  a  good  starting  guess  for  the  fit 


parameters. 

In  order  to  check  this,  we  try  fitting  several  different  temperature  profiles  to 
experimental  data  {(/)  =  1.15,  U  =  50  cm/s).  The  initial  profiles  along  with  the  initial 
guesses  for  their  defining  parameters  are  summarized  in  table  4-1.  Both  the  type  and 
magnitude  of  the  profiles  are  varied. 


Set  No 

Profile  type 

a4  (cm'"^-K) 

C4  (cm'"^-K) 

64  (K) 

1 

Flat 

0 

0 

500 

2 

Flat 

0 

0 

1000 

3 

Flat 

0 

0 

1500 

4 

C  shape 

2e4 

2e3 

500 

5 

C  shape 

2e4 

9e3 

500 

6 

C  shape 

2e5 

2e4 

800 

7 

C  shape 

2e5 

5e4 

800 

8 

Inverted  C  shape 

2e4 

0 

500 

9 

Inverted  C  shape 

5e5 

2e3 

800 

10 

Inverted  C  shape 

7e5 

le3 

800 

Table  4-1:  Initial  guess  choices  for  4*''  order  T  profile 


4.8. 1.  Profile  Comparisons 

i)  Flat  Temperature  Profile:  The  first  temperature  profile  fitted  is  a  fiat  profile.  We 
take  eq.  4-6,  set  and  C4  to  zero,  and  guess  three  different  values  for  (500,  1000,  1500 
K  see  table  above).  We  would  expect  the  solutions  to  converge  as  the  profiles  guessed  are 
qualitatively  the  same  and  only  have  different  starting  values.  The  centerline  gas 
temperatures  returned  for  the  converged  solutions  in  each  case  are  plotted  in  figure  4-17. 
The  figure  shows  that  the  solution  is  independent  of  profile  choice  for  all  but  one  point  in 
the  pre-flame.  Since  the  discrepancy  is  within  the  errorbar  for  the  temperature,  it  is  not 
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2 

significant.  The  R  ,  which  is  a  measure  of  the  goodness  of  the  fit,  was  greater  than  0.99 
for  all  these  points^. 


Figure  4-17:  Plot  of  centerline  temperature  as  a  function  of  axial  distance  for  three 
different  initial  guess  values  for  a  flat  temperature  profile  (^=  1,15,  U  =  50  cm/s) 

Figure  4-18  shows  the  corresponding  transverse  profiles  for  a  few  select  axial 
locations.  The  axial  locations  chosen  are  as  follows:  one  in  the  pre-flame,  one  at  the 
flame,  one  in  the  post-flame  close  to  the  flame,  and  one  far  downstream  from  the  flame. 
By  visual  inspection  of  the  spectra  (the  disappearance  of  CH4)  the  flame  location  was 
ascertained  to  be  14  mm  downstream  of  the  inlet  slot. 

The  figure  shows  that  the  profiles  deviate  very  little  from  the  starting  guess  shape.  As 
with  the  centerline  measurements,  the  different  starting  guesses  converge  to  the  same 
final  solution.  What  this  means  is  that  choosing  an  initially  flat  temperature  profile  causes 

^  value  of  1  implies  perfect  fit. 
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the  EM2C  program  to  treat  the  optical  pathway  as  one  block  instead  of  n  discretized 
blocks  and  the  temperature  returned  is  effectively  the  average  gas  temperature  between 
the  two  wafers. 

X  =  10  mm 

o-2r!  ^  n 

0.1-  I  I  ■ 

°  500  1000  1500 

X  =  19  mm 


500  1000  1500 

T(K) - . 

Figure  4-18:  Plot  of  temperature  as  a  function  of  transverse  coordinate  at  four  axial 
locations  for  four  different  initial  “C”  shaped  guesses  (^=  1,15,  U  =  50  cm/s) 

ii)  “C”  shaped  profile:  The  next  profde  we  consider  is  a  “C”  shaped  profde,  in 
which  the  wall  is  hotter  than  the  centerline  of  the  flow  as  might  be  the  case  in  the  pre¬ 
flame.  We  consider  the  four  initial  guesses  outlined  in  table  4-1.  Figure  4-19  shows  the 
corresponding  centerline  and  wall  temperatures.  Unlike  the  flat  temperature  profiles,  they 
do  not  converge  nicely  to  the  same  solution  at  all  axial  locations.  There  are  a  few  features 
to  note  in  these  plots.  First,  the  wall  temperature  is  always  higher  than  the  centerline  at  all 
axial  locations  (mimicking  the  starting  shape).  Second,  initial  guess  sets  5  and  6  overlap 
the  most.  From  the  shapes  of  the  temperature  profiles  evidenced  in  Chapter  3,  it  is  clear 
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that  this  initial  profile  shape  is  a  good  guess  for  the  pre-flame  region  and  not  a  likely 
shape  for  the  post-flame  region.  Accordingly  it  is  seen  that  three  of  the  sets  (4,  5  and  6) 
converge  to  the  same  solution  in  the  pre-flame. 
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Figure  4-19:  Plot  of  centerline  and  wall  temperatures  as  a  function  of  axial  distance 
for  four  different  initial  guess  values  for  a  “C”  shaped  temperature  profile  (^=  1,15, 
U  =  50  cm/s) 

The  very  high  wall  temperatures  seen  at  some  points  (near  the  flame)  is  probably  a 
result  of  the  shuffling  problem  identified  earlier  coupled  with  the  fact  that  near  the  flame, 
the  concave  upward  temperature  profile  is  no  longer  appropriate..  A  plot  of  the  transverse 
profiles  across  the  channel  gives  more  insight.  Figure  4-20  shows  transverse  profiles  at 
the  four  axial  locations  selected  previously  for  the  flat  profile  case. 
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Figure  4-20:  Plot  of  temperature  as  a  function  of  transverse  coordinate  at  four  axial 
locations  for  four  different  initial  “C”  shaped  guesses  =  1,15,  U  =  50  cm/s).  The 
dashed  lines  show  the  guesses  and  the  solid  lines  show  the  converged  solutions. 

The  figure  shows  that  the  shape  of  all  converged  solutions  matches  that  of  the  starting 
guess.  Guess  set  4  has  a  very  slight  curvature  and  the  solutions  (blue  lines)  corresponding 
to  this  guess  also  converges  with  a  slight  curvature  almost  like  the  flat  profile  guess. 
Initial  guess  set  5  (red  lines)  starts  with  the  same  wall  temperature  (500  K,  see  table  4-1) 
as  set  4  but  has  more  curvature  (about  100  K  between  wall  and  centerline).  Guess  sets  5 
and  6  have  the  same  curvature  but  start  with  different  wall  temperatures  (500  K  and  800 
K).  Guess  sets  6  and  7  start  at  the  same  wall  temperature  (800  K)  but  with  set  7  having  a 
strong  curvature. 

The  main  point  where  they  deviate  is  at  the  flame.  The  reason  for  this  could  be  that 
the  choice  of  fitting  using  a  least  squares  approach  cannot  guarantee  a  close  fit  over  all 
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spectral  wavelengths.  In  the  lower  absorbance  regions  (around  wavelengths  of  2250  cm’') 
when  the  least  squares  approach  tries  to  minimize  the  sum  of  the  squares  of  the  residual, 
the  residual  would  be  lower  than  at  a  higher  wavelength  primarily  owing  to  the  low 
value  of  the  absorbance  there.  A  plot  of  the  fits  (Figure  4-21)  at  the  corresponding  axial 
locations  shows  that  the  least  squares  may  not  be  the  best  approach  as  seen  from  the  fit  at 
X  =  14  mm. 


X  =  10  mm 


X  =  14  mm 


<D 


Wavenumber  (cm"^) 


Figure  4-21:  Plot  of  EM2C  fit  to  the  experimental  absorbance  as  a  function  of 
wavenumber  at  four  axial  locations  for  four  different  initial  “C”  shaped  guesses  (^  = 
1.15,  U  =  50  cm/s) 

2  2 

Since  R  is  greater  than  0.99  for  all  of  these  points,  it  would  be  prudent  not  to  use  R 

as  the  only  measure  of  fit  quality.  Rather,  a  convergence  criterion  should  be  imposed  at 
each  point  of  the  fit  to  ensure  a  good  fit. 
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iii)  Inverted  “C”  shape:  Finally  we  test  an  inverted  “C”  profile  in  whieh  the 
eenterline  is  hotter  than  the  wall  as  might  be  the  ease  in  the  post-flame.  Again,  we  get 
nieely  eonverged  solutions  with  R  greater  than  0.99  at  every  axial  point.  Again,  the 
solutions  follow  the  pattern  of  the  starting  guess  and  henee  the  centerline  is  always  hotter 
than  the  wall  (Fig.  4-22). 


Figure  4-22:  Plot  of  centerline  and  wall  temperatures  as  a  function  of  axial  distance 
for  four  different  initial  guess  values  for  a  “C”  shaped  temperature  profile  (^=  1,15, 
U  =  50  cm/s) 
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Figure  4-23:  Plot  of  temperature  as  a  function  of  transverse  coordinate  at  four  axial 
locations  for  four  different  initial  “C”  shaped  guesses  =  1,15,  U  =  50  cm/s).  The 
dashed  lines  show  the  guesses  and  the  solid  lines  show  the  converged  solutions, 

A  plot  of  the  transverse  profiles  (Fig.  4-23)  at  the  four  axial  loeations  shows  that  two 
of  the  guess  sets  (9  and  10)  that  have  a  pronounced  curvature  in  the  starting  shape 
converge  to  the  same  solution  in  the  post-flame  regions  (x  =  19  and  28  mm).  The  guesses 
do  not  converge  at  the  pre-flame  point  (x  =  10mm)  probably  owing  to  the  fact  that  this 
shape  is  a  bad  guess  for  the  pre-flame.  A  plot  of  the  spectral  fits  (Fig.  4-24)  at  these  axial 
locations  reveals  that  all  sets  fit  well  in  the  post- flame.  At  the  flame  (x  =  14  mm), 
convergence  using  the  least  square’s  method  does  not  seem  to  be  the  best  approach.  Set  9 
(red  circles)  has  the  best  fit  at  the  flame,  and  the  corresponding  transverse  temperature 
profile  is  the  closest  to  that  predicted  by  the  analytical  model  in  Chapter  3. 
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Figure  4-24:  Plot  of  EM2C  fit  to  the  experimental  absorbance  as  a  function  of 
wavenumber  at  four  axial  locations  for  four  different  initial  “C”  shaped  guesses  (^  = 
1.15,  U  =  50  cm/s) 


4.8.2.  Initial  guesses  based  on  analytical  model: 

The  preceding  subsection  has  demonstrated  the  need  to  have  good  guesses  for  the 
shapes  and  the  amplitudes  of  the  temperature  profile  in  order  to  produce  realistic  fits. 
Care  has  to  also  taken  to  visually  inspect  the  fits  or  develop  a  criterion  for  each  point  in 
the  fit  to  ensure  a  good  fit  instead  of  the  least  squares  fit  (especially  in  and  around  the 
flame)  criterion  of  R  .  To  accomplish  this,  the  previously  developed  2-D  analytical  model 
for  the  temperature  (Chapter  3)  can  be  used  as  a  basis  for  the  initial  guesses  of  the 
temperature  profile  shapes.  Recall  that  the  model  was  developed  for  two  different 
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velocity  profiles:  Plug  flow  and  fully  developed  Poiseuille  flow.  Schilichting  [ref]  gives 
the  entrance  length  (4),  for  the  development  of  the  velocity  profile,  as 

-^  =  0.16  Re ^  (4-28) 

where,  Red  is  the  Reynolds  number  based  on  channel  half  height.  The  Reynolds  number 
in  our  burner  ranges  between  16  and  40  depending  on  the  choice  of  the  average 
temperature  (at  300  K  and  1000  K)  for  evaluating  fluid  properties  like  viscosity,  density 
etc.  Using  this  value  for  Red  in  Eq.  4-28  gives  us  2.5  <  k/d  <  6.5.  In  comparison,  the 
interrogation  region  for  our  burner  starts  at  10  mm  or  l/d~l0.  Therefore,  we  expect  the 
flow  to  be  hydrodynamically  fully  developed  in  the  region  where  optical  measurements 
are  made.  Hence,  the  Poiseuille  flow  solution  is  selected. 


Figure  4-25:  Transverse  temperature  profiles  in  the  channel  predicted  hy  the 
analytical  model  using  the  Poiseuille  flow  solution,  a)  Pre-flame,  h)  Post-flame, 

Figure  4-25  shows  the  temperature  profiles  predicted  by  the  analytical  model  for  fully 
developed  flow  in  the  burner.  The  figure  shows  that  the  pre-flame  profiles  are  generally 
“C”  shaped  curves  with  a  slight  curvature  except  near  the  flame  where  a  double  peak 
shape  is  present.  These  off-center  peaks  are  the  result  of  the  velocity  profile  as  discussed 
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in  the  modeling  seetion.  The  post-flame  profiles  at  the  same  axial  distanees  from  the 
flame  (Fig.  4-25  b)  reveals  that  the  profiles  are  mostly  inverted  “C”  shapes  exeept  very 
near  the  flame  where  there  is  a  bit  of  flattening  of  the  eurvature  of  the  “C”. 

These  starting  shapes  are  used  as  initial  guesses  for  the  temperature  profiles  in  the  pre 
and  post-flame  regions.  Speeial  attention  is  paid  to  the  eonvergenee  between  EM2C  and 
measured  absorbanee  at  the  individual  points  (i.e.,  at  the  different  wavenumbers)  of  the 
fit  -espeeially  in  the  near  flame  region.  While  the  aeeuraey  of  the  profiles  predicted  by 
the  analytical  model  developed  in  Chapter  3  is  limited  by  the  assumptions  made  in  the 
model  (i.e.,  flame  sheet,  fiat  flame,  constant  fluid  properties  with  temperature  etc),  the 
shapes  of  the  temperature  profiles  from  the  model  still  give  a  better  starting  point  for  the 
guesses  than  arbitrarily  attempting  to  fit  the  profiles. 
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Chapter  5:  Experimental  Results 


In  this  chapter  the  absorbance  data  eollected  from  operating  the  burner  over  a  range 
of  equivalenee  ratios  and  mixture  veloeities  is  proeessed  using  the  EM2C  fitting 
proeedure  outlined  in  the  previous  ehapter  to  obtain  temperature  profiles.  The  main  aim 
of  the  experiment  is  to  measure  temperature  profiles  in  the  streamwise  (along  the  flow) 
and  the  transverse  (along  line  of  sight  between  the  two  plates)  directions  and  from  these 
determine  the  gas-wall  heat  transfer  and  then  infer  the  net  heat  reeireulation  to  the  pre¬ 
flame.  The  measured  heat  reeireulation  will  be  compared  with  the  predietions  of  the 
analytieal  model  developed  in  Chapter  3. 


5.1.  Model  Based  Temperature  Fit: 

In  ehapter  4  it  was  shown  that  the  fitting  routine  needs  a  good  starting  guess  in  order 
to  get  temperature  profiles  that  make  physieal  sense  (i.e.,  the  wall  being  hotter  than  the 
gas  in  the  pre-flame  and  viee-versa  in  the  post-flame).  It  was  also  shown  that  knowledge 
of  the  funetional  form  of  the  temperature  profile  was  neeessary  in  order  to  avoid  the  loss 
of  spatial  information  when  fitting  for  “integrated”  absorbanee.  Table  5-1  summarizes  the 
initial  guesses  for  the  temperature  profiles  in  different  parts  of  the  burner  with  referenee 


to  the  flame  loeation  (x/). 


S.No 

a4  (K/em"^) 

c4  (K/em"^) 

e4  (K) 

^location  (mm) 

1 

2.8e4 

le2 

520 

X  <  Xf  -3 

2 

-13e4 

-6e3 

579 

Xf  -3  <  X  <  Xf 

3 

-7e6 

-5e5 

650 

Xf  <  X  <  Xf  + 1 

4 

2e6 

le4 

600 

Xf +1  <  X  <  Xf+5 

5 

5e5 

le4 

550 

Xf+5  <  X 

Table  5-1:  Initial  guess  for  parameters  to  be  supplied  to  tbe  fitting  routine 
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Figure  5-1  shows  typical  centerline  and  wall  temperatures  for  one  particular  operating 


condition  {(/)=  1.15,  U  =  50  cm/s). 


Figure  5-1:  Centerline  and  wall  temperatures  as  a  function  of  downstream  distance 
(CH4/Air,  (|)  =  1,15,  U  =  50  cm/s) 

Figure  5-2  shows  transverse  temperature  profiles  for  a  few  axial  locations.  The  flame 
location  (x/)  is  15  mm  downstream  from  the  entrance  slot.  The  peak  gas  temperature  of 
aboyt  2000  K  occurs  near  the  centerline  of  the  passage  at  the  flame  location. 
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Figure  5-2:  Transverse  temperature  profiles  across  the  channel  at  select 
downstream  locations  (Flame  location  Xf  =  15  mm) 

Comparing  Figs.  5-1  and  5-2,  it  can  be  seen  that  near  the  flame,  where  the  wafer  is 
the  hottest,  the  flt  based  on  the  model  fails  to  give  solutions  that  are  physically 
meaningful.  The  reason  this  partieular  profile  is  suspeet  is  that  the  wall  temperature  is 
lower  than  the  preceding  point  in  the  pre-flame.  This  cannot  be  correct  and  the  reason 
can  be  inferred  from  figure  5-3  whieh  shows  an  IR  eamera  image  of  the  burner  operating 
at  the  same  condition  as  Figs.  5-1  and  5-2.  The  solid  vertical  rectangle  shows  the  zone  of 
interrogation.  We  are  able  to  make  measurements  starting  only  at  about  10  mm  from  the 
entrance  slot  (base)  of  the  burner. 
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342.8-0 


Figure  5-3:  IR  camera  image  of  the  burner  with  the  flame  stabilized  inside 
(CH4/Air,  ^  =  1.15,  U  =  50  cm/s). 

Returning  to  the  reason  for  poor  fitting  near  the  flame,  the  IR  image  shows  that  the 
hottest  spot  (near  the  flame)  is  nearly  340  °C.  Reeall  from  Fig.  4-13  (in  ehapter  4)  that  as 
the  wafer  temperature  approaehes  350  °C  and  higher  there  is  a  dramatie  reduetion  in  the 
transmissivity  (loss  of  50%  and  higher  as  eompared  to  room  temperature).  Therefore  it  is 
logieal  to  eonelude  that  the  poor  fitting  near  the  flame  probably  oeeurs  due  to  1)  loss  in 
transmissivity  of  the  wafer  and  2)  the  inability  of  the  least  squares  to  eonverge  uniformly 
at  all  wavelengths.  However,  away  from  the  flame  in  the  pre  and  post-flame  regions,  the 
fitting  proeedure  seems  to  work  well  and  we  get  physieally  reasonable  solutions.  Figure 
5-4  shows  the  fits  at  a  few  seleet  axial  loeations  to  illustrate  how  the  least  squares 
approaeh  does  not  give  good  solutions  near  the  flame.  An  attempt  to  fix  this  was  made  by 
eonstrueting  a  starting  guess  from  adjaeent  solutions  but  this  still  resulted  in  poor 
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matches  to  the  absorption  speetra.  A  possible  solution  to  this  problem  is  proposed  in 
Future  Work  (Chapter  6). 

X  =  10  mm  X  =  14  mm 


Figure  5-4:  Comparison  of  the  fit  to  the  experimental  data  at  different  axial 
locations 

The  problem  with  high  wafer  temperature  redueing  transmissivity  beeomes  more 
pronouneed  in  near-stoichiometrie  flames  (cj)  ~  1).  To  see  this,  eonsider  IR  images  (Fig. 
5-5)  of  the  combustor  operating  at  (j)  =  1. 
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a) 


b) 


c) 

Figure  5-5:  IR  images  of  the  burner  at  (j) 


d) 

1,0  for  four  different  burning  velocities: 


a)  U  =  55  cm/s;  b)  U  =  60  cm/s;  c)  U  =  65  cm/s  and  d)  U  =  70  cm/s 


The  figure  shows  that  the  wafer  temperature  exeeeds  400  °C  over  an  appreeiable 
portion  of  the  interrogation  region.  To  verify  that  transmission  is  indeed  poor  in  this 
region,  absorbanee  is  plotted  as  a  function  of  wavenumber  in  Figure  5-6  along  with  an 
attempted  fit  made  using  EM2C.  In  the  figure,  (j)  =  1  at  U  =  60  cm/s.  It  is  clear  that  the 
CO2  band  head  is  completely  suppressed  leaving  only  noise.  As  a  corollary,  the  plot  also 
illustrates  that  using  poor  quality  data  results  in  a  failure  of  the  EM2C  function  to  be  able 
to  fit  the  absorbance  data. 
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Figure  5-6:  Absorbance  vs  Wavenumber  ((j)  =  1,0,  U  =60  cm/s) 

While  the  relative  opaeity  of  Silieon  at  high  temperatures  in  the  2200-2300  em"' 
range  makes  it  difficult  to  get  accurate  line  of  sight  profiles  in  the  near  flame  (and  some 
post-flame)  points,  the  diagnostic  technique  can  still  be  used  to  make  some  first  order 
estimates  of  heat  recirculation  (as  will  be  described  later  in  §5.2)  based  on  the  average 
temperature  (which  is  a  function  of  the  integrated  line  of  sight  absorbance  only)  of  the 
gas  in  the  channel  and  the  pre-flame  profiles  which  have  converged  correctly. 

Figure  5-7  shows  a  2-D  temperature  distribution  in  the  pre-flame  and  post-flame 
regions.  The  power  of  the  technique  despite  some  of  its  short  comings  is  that  enables 
measurement  of  temperature  in  both  axial  and  transverse  directions  in  a  very  small 
volume. 
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Figure  5-7:  2D  temperature  distribution  in  the  channel  ((j)  =  1.15,  U  =  50  cm/s) 
a)  pre-flame  region;  b)  Post-flame  region 
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5.2.  Comparisons  of  Temperature  Profiles: 


The  burner  was  operated  over  a  range  of  equivalenee  ratios  and  burning  veloeities  for 
Methane/Air  mixtures  and  a  fewer  number  of  Propane/Air  mixtures.  Table  5-2  shows  the 
matrix  of  the  operating  eonditions  explored. 


S.No. 

Fuel/Air 

Eq.  Ratio  (^z5) 

U  (cm/s) 

1 

Methane 

0.80 

40 

2 

Methane 

0.80 

45 

3 

Methane 

0.80 

50 

5 

Methane 

0.85 

40 

6 

Methane 

0.85 

45 

7 

Methane 

0.85 

50 

8 

Methane 

0.85 

55 

9 

Methane 

0.90 

55 

10 

Methane 

0.90 

60 

11 

Methane 

0.90 

65 

12 

Methane 

0.90 

70 

13 

Methane 

1.0 

60 

14 

Methane 

1.15 

50 

15 

Methane 

1.15 

55 

16 

Methane 

1.15 

60 

17 

Propane 

0.80 

50 

18 

Propane 

0.80 

55 

19 

Propane 

0.85 

55 

20 

Propane 

0.85 

60 

21 

Propane 

0.90 

65 

Table  5-2:  Operating  conditions  for  data  collection  from  the  burner 

Collection  of  reliable  data  was  not  feasible  in  case  13  for  several  reasons;!)  The  high 
flow  rates  needed  to  push  the  flame  into  the  zone  of  interrogation  increase  the  heat 
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release  rate  driving  the  wall  temperature  above  350  °C  as  ean  be  seen  from  Fig.  5-7  b-d. 
This  eauses  a  dramatie  drop  in  transmissivity  (refer  Fig.  4-13)  whieh  drives  the 
absorption  signal  into  the  noise.  2)  At  lower  flow  rates,  (see  Fig.  5-7  a)  the  flame  is  not 
far  into  the  ehannel  to  make  any  reliable  gas  temperature  measurements  to  estimate  heat 
reeireulation. 

For  the  purposes  of  providing  a  simple  comparison  between  operating  conditions,  we 
define  the  average  gas  temperature  at  a  given  axial  location  as 

Tavg  (•^)  =  \l'  U iy)  T {x, y)  dy  (5-i) 


In  this  expression,  d  is  the  channel  half  height  and  U(y)  is  the  fully  developed  Poiseuille 
flow  velocity  profile  given  by 


U(y)  =  U  {2^ 
d 


yd  j 


); 


0<  y<2d 


(5-2) 


Earlier,  in  §  4.8,  it  was  shown  that  the  velocity  profde  will  always  be  fully  developed 
in  this  burner  over  the  region  of  interrogation.  Therefore  the  use  of  this  velocity  profile 
appears  justified 


5.2.1.  Constant  Equivalence  Ratio  dr. 

i)  Methane  Plots: 

The  first  comparison  we  make  is  between  different  flow  velocities  at  the  same 
equivalence  ratio. 
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a 


b 


^  d 

Figure  5-8:  IR  images  of  burner  at  four  different  flow  velocities  (CH4/Air;  (j)  =  0,85) 

a)  U  =  40  cm/s;  b)  U  =  45  cm/s;  c)  U  =  50  cm/s  and  d)  U  =  55  cm/s 


Figure  5-8  shows  IR  images  of  the  burner  at  four  different  flow  veloeities  at  the  same 
equivalenee  ratio  (cj)  =  0.85)  for  methane/air.  As  the  flow  veloeity  increases,  the  outer 
wafer  temperature  increases  and  the  flame  gets  pushed  further  up  into  the  channel.  Notice 
also  that  increasing  the  flow  rate  increases  the  fraction  of  the  total  wafer  area  that  exceeds 
350  °C  (beyond  which  we  have  poor  transmissivity).  This  is  because  increasing  the  flow 
rate  increases  the  total  heat  release.  Figure  5-9  shows  the  span-averaged  gas  temperature 
and  outer  surface  temperature  as  a  function  of  axial  position  in  the  burner  for  two 
velocities  (U  =  40  cm/s  and  U  =  50  cm/s).  As  expected,  the  surface  and  gas  temperatures 
increase  as  the  flow  rate  increases.  Increasing  the  flow  rate  also  causes  the  flame  to  move 
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farther  downstream.  This  is  beeause  more  physieal  spaee  is  required  to  preheat  the 
mixture  to  the  ignition  temperature. 


Figure  5-9:  Average  gas  and  surface  temperatures  as  a  function  of  x  for  two 
different  velocities  for  CHVAir  flame;  (j)  =  0,85 

ii)  Propane  Plots: 

We  make  eomparisons  of  propane  flames  at  the  same  equivalenee  ratio  ((j)  =  0.85)  for 
two  flow  veloeities  of  55  and  60  em/s.  Fig  5-10  shows  the  IR  images  of  the  burner 
operating  at  the  two  veloeities.  Propane  flames  are  more  stable  than  methane  flames  of 
the  same  (j)  and  U.  This  ean  be  seen  by  eomparing  the  relative  distanee  between  the  flame 
and  the  entranee  of  the  ehannel  in  Figs.  5-8  d  and  5-10  a.  The  eonsequenee  of  this  is  that 
as  we  inerease  the  flow  rate,  the  propane  the  flames  get  even  hotter  than  methane  (as  seen 
in  Fig.  5-10  b)  and  again  we  get  poor  data.  This  limits  the  range  of  exploration  for 
propane  flames  more  severely  than  for  methane  flames. 


173 


a)  b) 

Figure  5-10:  IR  images  of  burner  at  two  different  burning  velocities  (CsHg/Air;  (j)  = 

0,85)  a)  U  =  55  cm/s;  b)  U  =  60  cm/s 

A  plot  of  the  average  and  outer  wall  temperature  reveals  trends  similar  to  methane/air 
flames.  As  the  flow  rate  increases,  the  gas  and  wall  temperatures  increase  and  the 
temperatures  rise  in  the  post-flame  region.  The  pre-flame  points  also  have  poorer 
convergence  and  this  leads  to  some  ‘wrinkles’  in  the  axial  temperature  profile  in  the  pre¬ 
flame.  This  is  primarily  because  of  the  higher  overall  wafer  temperatures  with  propane 
flames.  Absorbance  measurements  in  the  pre-flame,  where  the  CO2  concentrations  are 
much  lower  than  the  post-flame  (~1%  vs.  6%),  are  impacted  to  a  greater  extent  by  the 
loss  in  transmissivity. 
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Figure  5-11:  Average  gas  and  surface  temperatures  as  a  function  of  x  for  two 
different  velocities  for  CsHg/Air  flame;  (|)  =  0,85 

5.3.  Heat  Recirculation  and  Flame  Speed: 

The  main  aim  of  the  temperature  measurements  is  to  use  the  profiles  to  estimate  the 
total  heat  reeirculation  that  occurs  in  this  burner.  Fig  5-12  is  a  sketch  of  the  axial 
temperature  profile  showing  the  zone  of  optical  interrogation  and  the  heat  fluxes  that  we 
are  interested  in.  It  shows  that  our  optical  technique  does  not  span  the  entire  space  over 
which  heat  recirculation  occurs  although  it  does  span  most  of  it.  While  the  temperature  at 
the  start  of  the  zone  of  interrogation  {Ti„)  is  already  higher  than  that  of  the  environment,  it 
is  still  possible  to  make  estimates  of  the  heat  recirculation  based  on  what  is  measured. 
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Figure  5-12:  Cartoon  sketch  of  heat  fluxes  used  in  calculation  of  heat  recirculation 
in  the  burner 

This  is  accomplished  as  follows:  An  an  energy  balance  on  the  flow  before  it  reaches  the 
region  of  interrogation  gives: 

=  !hCp(T,„-Tj  (5.3) 

In  this  expression,  Tin  is  the  temperature  as  the  flow  enters  the  optical  measurement 
region  and  Too  is  the  temperature  of  the  gases  entering  the  apparatus  (standard  conditions, 
Too  =  298  K).  We  get  the  second  term  on  the  left  side  of  Eq.  5-3  and  the  right  side  of  Eq. 
5-3  from  the  measured  temperature  profiles.  Therefore,  the  pre-heating  can  be 
caleulated.  From  Chapter  3,  we  recall  that  total  heat  reeirculation  is  defined  as  the  sum  of 


O  +k  A 

^ pre-heating  g  '^channel 


dx 
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the  heat  reeireulation  oceurring  via  gas  and  wall  respectively.  The  heat  recirculation  via 
the  wall  for  our  burner  is  given  by: 


zj  =  Ci  A- 

Re  c  ire ,  wall  ^pre—h  eating 


C Aflame 

J,,: 


dT 


dy 


dx 


wall 


(5-4) 


And,  the  heat  recirculation  via  the  gas  is  given  by: 


H  =k  A 

Kecirc,wall  g  '^channel 


dT 


dx 


flame 


(5-5) 


From  equations  5-4  and  5-5,  the  total  heat  recirculation  becomes: 


H  =  H  +  H 

Re circ,  total  Recirc,wall  Recirc,gas 


(5-6) 


The  total  heat  recirculation  is  calculated  for  each  operating  condition  of  the  methane/air 
flames  and  plotted  below  in  figure  5-13.  The  red  hashed  region  is  the  region  where  data 
could  not  be  obtained  owing  to  excessive  heating  of  the  wafer  and  subsequent  loss  of 
transmitted  signal. 


Figure  5-13:  Heat  Recirculation  as  a  function  of  (|)  and  U  for  CH^Air 
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The  solid  curve  on  the  floor  of  the  plot  marks  the  burner’s  blow-off  limit  and  the  dashed 
marks  the  flash-back  limit.  We  non-dimensionalize  the  total  heat  recirculation  by  the  heat 
recirculation  that  would  occur  in  a  freely  propagating  flame  (like  in  Chapter  3)  via 
condution  through  the  gas.  This  is  given  by: 


H 


* 

Re  circ 


H 


Ks  circ, total 


pS,.r„Cp{T„,  -TJd 


(5-7) 


This  definition  enables  us  to  plot  the  non-dimensional  heat  recirculation  as  a  function  of 
the  non-dimensional  flame  speed  for  the  experiment  and  compare  to  the  predictions  of 
the  analytical  model.  In  order  to  do  this,  however,  recall  that  the  analytical  model 
contains  Pcfree  that  needs  to  be  selected  in  order  to  make  the  comparison.  Pe/ree  depends 
upon  the  freely  propagating  flame  speed  and  a  characteristic  gas  temperature  (which  sets 
the  density,  thermal  conductivity,  etc.).  The  flame  speed  and  the  characteristic  gas 
temperature  should  be  functions  of  equivalence  ratio.  Table  5-3  gives  Si/ree,  the 
corresponding  adiabatic  flame  temperature  (calculated  using  CHEMKIN)  and  Pe/ree  for 
freely  propagating  flames  for  several  characteristic  gas  temperatures  over  the  range  of 
equivalence  ratios  used  for  data  collection  in  the  burner. 


S.No 

4' 

Sufree 

(cm/s) 

Tad 

(K) 

P  e,free 

(T=  300  K) 

P  e.free 

(T  =  1200  K) 

P  e,free 

(T  =  1500  K) 

1 

0.8 

29.33 

1863 

12.63 

4.00 

3.49 

2 

0.85 

33.53 

1917 

14.44 

4.58 

3.99 

3 

0.9 

35.64 

1955 

15.35 

4.87 

4.25 

4 

1.15 

40.71 

2021 

17.54 

5.56 

4.85 

5 

1.2 

38.84 

2011 

16.73 

5.30 

4.63 

Table  5-3:  Parameters  for  freely  propagating  conditions 
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The  table  shows  that  Pe/ree  depends  more  strongly  on  the  temperature  used  to  evaluate  it 
than  the  equivalenee  ratio. 


Figure  5-14:  Comparison  between  model  and  experiments  of  non-D  heat 
recirculation  as  a  function  of  S* l 

Figure  5-14  eompares  the  non-dimensional  flame  speed  and  heat  recirculation 
measured  in  the  burner  to  the  model  predictions  for  a  few  Pefree  values.  The  figure 
experimental  data  are  clustered  together  in  a  small  range  of  non-dimensional  flame 
speeds  between  1  and  2  which  makes  it  somewhat  difficult  to  draw  broad  conclusions  but 
it  appears  that  a  Pefree  of  around  3  correlates  most  of  the  points  within  the  experimental 
uncertainty.  While  the  experimental  uncertainty  is  relatively  large,  the  measurements 
confirm  the  fundamental  prediction  of  the  model  which  is  that  flame  speed  is 
proportional  to  the  degree  of  heat  recirculation. 
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Chapter  6:  Findings  and  Future  Work 


6.1.  Key  Contributions: 

This  dissertation  has  investigated  the  influenee  of  struetural  heat  eonduetion  on  flame 
speed  and  stabilization  in  mieroseale  eombustors.  First,  analytieal  models  were  developed 
for  the  heat  transfer  problem  of  a  laminar  flame  stabilized  in  a  ehannel  with  inereasing 
degrees  of  eomplexity  that  started  with  a  eonstant  temperature  wall,  added  eonjugate  heat 
transfer,  and  finally  added  the  speeies  transport  equations  to  the  eonjugate  heat  transfer 
model.  In  this  last  step,  the  flame  speed  beeame  an  eigenvalue  eontrolled  by  the  design 
parameters  of  the  eombustor.  Seeond,  a  non-intrusive  diagnostie  teehnique  based  on  FT- 
IR  speetroseopy  was  improved  to  enable  measurement  of  temperature  profiles  along  line 
of  sight.  This  new  tool  was  used  to  obtain  2-D  gas  temperature  profiles  inside  a 
simulated  mieroeombustor  that,  in  turn,  were  used  to  validate  the  model’s  predietions. 
The  key  findings  and  eontributions  that  have  resulted  from  eombined  analytieal  and 
experimental  approaeh  are  summarized  below. 

1.  Three  new  analytieal  models  of  inereasing  eomplexity  were  developed  for  the 
stabilization  of  a  flame  in  a  miero-ehannel.  In  all  eases,  the  flame  was  represented  as 
an  infinitesimal  sheet  of  heat  release  and  both  plug  and  Poiseuille  flows  were 
eonsidered.  The  models  in  inereasing  order  of  eomplexity  are: 

a.  Constant  wall  temperature,  no  speeies  equation 

b.  Conjugate  heat  transfer  and  heat  loss  to  environment,  no  speeies  equation 
e.  Conjugate  heat  transfer,  heat  loss  to  environment,  and  speeies  equation. 
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2.  The  analytical  models  confirmed  several  aspects  of  micro-scale  combustion  that  were 
known  already.  First,  heat  recirculation  from  the  post- flame  to  the  pre-flame  through 
the  structure  increases  burning  velocity  beyond  the  freely  propagating  flame  speed 
and  enhances  flame  stabilization.  Second,  the  flame  will  quench  if  too  much  heat  is 
lost  to  the  environment.  Third,  multiple  flame  speeds  are  possible  but  only  the 
solutions  corresponding  to  enhanced  flame  speeds  are  stable. 

3.  The  model  has  extended  our  understanding  of  micro-scale  combustion  in  several 
important  ways.  First,  the  three  principal  design  parameters  governing  flame 
stabilization  in  micro-combustors  are  identified:  The  structure  to  gas  thermal 
conductivity  ratio  {k),  the  flow  passage  half  height  to  wall  thickness  ratio  (r),  and  the 
Nusselt  number  for  heat  loss  to  the  environment  {Nu^.  Second,  the  net  heat 
recirculation  from  the  post-flame  to  the  pre-flame  is  the  main  physical  process 
responsible  for  flame  stabilization.  Changes  in  the  three  design  parameters  and  even 
the  fuel  are  important  only  to  the  extent  that  they  alter  the  net  heat  recirculation. 
Changing  the  design  in  a  way  that  increases  heat  recirculation  increases  flame  speed. 
Third,  the  model  shows  that  the  velocity  profile  influences  the  shape  of  the 
temperature  distribution  across  the  channel.  Near  the  flame,  the  competition  between 
the  convective  and  diffusive  terms  causes  the  peak  temperature  to  occur  on  either  side 
of  the  passage  center  line.  While  previous  work  attributes  these  off-center  peaks  in 
temperature  to  flow  straining,  the  models  developed  here  indicate  that  their  origin  is 
at  least  partly  thermal. 

4.  The  analytical  model  was  used  to  show  that  heat  recirculation  enables  micro-scale 
chemical  thrusters  to  achieve  higher  thrust/weight  than  conventional  thrusters. 
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5.  Transverse  temperature  profiles  in  micro-combustors  are  well  represented  using 
power  laws. 

6.  Key  improvements  were  made  to  a  novel  non- intrusive  diagnostic  technique  for 
measuring  temperature  and  chemical  composition  in  micro-combustors.  The 
technique  was  extended  to  permit  measurement  of  temperature  along  a  line  of  sight 
and  a  comprehensive  uncertainty  analysis  was  presented. 

7.  The  key  limitations  of  this  non- intrusive  technique  were  also  identified.  These  include 
the  need  for  a  good  initial  estimate  of  the  temperature  profile  and  the  un-suitability  of 
the  least-squares  approach  to  spectral  fitting. 

8.  The  diagnostic  technique  was  used  to  measure  net  heat  recirculation  in  a  simulated 
micro-combustor.  These  measurements  confirm  the  model  prediction  that  the  flame 
speed  is  a  linear  function  of  heat  recirculation. 

6.2.  Future  Work: 

This  dissertation  has  advanced  the  understanding  of  the  main  processes  responsible 
for  flame  stabilization  in  narrow  channels  and  has  successfully  shown  the  impact  of  heat 
recirculation  through  the  wall  of  the  combustor  on  the  flame  speed.  However,  much  more 
work  is  needed  to  suitably  model  and  make  measurements  of  the  combustion  process  in 
small  scale  combustors. 

First,  the  current  analytical  model  applies  to  an  infinitely  long  channel  whereas,  real 
combustors  have  finite  lengths.  Therefore,  a  new  model  incorporating  finite  length 
effects  would  be  more  capable  of  predicting  the  flame  speeds  more  accurately.  The 
temperature  profiles  from  the  model  could  also  serve  as  a  guide  for  fitting  routines  that 
make  line  of  sight  measurements. 
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Second,  the  techniques  used  to  fit  the  spectra  have  been  shown  to  be  deficient.  A  new 
fitting  procedure  based  on  a  genetic  algorithm  that  would  fit  each  and  every  wavenumber 
point  may  reduce  or  eliminate  these  problems. 

Third,  the  one-step  reaction  model  in  the  current  work  could  be  replaced  with  a  multi- 
step  model.  While  it  would  result  in  more  computational  expense,  it  would  give  more 
detailed  information  on  the  combustion  process  in  the  narrow  channels  and  could  serve  as 
a  tool  to  verify  the  accuracy  of  using  one  step  models. 

Fourth  and  finally,  the  utility  of  the  diagnostic  technique  would  be  greatly  improved 
by  replacing  the  current  system  of  mirrors  and  masks  with  optical  fibers. 


183 


Appendix  A:  Solution  and  Eigenequation 


The  solution  procedure  involves  separating  variables  in  the  two  coordinate  directions 
andy'^.  Thus, 

=  (A-1) 


Inserting  this  into  the  energy  equation  for  the  gas  leads  to 


-{X'-PeU\y^)X')  Y'' 


T 


(A-2) 


At  first  glance  it  appears  that  separation  of  variables  has  not  worked!  However,  taking 
note  of  the  fact  that  most  combustion/heat  transfer  problems  in  pipe  or  channel  flows 
have  an  exponential  solution  for  the  variable,  we  make  the  following  substitution: 


=  ±a 


(A-3) 


The  choice  of  the  plus  or  minus  is  based  on  whether  we  are  interested  in  pre-flame  or 
post-flame  solutions.  In  all  cases  the  non-dimensional  temperature  must  become  zero  at 
the  extremities.  Substituting  eq.  A-3  into  eq.  A-2  and  rewriting  gives 

Yl+[a^+aPeU\y^)y^=d  (A-4) 

Eq.  A-4  is  an  ODE  that  can  be  solved  after  a  transformation  of  variables.  It  has  a  general 
solution  of  the  form 

Yg{a,y^)  =  [cY^^{a,y^)  +  DY^^{a,y^))  (A-5) 

where 

Y^^{y^)  =  p[ypergeometric\Y\{a,y^)  or  Cos{a,y^)  (A-6a) 

Y^  2{y  )-HermiteH{a,y^)e^^  or  Sin{a,y^)  (A-6b) 
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The  choice  of  the  solution  depends  on  whether  we  solve  for  the  case  of  Hagen-Poiseuille 
flow  or  slug  flow. 

Next,  we  consider  the  boundary  and  interface  conditions  along  the  transverse  (y^) 
direction  in  order  to  relate  the  gas  and  structure  solutions.  Applying  eq.  2-6b  (symmetry 
condition)  to  eq.  A-1  gives: 

f;(0)  =  o  (a-10) 

This  gives 

Y^{0)  =  (CY^,ia,y^)  +  DY;^,{a,y^))=0  (A-11) 

and  leads  to 

Z)}^2fe0)=0  or  D=0  (A- 12) 

since 

T;i(a,0)  =  0  and  Y^^(a,0)^0  (A-13) 

The  Ygj  function,  fails  to  satisfy  eq.  A-10  and  hence  it  is  not  admissible  as  a  solution. 
This  is  a  property  of  that  particular  function  in  itself  namely  the  Sine  or  the  HermiteH 
function.  Therefore,  the  Ygj  (which  satisfies  eq.  A-10)  is  the  only  solution,  either  the 
Hypergeometric  or  Cosine  function.  At  this  point  we  do  not  know  what  these 
eigenvalues  (cfs)  are  but  we  explain  how  these  are  obtained  shortly. 

Isothermal  Wall  (Constant  Wall  Temperature): 

For  the  case  of  constant  wall  temperature  demanding  that  the  Yi  function  satisfy  eq.  2- 
11c  (the  other  boundary  condition  at  the  wall)  allows  us  to  obtain  the  eigenvalues  a„,  i.e. 

7,(a„,l)  =  0  (A-14) 
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Obtaining  the  eigenvalues  (a„)  in  the  ease  of  Plug  flow  is  easier  as  the  roots  of  the  Cosine 
funetion  are  well  known.  However,  in  the  ease  of  Hagen-Poiseuille  flow  the  eigenvalues 
(a„)  have  to  be  determined  numerieally  beeause  the  roots  of  the  HypergeometriclFl 
funetion  do  not  have  a  elosed  form. 

Conducting  Walls  (Conjugate  Heat  Transfer): 

For  the  ease  of  eonjugate  heat  transfer,  we  eannot  obtain  the  eigenvalues  at  the  gas-wall 
interfaee  as  we  allow  for  heat  transfer  through  the  wall.  Before,  we  attempt  to  solve  for 
the  energy  equation  of  the  wall  we  would  like  to  assert  that  the  X(x")  funetions  for  the  gas 
and  strueture  solutions  are  linearly  dependent  on  eaeh  other.  This  is  neeessary  to  mateh 
strueture  temperature  and  heat  flux  at  the  wall.  Therefore  for  simplieity  (and  without  loss 
of  generality)  we  assume  that  \hQX(x")  is  the  same  for  both  gas  and  strueture  i.e. 

X^{x^)  =  XXx^)  (A-15) 

Henee,  when  we  separate  variables  for  the  strueture’s  energy  equation  we  write  it  as 

^  +  ^  =  0  or  Y+aX=0  (A-16) 

Eq.  A-16  has  solutions  of  the  form 

Y^{a,y^)  =  EY^,{a,y^)  +  FY^,{a,y^)  (A- 17  a) 

where 

Y^^=Cos(ay^)  (A-17b) 

Tj  2  =  Sin{a  y^)  (A- 17  e) 

The  temperature  and  heat  flux  mateh  at  the  gas-strueture  interfaee  gives  (keeping  eq.  A- 
15  in  mind): 
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Y^ia,l)  =  Y^ia,l)  (A- 18  a) 

or 

CY^  ,{a,l)  =  EY^,{a,l)  +  FY^,{a,l)  (A-18  b) 

and 

Y^{a,l)  =  icY:{a,l)  (A-18c) 


or 

C 7;, (a,l)  =  ic(e 7;, (a,l)  +  E 7'^ (a,l))  (A-18  d) 


Solving  eqs.  A-18  b  and  d  as  a  system  of  linear  equation  relates  the  eonstants  E  and  F  to 
C.  The  link  between  the  temperature  fields  in  the  strueture  and  gas  arises  through  the 
relationship  of  E  and  F  to  C  and  also  beeause  the  X(x")  funetions  are  matehed.  Finally, 
we  invoke  the  last  boundary  eondition  (eq.  2-28d)  for  the  heat  loss  to  the  environment 
from  the  outer  surfaee  of  the  wall.  The  resulting  equation  is  termed  the  “eigenvalue 
equation” 


-7;(«,r)  = 


Nuj,  7^(a,r) 


(A-19) 


The  only  unknown  in  eq.  A-19  is  the  eigenvalue  (a).  The  linear  relationship  between  the 
eonstants  E  and  F  in  the  strueture ’s  solution  to  the  eonstant  C  of  the  gas  phase  solution 
ensures  that  the  above  equation  (after  rearrangement)  is  a  multiple  of  a  funetion  of  a  and 
C.  The  eigenvalues  of  the  problem  are  all  values  of  a  that  satisfy  eq.  A-19  and  do  not 
eorrespond  to  C=0,  (the  trivial  solution  to  the  PDF  where  the  temperature  is  zero 
everywhere).  Eaeh  eigenvalue  has  a  eorresponding  eigenfunetion  (T)  eaeh  of  whieh 
satisfies  the  ODE.  Consequently,  the  linear  eombination  of  all  these  funetions  is  the  final 
solution  to  the  PDE. 
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Appendix  B:  Orthogonality  Condition 


In  order  to  obtain  an  orthogonality  condition  we  go  back  to  the  ODE  eqn  for  Y  and  write 
it  out  for  two  different  eigenvalues  as  follows. 

C  +(«»  +“..  =0  (B-l  “) 

r,".  +(«;  =0  (B-i  b) 

multiply  the  first  equation  by  Y„  and  second  equation  by  Y^  and  subtract  one  from  the 
other  and  integrate  between  0  and  1  to  obtain,  i.e.  jeq.B-laxY  -eq.B-lbxY  dy* 


1  1 


(B-2) 


The  first  integral  can  be  integrated  by  parts.  This  gives 


1 

(r  Y  -Y'  Y  T+(a  -a){(a  +a  +PeUUy^))Y  Y  dy" 

V  ^5^  g^^^  \  fn  n/j\  m  n  Vy  //  g^m  g,n  p 


=  0 


(B-3) 


0 

Now  we  can  invoke  the  boundary  condition  at  the  centerline  (eq.  A- 10)  to  further 
simplify  eq.  B-3  as 


7'  Y 

g,n  g,m 


+«„  +PeU 


(B-t) 


Isothermal  Wall  (Constant  Wall  Temperature) 

In  the  simple  case  of  a  constant  wall  temperature  the  first  term  in  eq.  B-4  vanishes  owing 
to  the  boundary  condition  (eq.  A- 14  in  Appendix  A).  This  gives  the  orthogonality 
condition  as 
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i(a,  +a.  +PeU'(y'))Y^J,,dy'=0  (B-5) 

0 


Conjugate  Heat  Transfer: 

Revisiting  eq.  B-4,  we  can  use  the  relationship  between  the  Y  functions  of  the  structure 
and  gas  owing  to  the  interface  conditions  eqs.  A-18  a  and  c  and  rewrite  eq.  B-4  as 


1 

k{y'  Y  -Y'  Y  1  +(a  -a){ia  +a  +PeUUy^))Y  Y  dy^  =0  (B-6) 

V  5,m  s,n  s^n  s,m/\-^  \  m  n/j  \  m  n  Yy  //  g,m  g,n  P  V  '^/ 


In  order  to  replace  the  first  term  in  the  equation  eq  B-6  with  an  integral  we  look  at  the 
ODE  for  the  structure  (eq.  A- 16)  and  write  it  out  for  two  different  eigenvalues. 


Y!'  +  =  0 


s^m  m  s,m 


(B-7  a) 


Y"  +a^Y  =0 

s,n  n  s,n 


(B-7b) 


Applying  similar  manipulations  to  the  above  two  equations  as  done  earlier  we  can  get 


Y"Y.. 


T 

,  +  j  fe  -  ««  )Yy,n.Ys,ndy^  ^  0 


(B-8) 


Again  the  first  term  may  be  integrated  by  parts  to  give 


1 

r 

=  _ (7'  7  - 7'  F  \  +[{a^ 

\  s,m^  s,n  J  \  m 

1 


^)y  Y  dy^ 

n  )  s,m  s,n  y 

-a^)Y  Y  dy^  =  0 

n  /  s,m  s,n  p 


(B-9) 


Where  we  have  made  use  of  the  outer  wall  boundary  condition  (eq.  A- 19)  and  simplified. 
Using  the  result  from  eq.  B-9  in  eq.  B-5  gives  rise  to  the  final  orthogonality  equation  as 

1  k  +  a.  +  P^V'iy'  (a,  +  a,  )Y,J„dy'  =  0  (B-10) 

0  1 
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Note  that  this  orthogonality  condition  is  different  from  a  typical  Sturm-Liouville  type 
orthogonality.  This  is  because  the  ODE’s  for  the  Y  equations  have  the  eigenvalues  mixed 
up  in  a  non-linear  fashion  and  does  not  follow  the  Sturm-Liouville  pattern.  It  is  useful  to 
also  note  that  the  form  of  eq.  B-10  suggests  what  we  are  attempting  in  these  problems. 
The  integration  first  occurs  in  the  gas  phase  in  the  domain  of  operation  of  the  gas 
temperature  solution  (between  0  and  1)  links  up  with  the  structure’s  solution  at  the 
interface  (y^  =1)  and  integrates  from  there  on  to  the  outer  wall.  We  also  see  the  presence 
of  the  ratio  of  wall  thermal  conductivities  (/f)  multiplying  the  structure’s  part.  Thus,  we 
can  see  how  the  thermal  coupling  phenomenon  is  expressing  itself  mathematically  in 
these  equations  and  leads  to  the  physical  implications  seen  in  the  results  section. 
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Appendix  C:  Evaluation  of  Coefficients 


Isothermal  Wall: 

Substituting  eq.  2-13  (see  Chapter  2)  into  the  jump  conditions  (eqs.  2-6  a  and  2-6  e) 
evaluated  at  the  flame  {x"=0)  leads  to: 


00 


00 


m=\  m=l 


(C-l) 


(C-2) 


The  derivation  for  the  elosed  form  expression  for  the  eoeffieients  starts  with  performing 
the  following  manipulations  in  order  to  isolate  Ap. 

l(eq.C-  2)  X  (-  PeU^(y^  )  +  )pjy^  (C-3) 

This  leads  to: 

1  00  ,  , 

m=\ 

ri/  ^  ^ 

"  I  Wndy^  +  Z  I  (y^  )  +  mmdy^ 

m=\ 

m^n 

The  seeond  term  on  the  right  hand  side  ean  easily  be  identified  as  being  part  of  the 
orthogonality  eondition  (for  the  pre-flame).  Therefore,  we  ean  use  eq.  B-5  to  rewrite  eq. 
C-4  as  follows: 


\'XB,X-PeU'(y')  +  Xyi,„wJy' 

m=\ 

=  A, i(-PeU'{y')yX,  )Hdf  + 1'  y  (-  W 

m=\ 

m^n 


(C-5) 
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Simplifying  the  sum  on  the  RHS  gives: 

1  CO  ,  , 

\X'^.X-P'>v'(y')+KV.¥Jy' 

m=\ 

\  /  \  1  ^ 

=  /’el'*  (y* ) + 2A.  V.V:f*  +  £  X  (-  W 

m=\ 

Now  we  use  the  seeond  jump  eondition  and  do  the  following  manipulation: 

J'l 

/eq.C-l)x(p^dy^ 

This  gives: 

,00  ,  ,  CO 

=  lotdy^  -!o 

m=l  m=\ 


(C-6) 


(C-7) 


(C-8) 


The  faet  that  we  have  this  naturally  oeeurring  “jump”  eondition  in  terms  of  the  derivative 
of  the  temperature  aeross  the  flame  enables  us  to  aeeommodate  the  non-linearity  of  the 
eigenvalue  in  the  ODE  for  Y.  Inserting  eq.  C-4  into  eq.  C-6  and  rearranging  gives: 

dn\l[-PeU\y^)  +  2A„  )f>ldy^ 

=  \'jndy^  +  iBJl(-PeU\y^)  +  A„-fi^)^„if/Jy^ 

m=l 


The  term  multiplying  An  on  the  left  hand  side  of  eq.  C-9  ean  be  evaluated  as  ean  the  first 
term  on  the  right  hand  side.  The  seeond  term  on  the  right  hand  side  looks  messy  as  it 
eouples  the  eigenfunetions  and  eigenvalues  from  the  pre  and  post  flame  regions  in  a 
eomplex  fashion.  However,  if  we  eonsider  the  ODEs  for  the  Y  funetion  in  the  pre  and 
post-flame  regions 

A'  +  C-^,-Pet'*(j’*)k=0  (C-10) 

r:+te+A/’el'*(y*)k=0  (C-11) 

Following  the  same  approaeh  used  to  derive  the  orthogonality  eondition,  we  eompute  the 
following 
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(C-12) 


I,  {(eq-  C-10)xi(/^-(eq.C-ll)x(p^  )dy^ 

This  leads  to: 

iifnWn,  -<l>n¥"n,)dy^  +{K  +  P m)\l[- PeU\y^)  +  - P ^ =0  (C-13) 
Again  because  of  the  boundary  conditions  (eq.  2-10  and  2-11),  the  first  integral  (after 
integration  by  parts)  can  be  easily  shown  to  be  zero.  This  means  the  second  integral  is 
also  zero  (as  both  X  and  fX  are  positive  numbers).  Using  this  result  in  eq.  C-9  we  obtain  An 
(eq.  2-14  in  Chapter  2).  Bn  is  obtained  by  similar  algebraic  manipulations. 

Conjugate  Heat  Transfer: 

The  “jump  conditions”  (eqs.  2-31  a,  b  and  2-32  a,  b  from  the  main  text)  and  the 
orthogonality  condition  developed  in  the  previous  section  are  used  to  get  analytical 
closed  form  expressions  for  the  coefficients  A„  and  Bn.  Substituting  eqs.  2-33  and  2-34 
into  eqs.  2-31  a,  b  and  2-32  a,  b  leads  to 


¥g,,n 

m=l  m=\ 

(C-14  a) 

CO  CO 

=  1 

m=l  m=\ 

(C-14  b) 

CO  CO 

=  I^Pm¥s,m 

m=\  m=\ 

(C-15  a) 

00  00 

ILA^X^  -h  ¥s,m  =  0 

m=l  m=\ 

(C-15  b) 

Startingwith  ^eq.C -  14a><{^ PeU(y^ )  +  X^^^dy^  gives: 
0 
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(C-16) 


1  1  00 

A,\(-PeU'{y*)  +  X,)pl,dy'  (-PeU'(y')  +  X,)p,,<P,,dy' 

0  0  “=1 

m^n 

1  CO  ,  . 

=  [-P^U(y^)  +  KVg,m^g,ndy^ 

0  "1=1 

Rewriting  the  second  term  on  the  left  hand  side  using  the  orthogonality  condition 
developed  in  the  previous  section  (eq.  B-10)  and  then  adding  the  missing  term  in  the 
summation  transforms  eq.  C-16  to: 


f  I  r  ' 

A\\(-PeU'{y')  +  2X,)plJf  +  2X,k\  i,dy' 

VO  1  ) 

1  00  ^  oo 

-  j  Z  ^s,ndy^  “  Z  )^sA,ndy^  (C'  1  7) 

Q  m=\  \  m=\ 

=  (-PeU{y^)  +  X„)//^^^(pgJy^ 

0  "1=1 


Once  we  have  eq.  C-17,  the  next  set  of  algebraic  steps  becomes  obvious:  We  need  to 
rewrite  the  terms  in  the  second  line  of  the  LHS.  From  a  cursory  inspection  of  the 

remaining  ‘'jump  conditions"  we  can  see  that  we  need  the  following  steps. 

/ 

First  je^.C- 746  gives 

0 


1  00  i  i  00 

\Yjd„,K^g,n,^g,ndy^  =\^g,ndy^  -\YAPn,Pm¥g,„,^g,ndy^  (C'lS  a) 


Then  ^eq.C  - 1 5bx  (l)^^dy^  gives 

I 


^00  ^00 

^s,ndy^  =  -  jZ^".^".  ¥s,Js,ndy^  (C-18  b) 

1  m=\  \  m=\ 


r 

and  ^eq.C -  15ax (j)^pdy^ gives 
/ 
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‘00  ‘00 
\Y.^'n^s,m^s,ndy^  =  j  Z 


(C-18  c) 


Using  eqs.  C-18  a-c  in  eq.  C-17  and  rearranging  gives 


1  00  ,  .  /  Z"  00 

=  \-P^Uiy^)  +  K-PmFg,n.<l>gJy^  \Y.^SPn.-K)¥sJsJy 


(C-19) 


Inspeeting  eq.  C-19  we  find  that  we  ean  evaluate  all  the  integrals  on  the  LHS  onee  we 
know  the  eigenvalues  and  eigenfunetions.  Therefore,  if  we  are  able  to  simplify  the  RHS 
to  known  quantities  we  will  be  able  to  express  the  eoeffieient  An  in  terms  of  known 
quantities.  Before  embarking  on  the  manipulations  it  is  useful  to  note  that  the  term  on  the 
RHS  has  a  form  similar  to  the  orthogonality  eondition  developed  in  the  previous  seetion 
exeept  that  in  the  orthogonality  eondition  we  were  either  in  the  pre-fiame  and  post-flame 
region  only  (as  deeided  by  the  minus  or  plus  in  the  equation)  whereas  here  the  pre-fiame 
and  post-flame  funetions  and  eigenvalues  are  mixed  together.  We  look  for  the 
simplifieation  of  the  term  on  the  RHS  by  doing  the  following:  Start  by  writing  the  ODEs 


for  the  gas  phase  iny^  for  the  post-flame  and  pre-flame  regions: 

+  =0 

C.+kE-l„/’et/*(/)V„=0 

/ 

Then,  take  j {eq. C-20ax (f)^^  -eq.C-20bx ^y^  to  get 


(C-20  a) 


(C-21  b) 


1  1 

I  (wlJ,,  -  €,v„W  K)PeU'  (V-  =  o  (c-22) 


195 


This  can  be  further  written,  after  integrating  the  first  term  by  parts,  applying  the  boundary 
eondition  at  the  eenterline  on  the  eigenfunetions,  and  using  the  gas-strueture  interfaee 
eonditions,  as 

1 

m  +K)\^K  -  P mV  P^U  ))l/  Q  (C-23) 

0 

Writing  out  the  structure’s  ODE  in  a  similar  fashion 

(C-24a) 

C+^X=0  (C-24b) 

r 

and  performing  similar  manipulations,  (i.e.  j  {eq.  C  -  24a  x(f)^^-eq.C- 24b  xy/^„  )dy^  , 

1 

integrating  first  term  by  parts  and  applying  outer  wall  boundary  eondition)  gives: 

r 

-  (KJs,n  -  +  j  ^  =  0  (C-25) 

1 

Together  eqs.  C-23  and  C-25  give 

+jii^n  - P m)- P^U \y^))i/ =  0  (C-26) 

1  0 

Using  eq.  C-26  in  eq.  C-19  and  rearranging  finally  gives  the  expression  for^f^  (eq.  2-37 
in  the  main  text).  Bn  ean  also  be  obtained  by  performing  similar  algebraie  manipulations. 
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Appendix  D:  Table  of  Maximum  Temperatures 


Pe 

k 

t 

NuE 

Theta  Max 

Flow  type 

0.1 

10 

2 

10 

0.4182 

Plug 

1 

10 

2 

10 

0.3938 

Plug 

5 

10 

2 

10 

0.1951 

Plug 

10 

10 

2 

10 

0.1018 

Plug 

5 

0.1 

2 

10 

0.2001 

Plug 

5 

1 

2 

10 

0.1999 

Plug 

5 

100 

2 

10 

0.1924 

Plug 

5 

10 

1.5 

10 

0.1926 

Plug 

5 

10 

3 

10 

0.1997 

Plug 

5 

10 

2 

0.01 

0.2696 

Plug 

5 

10 

2 

0.1 

0.2576 

Plug 

5 

10 

2 

1 

0.2181 

Plug 

5 

10 

2 

100 

0.1911 

Plug 

0.1 

10 

2 

10 

0.418 

HP 

1 

10 

2 

10 

0.3786 

HP 

5 

10 

2 

10 

0.1662 

HP 

10 

10 

2 

10 

0.1435 

HP 

5 

0.1 

2 

10 

0.264 

HP 

5 

1 

2 

10 

0.2097 

HP 

5 

100 

2 

10 

0.158 

HP 

5 

10 

1.5 

10 

0.1634 

HP 

5 

10 

3 

10 

0.1708 

HP 

5 

10 

2 

0.01 

0.2654 

HP 

5 

10 

2 

0.1 

0.2489 

HP 

5 

10 

2 

1 

0.1979 

HP 

5 

10 

2 

100 

0.1606 

HP 

Table  of  0max  for  heat  transfer  solutions  at  various  design  cases 
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